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ABSTRACT 

Magnetohydrodynamic  power  generation  has  now  been  under  active  de¬ 
velopment  for  over  four  years,  but  there  has  not  yet  appeared  any  complete 
description  of  the  theory.  This  report  is  intended  to  close  this  obvious  gap. 
Most  of  the  theory  presented  herein  was  developed  the  author  and  personnel 
at  the  General  Electric^  Company.  Some  of  this  material  has  not  been  published 
previously;  those  results  which  have  been  published  are  referenced. 

The  topics  covered  are:  electrical  conductivity  in  MHD  generators, 
optimum  "seed"  ratio,  local  analyses  of  the  continuous  and  segmented  electrode 
geometries;  Hall  geometry,  helical  flow  geometry;  magnetically  induced  ioniza- 
tion;  polytropic  efficiencies ;  compressible  analyses  of  the  constant  velocity, 
temperature,  Mach  number,  pressure  and  cross-sectional  area  flows;  end 
losses;  AC  generation;  cycle  efficiencies;  and  a  summary  of  generating  exper¬ 
iments  at  the  General  Electric  Company  and  other  places.  Geometrieti  other 
than  linear  are  not  considered  herein;  the  most  important  of  those  omitted  is  the 
vortex  generator. 

This  report  has  been  written  as  Chapter  14  for  a  forthcoming  book,  Engi¬ 
neering  Magnetohydrodynamics.  References  to  other  chapters  refer  to  chapters 
in  that  book. 
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14.1  INTRODUCTION 


In  principle,  the  magnetohydrodynamic  (MHD)  generator  is  a  variation  of 
the  Faraday  generator,  in  which  the  solid  conductors  are  replaced  by  a  moving 
electrically  conducting  fluid.  Historically,  Faraday  recognized  the  possibilities 
of  this  method;  he  planned  to  use  the  interaction  of  the  motion  of  the  Thames 
River  with  the  magnetic  field  of  the  earth  to  generate  electrical  power;  however, 
his  ejqjeriments  yielded  only  sporadic  results.  ^ 

A  more  obvious  source  of  electrical  power  is  fuel.  The  conversion  is 
usually  accomplished  by  combusting  a  chemical  or  fossil  fuel  with  an  oxidizer, 
which  converts  the  stored  chemical  ener;gy  of  the  fuel  into  thermal  energy  of  the 
products  of  combustion.  This  thermal  energy  is  then  usually  converted  directly 
into  mechanical  energy  in  an  internal  combustion  engine  or  gas  turbine ;  or  the 
thermal  energy  may  be  transferred  to  a  different  working  fluid  from  which 
mechanical  power  is  extracted  in  a  closed  cycle,  such  as  a  steam  or 

gas  cycle. 

In  the  case  of  nuclear  energy,  this  is  normally  converted  into  thermal 
energy  in  the  fuel  rods,  which  is  transferred  to  a  cooling  fluid,  and  then  converted 
into  mechanical  energy.  Finally,  the  mechanical  energy  is  converted  into  elec¬ 
trical  energy  by  means  of  an  alternator,  in  either  of  these  two  systems. 

Now,  in  both  cases,  if  the  primary  working  fluid,  that  is,  the  combustion 
gases  or  nuclear  reactor  coolant,  can  be  made  electrically  conducting  and  can 
have  part  of  its  thermal  energy  converted  into  directed  kinetic  energy,  then  the 
working  fluid  can  be  passed  directly  into  a  magnetic  field,  and  there  convert 
its  kinetic  energy  into  electrical  energy  directly,  thus  eliminating  the  "mechanical" 
intermediary  step. 

The  potential  advantage  of  MHD  power  generation  is  based  on  utilizing 
higher  temperatures  of  the  working  fluid  than  are  comp>atible  with  turbines.  Ihis 
may  lead  to  either  higher  overall  thermal  efficiencies,  or  in  the  case  of  "space 
power"  to  a  higher  reject  temperature  and  thus  a  reduction  of  the  radiator  area. 


1 


In  comparison  to  the  usually  wire-wound  generator,  the  MHD  generator  is 
more  complex  to  analyze,  because  the  electrical  conductivity  and  gas  velocity 
are  variable  properties  and  the  gas  is  not  constrained  to  move  in  an  exact  pre¬ 
scribed  path.  In  one  respect,  however,  the  MHD  generator  is  simpler  because 
the  magnetic  field  due  to  the  induced  currents  is  negligible. 

In  comparison  to  the  analysis  of  compressible  flow  in  a  duct,  the  flow  in  the 
MHD  generator  is  also  more  complex  because  the  Lorentz  force  and  ohmic  heat¬ 
ing  vary  over  a  cross-section,  which  changes  the  boundary  layer  growth. 

Also,  the  Hall  effect  can,  imder  certain  circumstances,  introduce  strong  cross- 
flow  components.  To  make  the  analysis  tractable,  we  will  therefore  make 
suitable  simplifying  assumptions  as  needed. 

Tbe  basic  principle  of  the  MHD  generator  is  that  the  electrical  power 
E  •  J  can  be  made  negative  by  the  presence  of  the  magnetic  field;  that  is,  since 
j  =  0  •  (E  +  V  X  B),  J  can  have  a  component  in  a  direction  opposite  to  E.  The 
total  generated  power  in  an  MHD  generator  is  the  integral  of  E  •  ^  over  the 
entire  volume  of  the  generator,  including  the  inlet  and  exhaust.  Since  E  =  -V<P  - 
( 3  A/  9t),  this  can  be  written  as : 


P=-  /E-JdV^  /j-V  <p  dy+ •  J  d^ 


or,  since  V*  (<p^  =  <pV’^  +  J*V(p, 


P=  /v-  (<pJ)d^  f  <pV-  Jd  .  J  d 


Finally,  the  first  integral  is  converted  to  a  surface  integral;  and  in  the  second 
integral  V*^  =  -3p  /9t,  where  o  is  the  excess  charge  density,  so  that 


p=  /0j.  ds+  fipjf  diA  jd^" 


(14.1) 
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In  an  MHD  generator  in  which  the  flow  is  steady,  the  last  two  terms  of 
(14. 1)  are  zero.  In  the  first  term,  the  normal  component  of  the  current  on  all 
insulator  surfaces  is  zero.  Thus,  the  power  appears  as  current  to  the  electrodes 
if  the  potential  on  the  different  electrodes  is  not  the  same.  The  potential  dif¬ 
ference  between  electrodes  is,  of  course,  the  ratio  of  the  electrode  current  to 
the  load  resistance. 

In  an  electrodeless  induction  type  generator,  the  first  term  of  (14. 1)  is 
zero.  Also,  it  is  common  to  design  it  in  such  a  way  that  the  potential  is  zero. 
Thus,  the  power  delivered  to  the  field  coils  comes  from  the  last  term. 

The  second  term  in  (14. 1)  corresponds  to  capacitive  coupling  with  the  flow 
and  corresponds  to  electrodeless  electrohydrodynamic  power  generation. 

14.  la  Electrical  Conductivity  in  MHD  Generators 

To  make  the  working  fluid  electrically  conducting,  there  are  three  general 
methods  as  follows; 

(i)  Non-Thermal  Ionization 

For  a  working  fluid  which  is  normally  unionized,  ionization  can  be 

achieved  by  electron  beams,  d.c.  orr.f.  discharges,  fission  product 

ionization,  etc.  However,  the  applicability  of  any  of  these  methods  for 

MHD  power  generation  has  yet  to  be  demonstrated.  Historically,  electron- 

beam  ionization  was  the  first  method  tried  with  combustion  gases,  but 
2 

was  not  successful.  However,  in  pure  or  noble  gases,  where  the  re¬ 
combustion  rate  for  electrons  is  smaller,  this  method  may  be  proved 
useful. 

In  combustion  flames,  there  is  always  a  small  amount  of  chemioniza- 
tion  as  a  result  of  the  chemical  reactions.  Although  this  can  be  enhanced 
somewhat  by  the  addition  of  small  amounts  of  metals,  generally  the  degree 
of  ionization  is  insufficient  to  be  useful  for  an  MHD  generator,  and  further¬ 
more,  there  are  usually  a  large  proportion  of  negative  ions  present,  so 
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that  the  number  of  free  electrons  present  is  very  small.  It,  therefore, 
appears  that  non-thermal  ionization  can  be  achieved  only  in  special  cases, 
see  Section  14. 3e. 

(ii)  Thermal  Ionization 

Thermal  ionization  refers  to  thermochemical  equilibrium  such  that 
the  thermal  energy  of  the  gas  is  an  appreciable  fraction  of  the  ionization 
potential  of  one  of  its  constituents;  then  Saha's  equation  predicts  that  a 
high  number  density  of  electrons  will  be  present.  Now,  in  ordinary  flames, 
which  burn  a  fuel  with  an  oxidizer,  the  flame  temperature  is  rarely  high 
enough  to  cause  appreciable  thermal  ionization  of  the  products  of  combustion, 
mainly  because  of  the  relatively  high  ionization  potential  of  the  products 
of  combustion.  Fortunately,  a  number  of  elements  have  lower  ionization 
potentials,  see  Table  14. 1.  If  a  small  fraction  of  alkali  metal  "seed"  is 
added  to  a  combustion  flame,  then  sufficient  ionization  may  be  obtained. 
Unfortunately,  in  flames  of  fossil  fuels  with  air,  the  flame  temperature 
is  too  low  to  ionize  even  the  "seed" ;  the  air  must  be  either  preheated  or 
enriched  with  oxygen  to  obtain  the  required  temperatures  of  4000  -5000°F.  ^ 

For  working  fluids  which  are  heated  by  solid  fuel  rods  in  a  nuclear 
reactor,  it  is  unlikely  that  sufficiently  high  temperatures  for  thermal  ioni¬ 
zation  can  be  achieved  in  the  near  future;  hence,  techniques  for  non-thermal 
ionization  will  have  to  be  developed  for  this  application.  This  problem 
will  be  alleviated  if  high -temperature  gaseous  nuclear  reactors  are  de¬ 
veloped.  The  application  of  MHD  generators  to  controlled  nuclear  fusion 
has  not  been  analyzed. 

(iii)  Liquid  Metal  Conductors 

Even  if  the  gas  is  non-conducting,  there  is  a  possibility  that  it  may 
be  mixed  mechanically  with  liquid  metals  as  a  foam  or  droplets.  If  the 
liquid  metal  is  foamed  with  the  gas,  then  the  entire  mixture  may  be  passed 
through  the  generator.  If  the  liquid  metal  is  mixed  with  the  working  gas  as 
droplets,  then  the  mixture  must  be  expanded  to  accelerate  the  drops;  that 
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is,  the  thermal  energy  of  the  gas  must  be  converted  into  directed  kinetic 
energy.  This  directed  kinetic  energy  must  then  be  exchanged  with  the 
liquid  metal  droplets.  Next,  the  kinetic  energy  of  the  liquid  metal  drop¬ 
lets  must  be  converted  into  a  pressure  head  in  a  seperator,  since  this 
pressure  head  must  be  used  to  flow  the  liquid  metal  throu^  the  generator 
against  the  Lorentz  force.  This  technique  has  not  yet  been  demonstrated 
experimentally . 

14.  lb  MHD  Generator  Geometries 

There  are  also  several  different  MHD  generator  geometries  presently 

under  consideration.  The  chief  requirement  for  any  geometry  is  that  there  be 

a  component  of  the  gas  velocity  which  is  not  parallel  to  the  magnetic  field,  so 

that  a  Faraday  electric  field  v  x  B  is  created.  This  field  will  then  cause  electric 

currents  and  fie^s  in  the  plane  normal  to  the  magnetic  field.  The  second 

requirement  is  that  there  be  electrodes  in  this  plane  which  collect  the  current 

xre, 

for  the  electrical  load.  The  simplest  geometries  for  accomplishing  this  the 

A 

linear  (together  with  its  many  variations),  vortex,  and  radial  outflow,  see 
Figure  14.1.  These  are  described  qualitatively  below,  and  the  linear  geometry 
is  examined  more  thoroughly  in  the  remaining  sections  of  this  chapter. 

(i)  Linear  MHD  Generator  -  "Duct  or  Chaiinel"  Geometry 

The  simplest  MHD  generator  geometry  is  the  "linear"  geometry  in 
which  the  gas  flows  through  a  linear  duct  or  channel,  see  Figure  14. 1. 

The  magnetic  field  is  at  right  angles  to  the  gas  flow  velocity,  which  induces 
a  Faraday  electric  field  at  right  angles  to  both  the  flow  velocity  and  magnetic 
field.  If  suitable  electrodes  are  placed  on  either  side  of  the  channel  and 
connected  through  an  electrical  load  or  resistance,  then  current  will  flow 
through  the  gas,  electrodes,  and  load. 

In  addition  to  the  Faraday  electric  field,  the  magnetic  field  causes 
a  Hall  current  to  flow  in  the  direction  of  fluid  flow;  that  is,  the  current 
does  not  flow  strai^t  across  the  duct.  This  phenomena  can  be  prevented 
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Figure  14.  1.  MHD  Generator  Geometries 
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by  employing  "segmented  electrodes"  which  is  a  variation  of  the  linear 
generator.  This  allows  a  Hall  electric  field  to  develop  in  the  flow  direction. 

In  still  another  variation  of  the  linear  generator,  the  Faraday  current  is 
shorted,  and  the  Hall  current  is  allowed  to  flow  throu^  the  electrical  load. 
Other  variations  are  also  possible,  see  Section  14.3. 

(ii)  Vortex  MHD  Generator 

In  addition  to  the  linear  geometry  with  its  many  variations,  there  is 
also  interest  in  the  vortex  or  spiral  geometry  in  which  the  gas  is  intro¬ 
duced  tangentially  into  a  cylindrical  geometry,  and  withdrawn  along  the 
surface  of  an  inner  coaxial  cylinder,  see  Figure  14.1.  The  magnetic 
field  is  in  the  axial  direction,  and  the  inner  and  outer  cylinder  are  the  two 
electrodes.  Under  open  circuit  conditions,  free -vortex  flow  should  re¬ 
sult;  that  is,  the  tangential  velocity  should  be  inversely  proportional  to 

the  radial  position.  Under  load,  the  tangential  is  poetically  constant  in 

A 

the  radial  direction.  When  the  inner  cylindrical  diameter  is  much  smaller 
than  the  outer  cylinder,  the  gas  makes  several  revolutions  in  the  generator; 
thus,  this  geometry  permits  a  longer  magnetic  interaction  length;  or 
alternately,  for  a  given  interaction  length,  the  vortex  generator  has  a 
more  compact  magnetic  field.  On  the  other  hand,  the  Hall  currents  tend 
to  flow  in  the  tangential  direction;  thus,  to  prevent  the  Hall  reduction  in 
electrical  conductivity,  the  gas  pressure  must  be  sufficiently  high  that 
03  T  <  1  at  the  exit.  (Note  that  when  the  inner  radius  is  almost  equal  to 

G 

the  outer  radius,  the  vortex  geometry  is  very  similar  to  the  linear  geometry. ) 
Up  imtil  the  present  time,  the  vortex  geometry  has  not  been  demonstrated 
experimentally,  and  will  not  be  considered  in  this  chapter. 

(iii)  Radial  Outflow  MHD  Generator 

A  variation  of  the  vortex  generator  is  one  in  which  the  gas  is  injected 
radially  outward  from  the  inner  cylinder.  In  this  case,  the  Faraday  current 
flows  tangentially  and  the  Hall  current  flows  radially;  the  latter  interacts 
with  the  magnetic  field  to  rotate  the  flow,  so  that  flow  becomes  a  spiral 
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outward.  This  variation  is  essentially  the  same  as  the  Hall  geometry, 
with  the  duct  bent  so  that  the  Lorentz  force  caused  by  the  Hall  current  is 
equal  to  the  centrifugal  force  in  the  fluid. 

(iv)  Other  Considerations 

All  of  the  geometries  considered  above  generate  dc  power.  However, 
by  switching  the  magnetic  field  or  by  using  a  traveling  magnetic  field,  it 
may  be  possible  to  generate  ac  power.  Although  this  has  not  yet  been 
demonstrated  experimentally,  the  general  principles  for  electrodeless  ac 
"induction"  generators  also  will  be  covered  in  this  chapter. 

In  order  that  the  ionized  fluid  behave  as  a  continuum  conductor  in 
an  MHD  generator,  the  Debye  length  must  be  much  smaller  than  the  size 
of  the  apparatus:  otherwise,  charge  separation  will  occur. 

To  prevent  ion  slip,  the  magnetic  field  must  not  be  so  large,  nor  the 
density  so  low  that  Tj  exceeds  unity.  This  will  be  explained  later. 

The  pressure,  conductivity,  gas  velocity,  and  magnetic  field  must  be 
such  to  make  the  interaction  len^^th  L,  defined  as 


L  = 


a  u  B 


(14.2) 


equal  to  the  length  of  the  path  of  the  gas  through  the  generator,  where  p 

2  ^ 

is  the  inlet  static  pressure,  and  o,  u,  and  B  are  average  values  in  the 
generator.  This  condition  insures  that  an  appreciable  fraction  of  the  inlet 
enthalpy  is  converted  into  electrical  energy. 

Finally,  if  electrodes  are  used,  the  current  density  from  the  cathode 
must  equal  or  exceed  the  current  density  in  the  gas;  otherwise,  the  effec¬ 
tive  electrical  resistivity  of  the  gas  will  be  increased. 
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14. 2  OPTIMUM  SEEDING  RATIO 


As  indicated  in  Section  14. 1,  to  make  a  combustion  gas  (or  working  fluid) 
electrically  conducting,  it  is  necessary  to  add  to  it  a  small  amount  of  seed 
material  which  has  a  low  ionization  potential,  see  Table  14.1.  There  are  two 
ways  in  which  this  may  be  done,  either  as  a  pure  element,  or  as  a  compound. 

To  determine  the  degree  of  ionization,  Saha's  equation  may  be  used;  thus,  it 
is  generally  necessary  to  determine  both  the  depletion  of  electrons  due  to  elec¬ 
tron  attachment  to  gas  atoms,  and  the  depletion  of  seed  atoms  due  to  chemical 
combination  with  the  other  atoms  present  in  the  gas.  This  generally  requires 
the  simultaneous  solution  of  Saha's  equation,  the  electron  attachment,  and  the 

A- 

chemical  equilibrium  equations.  This  may  be  represented  as  follows: 

S  <  ^  S  +  e 
e  +  A  <— 1 A 
S  +  B  ^Z»SB 
C  +  B^  D 
etc. 

where  S  is  a  seed  atom,  and  A,  B,  C,  D,  etc.,  are  gas  molecules. 

At  temperatures  of  about  3000°K,  all  reactions  except  the  first  generally 
can  be  neglected;  thus,  simplified  equations  can  be  obtained  for  the  optimum 
seed  concentration.  Generally,  the  upper  limit  for  the  seed  concentration  can 
be  obtained  by  two  different  sets  of  approximations,  as  follows: 

14. 2a  Isothermal  Optimum  Seeding 

The  temperature  of  the  mixture  is  assumed  to  be  unaffected  by  the  addition 
of  seed  material.  This  is  true  in  combustion  gases  if  the  seed  percentage  is 
very  small,  or  for  the  working  fluid  of  a  closed  cycle,  if  the  seed  plus  working 
fluid  is  heated  to  the  same  uniform  temperature.  From  Chapter  V,  the  scalar 
resistivity  of  the  gas  is  given  by: 
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TABLE  14.  1.  MOLECULAR  WEIGHTS  AND  IONIZATION  POTENTIALS 


(From  33rd  Ed. ,  Handbook  of  Chemistry  and  Physics,  Chemical  Rubber 


Publishing  Co. ) 


Gas 

Molecular  Weight 

Noble  Gases 

Helium 

4.03 

Neon 

21  .*83 

Argon 

39.4 

Krypton 

83.7 

Xenon 

130.2 

Common  Gases 

H 

1.008 

”2 

2.016 

N 

14.008 

28.016 

O 

16.000 

'^2 

32.000 

CO 

28.01 

40.02 

NO 

30.008 

Metal  Vapors  (Atoms) 

Lithium 

6.940 

Sodium 

23.00 

Aluminum 

26.97 

Potassium 

39.  10 

Calcium 

40.8 

Rubidium 

8.5.48 

Cesium 

132.91 

Barium 

1.37.36 

Mei’cury 

200.61 

Ionization  Potential  (electron-volts) 


I 

II 

24.46 

54.  14 

21.47 

40.9 

15.68 

27.76 

13.93 

26.4 

12.08 

21.  1 

13.53 

- 

15.6 

- 

14.48 

29.47 

15.51 

- 

13.55 

34.93 

15.51 

14.  1 

14.4 

9.5 

5.363 

75.26 

5.  12 

47.06 

5.96 

18.74 

4.318 

31.66 

6.09 

11.82 

4.  16 

27.36 

3.87 

14.  8 

5.  19 

9.95 

10.. 39 

18.65 

10 


(14.3) 


■n 


m  <c  > 
e 


2 

e 


n  Q 
ns  es 

n 


'el 


where  n  is  the  number  density  of  neutral  seed  atoms;  Q  is  their  average 
ns  6S 

cross  section;  n  is  the  number  density  of  the  rest  of  the  gas  atoms  and  Q  is 
g  ^  ^  eg 

their  average  cross  section;  tj^  ^  is  the  electron-ion  resistivity,  which  from 
Chapter  V  is : 


n 


e  I 


=  65.3  X  T 


-3/2  i 


1.24x10  T 


■1  ^3/2 


In  n 

e 


(14.4) 


From  Chapter  VI,  the  number  density  of  electrons  for  thermal  ionization 
of  a  single  species  is  given  by: 


, _  2^/4-  -T  //T 

n  V2G  (2TTm  k  T/h  )  e 

e  f  ns  e 


(14.5) 


where  Tj  =  e  E./k  and  E.  is  the  ionization  potential.  Substitution  of  (14.4,  14.5) 
into  (14.3)  yields : 


n  = 


m  <  c  > 
e  e 


e^  VTg 


T 

3/4 

ZT 

e 


ns  es 


1/2 

1 

ns 


65.3  T 


-3/2 


-  ^.n  n 
4  ns 


+  other  terms 


(14.6) 


where  the  other  terms  do  not  depend  on  n  . 

ns 


1  1 


To  obtain  the  optimum  seed  atom  concentration,  (14. 6)  must  be  differentiated 

with  respect  to  n  ,  and  the  resulting  expression  set  equal  to  zero.  TTie  result 
ns 


is: 


n 

ns 

n 

g 


iL±_iI 

(1  -  6) 


Q 


Q 


es 


(14.7) 


where 


e  =  '^ei/^'^el  ^  (14.8) 

and  where  cr  a  ,,  and  ^n  A  are  defined  in  Ch.  V. 
en  el 

Equations  (14.7,  14.8)  must  be  solved  iteratively,  by  first  assuming 

6  =  0,  then  obtaining  n^^  from  (14.7).  Then  n^  is  obtained  from  Saha's  equation 

and  6  is  calculated  from  (14.8).  This  value  is  then  substituted  back  into  (14.7) 

and  a  new  value  of  n  is  calculated,  etc.  Tlie  total  required  seed  n'  is  the  sum 
ns  ’  s 

of  the  final  value  of  n  and  n  . 

ns  e 

This  simplified  procedure  works  well  for  mixtures  of  pure  seed  and  pure 
noble  gases,  see  Figure  14.2,  where  the  conductivity  of  a  mixture  of  cesium 
and  argon  has  been  plotted  against  mixture  ratio.  For  small  mole  fractions  of 
cesium,  the  electrical  conductivity  increases  rapidly  with  increasing  cesium 
fraction.  Beyond  the  maximum,  the  conductivity  decreases  sli^tly  because  the 
cross-section  of  cesium  is  larger  than  that  of  argon.  Note,  however,  that  the 
value  of  decreases  monotonically,  because  depends  on  the  average 

cross-section,  which  increases  with  increasing  cesium  fraction.  Some  average 
values  of  the  cross  section  are  shown  in  Table  14. 2 

14.2b  Isoenergetic  Optip;um  Seeding 

Equations  (14. 7,  14. 8)  can  be  used  with  combustion  gases  only  if  the  seed 
atom  cross-section  is  much  larger  than  the  average  cross-section  for  the  combus¬ 
tion  gas;  for  otherwise.  (14.7  predicts  a  large  fraction  of  seed  atoms,  which  will 


a  (mho/m ) 


t 


Figure  14.  2. 


Electrical  Conductivity  and  Hall  Parameter  in  a  Mixture  of 


Cesium  and  Argon  at  1  Atm.  Pressure  and  2500‘^K. 

Q  =  3.5  X  10~^^cm^,  Q  =  2  x  10"^”^  cm^.  B  =  10,  000  gauss 
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TABLE  14.2.  AVERAGE  VALUES  OF  ELASTIC  ELECTRON  CROSS-SECTIONS 
WITH  VARIOUS  ATOMS  AND  MOLECULES 


Elastic 

Electron  Temperature 

Cross-Section, 

0 

CO 

0 

Particle 

in  10"^^  cm^ 

Noble  Gases^ 

A 

0.02 

2-3.5 

.  11 

7.7 

2.38 

50 

He 

0.54 

0-30 

Ne 

.  142 

0-8 

2 

Alkali  Metals 

K 

1.0 

2.4 

3.6 

1.1 

3 

Combustion  Gas  Constituents 

®2 

.31 

3-4 

N2 

.6 

3-4 

H2 

.85 

3-4 

CO 

.81 

3-4 

C02 

1.2 

3-4 

HgO 

3.6 

5 

N 

1.4 

3-4 

0 

1.7 

3-4 

1.  S.  Brown,  Basic  Data  of  Plasma  Physics,  Wiley,  1959. 

2.  G.  J.  Mullaney  and  N.  Dibelius,  ARS  Journal,  31,  1575  (1961);  G.  J. 
Mullaney,  P.H.  Kydd,  N.  Dibelius,  J.  Appl.  Physics,  32,  668  (1961). 

3.  Data  Compiled  by  W.  Chinitz,  L.  Eisen,  R.  Gross,  Am.  Rocket  Society 
preprint  706-58,  1958. 
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lower  the  flame  temperature  by  dilution  and  therefore  violates  the  assumption 
that  the  gas  temperature  remains  constant.  For  this  latter  case,  the  analysis 
must  include  the  effect  of  the  seed  atom  dilution  on  the  gas  temperature. 

To  analyze  this  case,  the  following  assumptions  are  made: 

i)  the  percentage  of  seed  is  small  and  the  degree  of  ionization  is  small; 

ii)  the  seed  atoms  and  electrons  do  not  react  chemically  with  the  rest  of 

the  gas.  Let  the  initial  enthalpy  of  the  gas  be  h  ;  and  the  initial 

og 

enthalpy  of  all  the  added  seed  be  h  ..  It  is  assumed  that  the  seed  is 

os* 

added  as  an  unionized  element.  Then  conservation  of  enthalpy  yields: 


nh  +nh=nh  +nh 
s'  s’  g  g  g  og  s  os' 


(14.9) 


If  the  seed  atoms  are  atomic,  then  the  enthalpy  per  seed  particle  is: 


h  =  I  k  T  +  ak  T, 
s'  2  I 

where  a  is  the  degree  of  ionization  defined  as : 


(14.10) 


n 

e 


n,  +  n 
1  s 


n 

e 


(14.11) 


and  where  electronic  exitation  has  been  neglected.  The  enthalpy  of  the  combustion 
gas  is  ~  (N  +  2)  k  T;  where  N  is  the  average  degrees  of  freedom  of  the  gas; 
thus  the  ratio  of  specific  heat5of  the  combustion  gas  is  y  =  (N  +  2)/N.  Finally, 
it  may  be  assumed  that  the  seed  material  has  zero  enthalpy  when  added  to  the 
flow.  The  temperature  of  the  gas,  after  addition  of  the  seed  material,  is  there¬ 
fore  given  by : 


T 


-Z-T 
y-1  0 


_)L 

y-1 


(14. 12) 
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where  0  =  n  /n  is  the  mole  fraction  of  added  seed.  Hie  upper  limit  for  0  can 
s*  g 

be  obtained  by  assuming  a.  «  1.  Then  for  0  «  1,  (14. 12)  reduces  to: 


(14.13) 


and 


^  ^  _  5.(y_-l)_  ^ 

3  0  2‘y  o 


(14.14) 


where  T  is  the  flame  temperature  without  seed.  Next,  (14.  6)  is  differentiated 

°  -3/4 

with  respect  to  0  using  (14. 14)^ and  equated  to  zero,  but  the  terms  T  and 

T)  ,  may  be  neglected  in  comparison  to  the  exponential  term.  The  result  is: 
el 


i  .-3/2 


+ 


+ 


LLlrJl 

2y 


T,  T 
1  o 


(14.15) 


Finally,  if  is  of  the  order  of  and  0  «  1,  then  may  be  neglected. 
Also,  after  the  seed  addition,  T  ~  T^  if  0  is  small.  Then,  solving  (14. 15)  for  0 
one  obtains: 


J>Yyri) 

y 


T 


o 


(14.16) 


Now,  in  combustion  gases  y  ~  1. 2;  thus,  the  final  expression  for  the 
maximum  seed  concentration  is 


T  k  T 


(14.17) 
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As  an  example,  consider  a  flame  temperature  of  3000°K  so  that  T  is 

o 

about  0. 25  volts;  the  ionization  potential  of  cesiiun  or  potassium  is  about  four 
volts;  thus,  the  mole  fraction  of  pure  seed  should  be  less  than  1/16  or  6%.  It 
therefore  appears  that  the  seed  fraction  for  combustion  gases  must  be  small; 
in  fact  up  to  the  present  time  concentrations  of  only  0. 1%  to  4%  have  been  con¬ 
sidered.  Of  course,  another  factor  in  keeping  the  seed  ratio  small  is  the  cost 
of  the  seed  material.  Also,  if  it  is  desired  to  minimize  the  Hall  parameter, 
seed  concentrations  greater  than  6%  could  be  used,  since  electron  cross-section 
of  ions  is  very  large. 
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14. 3  LOCAL  INVISCID  ANALYSIS  OF  LINEAR  MHD  GENERATORS 

The  local  analysis  of  an  MHD  generator  refers  to  the  performance  between 

two  adjacent  cross-sections  in  the  downstream  direction  in  a  linear  geometry. 

The  advantage  of  a  local  analysis  is  that  compressibility  and  area  changes 

generally  can  be  neglected,  if  the  two  cross-sections  are  sufficiently  close.  It 

is  also  usual  to  neglect  viscous  and  thermal  losses  to  the  walls,  in  comparison 

c 

to  the  Lorentz  force  and  electrical  energy  extration,  respectively.  This  assump- 

A 

tion  is  valid  if  the  viscous  Reynolds  number  of  the  flow,  based  on  the  hydraulic 
diameter,  is  very  large,  compared  to  the  magnetic  interaction  parameter.  The 
effect  of  these  losses  is  discussed  in  Section  14.4.  In  addition,  the  magnetic 
Reynolds  number  in  most  MHD  generators  will  be  small  so  that  the  induced 
magnetic  field  may  be  neglected.  With  these  two  assumptions  the  velocity  U 
is  taken  constant  in  the  x  direction  and  the  magnetic  field  in  the  z  direction  is 
taken  to  be  constant,  see  Figure  14.3a.  The  induced  electric  field  v  x  B  is  then 
UB  in  the  negative  y  direction.  Then  from  Chapter  V,  ohm’s  law,  including 
the  effects  of  ion  slip,  becomes : 

j  - - 2 - 

X  «  ,2  „  2 

(1  +  )  +8 
I  e  e 


v  2  2 

^  (I  +  B.B)  +b^ 

1  e  e 


(14. 18) 


[d  ^  e, 


8  )  (E  -  U  B)  +  8  E 

e  y  6  X 


2 

where  a  is  the  scalar  electrical  conductivity  neT/m,B.  =  a;T.  ,  and  6  = 

e  eel  Iln  e 

03  T  . 

®  ®  layeri 

In  the  absence  of  velocity  or  thermal  boundary  both  E  and  E  may  be 

'A  X  y  •' 

considered  constant  (although  more  generally  they  may  be  a  function  of  x).  Then 
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14.  3a. 


Schematic  Diagram  for  "Local  Analysis" 
of  Linear  MHD  Generators 
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all  variables  are  constant  across  a  cross-section,  so  that  the  electric  field  and 
current  in  the  direction  of  the  magnetic  field  is  zero.  * 

We  will  next  use  this  simple  model  to  analyze  the  three  most  common 
t3rpes  of  linear  generators:  continuous  electrodes,  segmented  electrodes,  and 
the  Hall  generator,  see  Figure  14. 3b^  and  a  variation  of  the  segmented  electrode 
geometry,  in  which  the  flow  is  helical  through  an  annulus. 

14.3a  Continuous  Electrodes 

With  the  electrodes  continuous  along  each  side  of  the  generator  with  a 
different  constant  potential  on  each,  no  axial  electric  field  can  develop  and 
therefore  =  0.  Then  from  (14. 18),  the  transverse  component  of  the  current 
is; 


^"1 - T  [(l-e^Bj)(E  -UB) 

e.i  +  B  •-  >  J 


(1+6  6J  +  B 
el  e 


(14. 19) 


The  open  circuit  condition  implies  that  j  =  0;  thus,  from  (14. 19),  the 
open  circuit  electric  field  is  E  (open  circuit)  =  UB.  Under  short  circuit  condi¬ 
tions,  E  =  0  (except  for  the  finite  electrical  resistance  of  the  electrodes  and 
y 

external  connection).  Thus,  in  general,  0  <E  <  UB.  It  is  common  to  express 

y 

the  electric  field  under  load  as  a  fraction  K  of  the  open  circuit  electric  field; 
thus,  for  the  continuous  electrode  geometry: 

E 

K  =  (14.20) 


where  for  generation  of  electrical  power,  0  <  K  <  1.  Note  that  j  is  then 
negative,  and  the  Lorentz  force  j  x  B  is  in  the  negative  x  direction;  that  is. 


♦Velocity  boundary  layers  on  the  insulators  cause  the  induced  field  v  x  B  within 
the  boundary  layer' to  be  less  than  that  at  the  centerline.  Under  these  conditions, 
the  current  in  the  boundary  layer  may  be  in  the  opposite  direction  as  that  in  the 
main  flow  and  thus  constitute  a  "leakage"  current.  However,  if  the  electrical 
conductivity  is  constant,  then  the  average  axial  velocity  across  a  cross-section 
may  be  used  instead  of  U. 
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CONTINUOUS  ELECTRODES 


HALL  GENERATOR 


Figure  14.  3b.  Linear  Generator  Geometries  Having 
Rectangular  Cross  Sections 
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opposite  to  the  flow  direction.  Thus,  in  the  generator,  the  Lorentz  force  tends 
to  retard  the  flow.  (If,  on  the  other  hand,  the  electrodes  were  attached  to  a 


voltage  source  so  that  E  >  UB,  that  is,  K  >  1,  then  j  is  positive  and  the 

y  y 

Lorentz  force  tends  to  accelerate  the  flow.  If  the  polarity  of  the  inaposed  voltage 
is  reversed  so  that  <  0,  the  Lorentz  force  will  retard  the  flow;  the  energy 
added  to  the  flow  from  the  voltage  source  liien  increases  the  thermal  energy  of 
the  gas.) 

In  addition  to  the  transverse  component  of  the  electric  current,  there  is 
also  a  component  in  the  downstream  direction  given  by: 


a  6^  (1-K)  UB 


(14.21) 


where  use  has  been  made  of  (14.20).  Thus,  the  current,  as  it  flows  across  the 

generator,  also  tends  to  flow  downstream.  Actually,  it  is  the  current  in  the 

downstream  direction  (14. 21)  which  is  responsible  for  the  Lorentz  force  on  the 

gas.  This  may  be  seen  by  simplifying  the  problem  somewhat:  let  oJt  T,  =  8.  =  0, 

I  in  1 

so  that  the  average  ion  motion  is  exactly  equal  to  the  gas  velocity  The 
current  given  by  (14. 21)  is  the  total  downstream  current;  if  the  average  ion 
velocity  is  equal  to  the  gas  velocity,  then  with  Bj  =  0  (14.21)  is  also  the  expres¬ 
sion  for  the  electron  conduction  current  relative  to  the  gas  motion,  since  the 
gas  is  electrically  neutral.  Thus,  the  diffusion  velocity  of  the  electrons  is  in 
the  downstream  direction  relative  to  the  gas,  that  is: 


Te^  =  „  u.^t^(1-K)UB 

('  ■*•  fe) 


(14.22) 


The  drj^  force  exerted  on  the  gas  by  the  electrons  is 
with  the  use  of  (14.22)  this  becomes: 


F  =  n  m 
X  e 


F 

X 


-g  (1-K)  UB 
e 


(14.23) 
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This  is  identical  to  the  Lorentz  force  jy®>  calculated  from  (14.19)  when  Bj  is 
neglected. 

(i)  Power  Density,  Continuous  Electrodes 

The  power  generated  per  unit  volume  is  P  =  -  E  •  j  which  becomes 


just  E^  since  E^  ~  0.  From  (14. 19),  (14.  20),  this  becomes: 


P  = 


aK  (1-K)  (1  +  6  Bj) 

2  2 


(14.  24) 


Equation  14.  24  can  also  be  written  as : 


P  = 


n  m  U  K  (1-K) 
6  6 


2  2 
9  (1+8  —  ) 

e  e  u 

_  e 

u  2 

2  2 
(1+0  — )  +  6 

eu  e 

e 


(14.25) 


where  U,,  M  are  the  electron  and  ion  mobilities  e  T,  /m„  e  t  /m  , 

1  e  In  1  I  e  e 

respectively.  From  (14.  25)  it  is  also  obvious  that  0  <  K  <  1  for  power 
generation.  The  coefficient  of  (14. 25)^n^  m^  v^bas  a  simple  interpre¬ 
tation:  consider  that  after  each  collision  with  a  heavy  particle  the  electron 

acquires  the  gas  velocity  U  in  the  downstream  direction.  But  then  the 
+h«! 

magnetic  field  turns  electron  velocity  in  the  +y  direction  parallel  the 

A 

electric  field,  so  that  its  kinetic  energy  is  transformed  into  potential 

energy.  If  the  electron  is  then  hit  again  by  a  heavy  particle,  it  again 

acquires  a  velocity  U  in  the  downstream  direction.  The  rate  at  which  an 

electron  acquires  kinetic  energy  from  the  gas  is  thus  approximately 

(  ^  m  V  \  the  total  rate  at  which  all  electrons  acquire  energy  from 
\2  e  e/  ^ 

the  gas  is  therefore  ~  r  n  m  V/T  . 

2  e  e  /  e 


This  same  result  can  be  obtained  from  a  different  point  of  view : 
between  collisions,  the  magnetic  field  turns  an  electron  in  the  y-direction 


a  distance  of  order  of  magnitude  equal  to  the  larmor  radius,  which  may 


be  taken  as  approximately  thereby  acquiring  a  potential  energy 

due  to  the  stationary  electric  field  equal  to  EU/oJ^.  The  total  rate  at 
which  the  electrons  acquire  potential  energy  is  then^^  EU/oJ^  Since 


E  Oi  UB,  and  o)  =  eB/m  ,  this  rate  is  again  equal  to 
e  e 

agrees  with  the  previous  argument. 


which 


The  factor  on  the  ri^t  side  of  (14. 25)  represents  the  effects  of  the 

Hall  reduction  in  electrical  conductivity  and  ion  slip  for  this  geometry, 

and  has  two  asymptotes.  The  first  asymptote  is  obtained  by  taking  or 

equivalently  or  9^  equal  to  zero,  which  is  equivalent  to  staling  that  the 

2 , 

ion  velocity  is  equal  to  the  gas  velocity.  This  factor  then  becomes  9  / 

2  ® 
(1+8^  ).  Thus,  for  a  given  value  of  r^,  the  power  P  initially  increases 

quadratically  with  the  magnetic  field,  since  B  =  n  B.  This  is  shown  in 

0  0 

Figure  14.4.  As  the  magnetic  field  is  increased  further,  the  first 

2  2 

asymptote  is  reached  8  /(I  +  6  ) - >1.  Thus,  the  maximum  energy 

®  ®  /  2 
transfer  or  power  generation  due  to  the  electrons  is(n  m  U  v 

\  e  e 

of  course,  by  the  loading  factor  term  K  (1-K).  In  Figure  14,4  the  inflec 


modified. 


tion  point  near  unity  represents  this  first  asymptote.  Since  magnets  gen¬ 
erally  consume  power  and  are  expensive  or  heavy,  it  is  desirable  to 

operate  at  a  point  somewhat  below  this  point;  suppose  the  power  level  is 
2 

set  at  P  T  /n  m  U  K  (1-K)  equal  to  80%  of  the  first  asymptote,  then 
2  ®  2  ^  ^ 

9  /(I  +  8  )  =  0.  8,  so  that  8  =  =  2.24.  Thus,  a  magnetic  field  which 

0  0  0 

corresponds  to  a  value  of  the  Hall  parameter  greater  Hun  two  will  not  lead  to 
substantially  greater  power  generation  by  electrons. 


(ii)  Magnetic  Interaction  Length 


With  =  0,  the  interaction  length  can  be  written  in  terms  of  the 
electron  mean  free  path  as  follows : 

2 


L  = 


-El. 


1  +  8 


aUB 


8 


Vm 

-S. 

m 


e 


(14.26) 


where  9?^is  the  gas  Mach  number. 
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Again,  as  9^  is  increased  from  zero,  the  interaction  length  decreases 

and  then  becomes  asymptotic.  Since  the  degree  of  ionization  will  be  of 

-2  -4  -2 

the  order  of  10  to  10  ,  the  square  root  of  the  mass  ratio  is  about  10  , 


and  the  Mach  number  is  of  order  of  magnitude  unity. 


the  interaction 


4  6 

length  of  the  generator  is  10  to  10  electron  mean-free  paths,  multiplied 
2  ,  2 

by  (1  +  6  )/6  For  constant  electron  mean -free  path,  increasii^  the 

G  6 

magnetic  field  decreases  the  interaction  length,  but  values  of  6^  beyond 

about  two  will  not  lead  to  any  appreciable  further  reduction.  On  the  other 

hand,  with  the  magnetic  field  constant,  if  6^  is  increased  by  decreasii^ 

the  gas  pressure,  then  the  electron  mean-free  path  will  increase,  which 

increases  the  interaction  length.  This  may  be  seen  by  substitution  of 

X  =m  <C  >8/eB  into  (14. 26).  The  interaction  length  then  varies  as 
e  2®  ®  ® 

(l  +  6  /8  ,  which  has  a  minimum  at  8  =1. 

0  0  0 

With  8  non-zero,  (14.25)  predicts  a  second  asymptote  when  8  >  H  /fi. 

1  0  0  1 

for  the  right  hand  factor,  equal  to  the  mobility  ratios  M^/Mj  which  is 

generally  much  greater  than  unity.  This  is  shown  on  the  right  side  of 

Figure  14.4  for  very  large  values  of  co  t  .  The  power  extraction  is  then 

2  ®  ® 

of  the  order  of  P  ~  n^  m^  U  /t  which  shows  that  for  very  large  values 
of  cOg  T^,  the  positive  ions  are  responsible  for  the  power  generation.  Be¬ 
cause  the  mobilities  vary  roughly  as  the  inverse  square  root  of  the  mass 

ratio,  the  asymptotic  interaction  length  for  8  large  becomes  L  =  n  X  /n  ‘’bn  . 

g  e  e  e  " 

However,  values  of  8^  larger  than  10  cannot  be  obtained  except  by  reducing 

the  pressure,  which  again  increases  the  interaction  length  because  of  the 

dependence  on  the  mean-free  path.  For  example,  consider  (i  /u.  =  200? 

4  ®  ^ 

from  Figure  14.4,  values  of  w  t  greater  than  2  x  10  are  required  to 

e  e 

approach  the  asymptote,  while  the  first  asymptote  is  reached  when  co  t  >  2. 

4  e  e 

Thus,  if  one  lowers  the  density  to  obtain  2  x  10  ,  then  the  electron  mean- 

4 

free  path  has  been  increased  by  a  factor  of  10  ,  so  that  the  interaction 
length  is  increased  by  two  orders  of  magnitude.  It,  therefore,  appears 
imlikely  that  an  MHD  generator  can  be  designed  to  utilize  the  ion  current 
with  continuous  electrodes. 
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(iii)  Efficiency,  Continuous  Electrodes 

From  (14. 24),  it  is  obvious  that  the  power  generation  density  is  a 
maximum  when  K  =  ^  which  obviously  corresponds  to  matched  impedance? 
that  is,  when  the  load  resistance  is  equal  to  the  generator  internal  resistance. 
However,  for  some  applications,  it  may  not  be  desirable  to  operate  at 
maximum  power  density,  but  instead,  closer  to  maximum  efficiency.  For 
..  the  purpose  of  a  local  analysis,  a  local  conversion  efficiency  may  be  de¬ 
fined  as  the  ratio  of  the  generated  power  to  the  flow  work  which  is  required 
to  overcome  the  Lorentz  force,  or: 


(14.27) 


In  the  present  geometry,  the  denominator  of  (14.27)  is  U  B,  while  the 
numerator  is  E^  thus,  the  local  efficiency  for  this  geometry  is  E^/UB  =K, 
from  (14. 20). 

The  local  efficiency,  as  defined  by  (14. 27),  is  also  equal  to  the  local 
adiabatic  or  polytropic  efficiency  (also  called  the  "stage"  efficiency)  when 
(i-K)  is  small  compared  to  unity,  see  Section  14.4.  Thus,  if  a 
high  polytropic  efficiency  is  required,  it  is  necessary  to  operate  at  values 
of  K  near  unity. 

Since  the  load  voltage  is  KUBh,  and  the  load  current  per  unit  length  of 
generator  is  (1-K)  UBb  where  is  the  "effective"  electrical  conductivity; 
the  external  load  resistance  per  unit  length  of  generator  is,  therefore. 


R  _Y _ Kh_  _  _K_  R 

L  J  (1-K)  c  b  1-K  i 
e 


(14.28) 


where  R.  is  the  effective  internal  resistance  of  the  generator,  h/a^bj 
Solving  (14.28)  for  the  value  of  K,  one  obtains: 


K 


1 


1  +  R./R 
1 


L 


(14.29) 
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which,  for  this  geometry  is  also  equal  to  the  conversion  efficiency.  One 
therefore  obtains  the  familiar  result  that  the  conversion  efficiency  in¬ 
creases  as  the  ratio  of  external  to  internal  load  increases. 


(iv)  Voltage  Control,  Continuous  Electrodes 


In  general,  the  continuous  electrode  geometry  will  be  operated  such 
that  K  1;  since  the  power  generated  is  a  parabolic  function  of  K, 
this  means  that  operation  will  be  on  the  falling  side,  see  Figure  14.5, 
This  means  that  as  the  generator  is  unloaded  (load  resistance  increases) 
the  power  to  the  remaining  load  decreases  although  the  load  voltage  in¬ 
creases  slightly.  If  the  magnetic  field  is  kept  constant,  then  the  effective 
electrical  conductivity  and  therefore  R.  will  remain  constant,  so  that  the 
local  conversion  efficiency  is  still  given  by  (14. 29). 


If  it  is  desired  to  keep  the  load  voltage  constant  as  the  load  is  changed, 
this  may  be  accomplished  by  varying  the  magnetic  field,  but  this  will 
change  the  internal  resistance.  For  example,  if  the  load  resistance  is 
increased,  the  load  voltage  will  increase;  to  keep  the  load  voltage  constant, 
the  magnetic  field  is  then  decreased.  However,  this  reduces  the  Hall 
parameter  P  ^  and  therefore  increases  the  effective  electrical  conductivity, 
thereby  simultaneously  decreasing  the  internal  resistance.  One,  therefore, 
expects  that  the  local  efficiency  will  increase  more  than  if  the  magnetic 
field  were  held  constant.  It  is  easily  shown  that  the  load  parameter  and, 
hence,  efficiency  is  given  by: 


K 


Z 


n  2  2 

-  4  (1  +  Z)  f  K 
o  o 

T(TTz) 


(14.30) 


where  ‘  is  the  initial  value  of  P  ;  8  was  assumed  to  be  zero,  K  is  the 
o  el  o 

initial  value  of  the  loading  factor,  and  Z  is  defined  as: 


Z 


K 

o 

1  -  K 

o 


(1  +  8^  ) 


(14.31) 
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GENERATED  POWER  DENSITY 


Figure  14,  5.  Local  Power  Density  For  Several  Generator 
Geometries  as  a  Function  of  the  Ratio  of 
Load  Voltage  to  Open  Circuit  Voltage 
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where  R,  is  the  new  load  resistance,  and  Rt  was  the  initial  value  of  the 
L  ■‘^o 

load.  For  example,  consider  K  =0.6,  Rt  /Ri  =  1.5,  6  =2.  If  the 

o  -^o  o  o 

load  resistance  is  doubled,  while  the  magnetic  field  is  decreased  to  keep 
the  load  voltage  constant,  then  K  becomes  0. 83.  On  the  other  hand,  if 
the  magnetic  field  were  kept  constant,  doubling  the  load  would  only  increase 
the  value  of  K  to  0 . 75 . 

In  addition  to  the  Lorentz  force,  which  acts  to  retard  the  flow,  there 
is  a  component  of  the  Lorentz  force  in  the  direction  of  the  cathode  which 
is  the  product  of  ®g/(l  +  ^j)  and  the  retarding  Lorentz  force.  This  is 

balanced  by  a  transverse  pressure  gradient;  if  the  magnetic  interaction 
parameter,  based  on  the  channel  width  is  large,  this  transverse  pressure 
gradient  codld  become  quite  large,  and  would  alter  the  density  or  tempera¬ 
ture  across  a  cross-section. 


14.3b  Segmented  Electrodes 

In  the  previous  section,  it  was  sjiown  that  with  increasing  magnetic  field, 
the  generated  power  density  due  to  electrons  becomes  asymptotic  to|n 


m  U 
e  e 


K  ( 1-K)  because  of  the  Hall  reduction  in  electrical  conductivity  in  the  direc¬ 
tion  of  the  electric  field.  To  remedy  this  situation,  the  segmented  electrode 

5  6 

generator  has  been  devised  ’  in  which  each  opposite  pair  of  electrodes  is 
connected  to  a  single  load,  see  Figure  14. 3.  In  this  arrangement,  the  current 

flow  is  essentially  transverse  to  the  flow,  and  no  net  Hall  current  exists  in  the 

A 

downstream  direction.  Since  the  gas  is  electrically  neutral,  this  means,  of 
course,  that  the  electron  and  ion  velocities  in  the  downstream  direction  are  equal 
to  each  other  and  are  also  equal  to  the  gas  velocity  in  the  downstream  direction, 
in  the  absence  of  ion  slip. 


With  the  Hall  current  equal  to  zero,  (14. 18)  predicts  that  an  electric 
field  will  develop  in  the  axial  direction  given  by 

B 

irre  <Ey-uB) 

el  ■’ 


E  = 

X 


(14.32) 
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so  that  the  current  in  the  y  direction  becomes; 


1  +  3 


(E  -  UB) 


e  I 


(14.33) 


As  in  the  pT'evious  section,  the  open  circuit  voltage  corresponding  to  =  0 
is  UB;  we  therefore  take  the  electric  field  under  load  as  KUB.  The  Lorentz 
force  i  X  B  now  acts  only  in  the  negative  x  direction  so  that  it  acts  only  to  re- 
tard  the  flow.  This  force  is  exerted  purely  on  the  ions  and  not  the  electrons, 
as  shown  in  the  following  argument.  The  electron  velocity  with  respect  to  the 
generator  may  be  resolved  into  two  parts,  one  due  to  the  electron  drift  motion 

E  X  B  and  the  other  in  the  downstream  direction  due  to  collisions  with  the  gas,  v 
~  c 


V 

~e  ^2 


+  X  V 

~  e 

c 


(14.34) 


see  Figure  14.6.  In  the  y-direction,  (14.34)  becomes  with  the  use  of  (14.32): 

V  =  -  E  b"^  =  (1-K)  CO  T  U/(l  +  3  5  )  (14.35) 

e  X  e  e  e  1 

y 

in  agreement  with  (14.33).  The  electron  velocity  in  the  downstreamdirection  is  equal 

ouiiih 

to  the  ion  velocity,  for  6  8  «  1  is  the  gas  velocity*,  thus,  from  (14.34); 

A  el 

V  =  E  B  +  V 
e  y  e 


or 


U 


KU  +  V 

e 


c 


(14.36) 
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ANODES 


CATHODES 


Figure  14.  6.  Electron  Motion  and  Electric  Fields  in  a 
Segmented  Electrode  Geometry 
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so  that  the  velocity  increment  in  the  downstream  direction  is  v  =  ( 1-K)  U . 

0 

Q 

Now,  the  total  force  on  electrons  in  the  x  direction  is  given  by: 


i’  =  -  ne  (E  +  V  B) 
e  X  e 

X  y 


(14.37) 


Substitution  of  E  from  (14.32)  (with  3  3  neglected)  and  v  from  (14.35)  shows 

X  0  1  6y 

that  the  net  axial  force  on  the  electrons  is  zero.  On  the  other  hand,  the  force 
on  ions*  in  the  x  direction  is  Fj  =  n  eE  ;  use  of  (14. 32)  then  yields  Fj  = 
ne  3  (K-1)  UB/(1  +  3  6  )  which  is  identical  to  j  x  B  as  obtained  from  (14.33). 

0  0  1  /W 

(i)  Power  Density  and  Efficiency,  Segmented  Electrodes 


The  power  density  is  -  E  •  j  =  -  E  j  ,  or: 

~  y  y 


K  (1-K)  g  U^B^ 

1  +  33 

e  I 

n  m  K(l-K)  3  ^ 

0  0  0 

’'e  1  +  3  ^  ti  M 

e  I  e 


(14.38) 


As  in  the  case  of  continuous  electrodes,  the  power  density  initially  increases 

quadratically  with  increasing  magnetic  field  (since  =  e  B/m^);  and  the 

last  factor  of  (14. 38)  becomes  asymptotic  to  M  /M  ,  see  Figure  14.4. 

0  1 

However,  for  intermediate  values  of  3^,  no  first  asymptote  appears  in  the 

power  density  as  was  the  case  for  continuous  electrodes.  Thus,  for 

1  3^  <  (#^g/Mj),  the  generated  power  of  the  segmented  electrode  geometry 

greatly  exceeds  that  for  the  continuous  electrode  geometry.  The  generated 

power  comes  from  the  ions  for  this  geometry;  in  fact,  the  electrons  actually 

*It  has  been  assumed  that  all  ions  are  positively  charged.  If  some  have  negative 
charges,  the  force  is  e  (nj^  -  nj")  Ej^.  But  nj'*'  -  nj“  =  Ug  in  an  electrically  neutral 
gas,  therefore,  this  same  expression  applies. 
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decrease  the  power  generation.  The  power  is  obtained  from  the  ions  by 
their  motion  against  the  x-component  of  the  electric  field; 


Pj  =  -  n^  e  U  e  (1-K)  U^B/(1  +  3^  3^) 


(14.39) 


The  power  generated  by  the  electrons  is; 


P  =  n  e  V  .  E 

2 

Now,  from  (14. 34)  the  electron  velocity  consists  of  two  parts;  E  x  B/B 

2  ~  ~ 

and  i  Vg  .  The  drift  velocity  E  x  B/B  generates  zero  power  since  E  x 

B  •  E  =  0;  while  the  part  due  to  v  in  the  x-direction  jdelds 

c 

=  -n^  e  (1  -  K)^U^  B/(l  +  6^3j)  (14.40) 


so  that  the  electrons  actually  consume  power;  that  is,  power  must  be 
expended  in  order  for  the  electrons  to  move  with  the  gas  velocity  U.  The 
sum  of  (14. 39),  (14.40)  is  identical  to  that  given  by  (14.  38). 


The  local  efficiency  is  still  equal  to  K,  since  =  E  j  /UBj  .  Thus, 

■*'  y  y  y 

the  increase  in  power  density  has  been  obtained  with  no  sacrifice  in  effi¬ 
ciency.  Since  the  power  density  is  still  a  maximum  when  K  =  -^  ,  this 
geometry  will  also  be  operated  on  the  downward  side  of  the  power  curve, 
see  Figure  14.5. 


14. 3c  Hall  Generator 

When  the  electrodes  are  segmented  it  was  shown  in  the  previous  section 
that  an  electric  field  develops  in  the  flow  direction,  see  (14.32).  This  electric 
field  is  maximum  when  the  opposite  electrodes  are  short-circuited;  that  is,  when 
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E  =  0 .  If  the  electrode  pair  at  the  inlet  is  connected  through  a  load  to  the  pair 
at  the  exit,  a  current  should  then  flow  through  the  load.  With  =  0,  the 
current  density  in  the  +x  direction  is  obtained  from  (14, 18); 


Jx  ' 


(1  +  3  E  +  3  UB 
e  I  X  e 


(14.41) 


When  open-circuited,  the  axial  electric  field  is  -  3^UB/(1  +  under  load  the 

axial  electric  field  will  be  some  fraction  K  of  this  value,  E  =  -  3  KUB/(1  +  3  3  ), 

H  X  e  H  e  I 

so  that;  ^ 

a3^(l-KjUB 


Jx  = 


(14.42) 


Note  that  the  definition  of  is  different  than  K  for  the  continuous  and  segmented 
electrode  geometry. 

(i)  Power  Density  Hall  Generator 

The  generated  power  density  is; 


j 


X 


a  2  2  2 

K  (1-K)  a  3  u  B 

M  H  6 


n  m  U^K(l-K) 
e  e  H  rt 

T 

e 


(14.43) 


In  this  geometry,  as  the  magnetic  field  is  increased  from  zero,  the  power 
increases  as  the  fourth  power  of  the  magnetic  field,  in  contrast  to  the 
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quadratic  dependence  of  the  continuous  and  segmented  electrode  geometries. 

3  -2 

An  as’rmptote  is  reached  when  3  ^  M  /M,,  equal  to  (M  /M,)  3  ,  which 

e  el  e  I  e 

is  considerably  smaller  than  the  asymptote  of  obtained  in  the  con¬ 
tinuous  and  segmented  geometry  generators.  However,  for  2  3 

/  1/2  ® 

(M  /M^)  ,  the  power  density  approaches  that  of  the  segmented  electrode 

0  X 

geometry. 

The  retarding  Lorentz  force,  j  B  becomes; 


F 

X 


_T  T-.2 

_  -  g  UB 

“1  +  3  3^ 

e  I 


1  +  3’  ^  K 

e  H 


1  +  3 


,  2 


(14.44) 


where  3^  is  an  effective  Hall  parameter  for  electrons  defined  as: 


3» 

e 


1  +  33 

e  I 


(14.45) 


Note  that  the  retarding  Lorentz  force  is  a  minimum  when  the  load  is  short- 
circuited;  that  is,  when  K  =  0.  This  is  in  sharp  contrast  to  the  previous 
geometry,  where  the  retarding  force  is  a  maximum  when  short-circuited. 
This  effect  is  again  due  to  the  force  on  the  ions;  in  the  segmented  electrode 
generator,  the  axial  electric  field  is  a  maximum  when  short-circuited;  in 
the  Hall  geometry  when  K  =  0,  there  is  no  axial  electric  field  and  there¬ 
fore  no  net  axial  force  on  the  ions.  This  may  be  seen  as  follows;  the 
axial  force  on  the  ions  is  (neglecting  ion  slip); 


F,  =  ilcE  =  -  n  e  3  UBK 
lx  1  X  e  e  H 


The  electron  current  relative  to  the  gas  is; 


j 


ex 


+  3  UB 

e 


(14.46) 


(14.47) 


3b 


so  that  the  electron  velocity  relative  to  the  gas  is 


-  e  T  B  UB 

e  e 


(1  +  3  )  m 


(1-K) 

H 


(14.48) 


The  motion  of  the  electrons  relative  to  the  gas  exert  a  friction  drag  on  the 
gas  given  by; 


^  T^~\r  P  TTT> 

e  e  V  e  neP  UB 

=-  - T - ^  (1-^  <14.49) 

Sx  ’e  1  +  B  2  H 

e 

The  sum  of  (14.46,  14.49)  then  yields  the  total  force  on  the  gas  heavy 
particles.  This  expression  is  identical  to  the  expression  (14.44),  neglect¬ 
ing  ion  slip. 

(ii)  Efficiency^  Hall  Generator 

The  local  efficiency  of  the  Hall  generator  is  the  ratio  of  the  power  (14.43) 
to  the  flow  work  U  when  is  given  by  (14. 44).  Thus, 

3’  ^  K  (1  -  K) 

^  (14.50) 

^  1  +  P' 

e  H 

2  2 
For  small  values  of  3  ’  ,  the  efficiency  is  proportional  to  3  >  and  is 

1  ®  2  ® 
maximum  when  ^  large  values  of  3^  ,  the  efficiency  becomes; 

g.  (14.51) 

e  i.  ^ 

that  is,  the  efficiency  is  a  maximum  when  K^is  small,  that  is,  when  close 
to  short-circuited.  For  arbitrary  values  of  3^  ,  as  K  is  increased,  the 
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efficiency  at  first  increases,  then  decreases,  as  shown  in  Figure  14.7. 

Thus,  for  any  given  value,  there  is  a  value  of  which  yields  the  maximum 
efficiency.  This  value  is  obtained  by  differentiation  of  (14.50)  with  respect 
to  K;  the  result  is; 


6' 

e 


2  K. 


m 


1  =  0 


which  yields 


K 

m 


Also, ^ from  (14.52); 


1-2  K 

m 


(14.52) 


(14.53) 


(14.54) 


Substitution  of  (14.54)  into  (14.50)  yields  V 

Lj 

m 


in  terms  of  K  ; 

m 


V  =  1  -  2  K  (14.55) 

mm 


Substitution  of  from  (14.55)  into  (14.  54)  yields  an  alternate  form 
for  the  maximum  local  efficiency; 


4  T) 


L 

m 


m 


(14.56) 


Thus,  3^  must  be  large  in  order  for  the  local  efficiency  to  be  close  to  unity. 
Figure  14. 8  shows  the  maximum  efficiency  as  a  function  of  Sj  given  by 
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LOCAL  STATIC  EFFICIENCY 


LOADING  FACTOR,K^ 

fl 


Figure  14.  7,  Local  Static  Efficiencies,  Hall  MHD  Generator, 

for  Various  Values  of  6  '  =  co  T  /(^  +05  T  05t  Tt  ) 

e  eel  eel  In 
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Figure  14.  8.  Maximum  Liocal  Static  Effic 


(14.56).  Note,  however,  that  as  the  magnetic  field  is  increased  from 

zero,  that  is,  as  3  is  increased,  3*  first  increases,  reaches  a  maximum , 
e  e 

and  then  decreases.  The  maximum  value  of  3  is  obtained  by  setting 

®  2 

equal  to  zero  the  derivative  with  respect  to3  of3/(l  +  3  M  M^).  The 

0  0  0  G  X 

result  is: 


b;  .  ,14.57) 

max  ’ 

Thus,  if  the  mobility  ratio  is  100,  the  maximum  value  of  3^  is  10  and  the 
maximum  local  efficiency  is  slightly  over  80%. 

(ili)  Comparison  With  Segmented  Electrodes 

It  was  shown  in  the  previous  section  that  the  segmented  electrode 

generator  has  a  higher  power  density  than  the  continuous  electrode  generator 

for  1  <  3  /M  ,  for  the  same  value  of  local  efficiency.  It  is  therefore 
0  0  1 

desirable  to  compare  the  power  densities  of  the  Hall  generator  to  that  of 
the  segmented  electrode  generator  for  the  same  values  of  V  From 
(14.38,  14.43), 


p 

Hall 

Seg 


8« 

e 


(14.58) 


where  K  =  K  ( ^7  ),  given  by  (14. 50).  This  ratio  is  shown  in  Figure  14.9 
as  a  fimction  of  the  parameter  3^  .  It  is  seen  that  for  the  same  efficiency, 
the  power  density  of  the  Hall  generator  approaches  that  of  the  segmented 
electrode  generator.  Of  course,  the  latter  generator  can  be  operated 
at  any  local  efficiency,  while  the  maximum  efficiency  of  the  Hall  generator 
depends  on  the  value  of  3g  ;  this  maximum  efficiency  is  shown  as  the  dashed 
line  in  Figure  14.9. 
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SEGMENTED 


Figure  14.  9. 


PoZ^r^nr"  Generator  to 

For  V.  ^  f  Segmented  Electrode  Generator. 
For  Various  Values  of  ^  -  //i  ,  . 


4Z 


Since  the  Hall  current  is  allowed  to  flow,  a  Lorentz  force  will  exist 
in  the  -y  direction,  as  in  the  case  of  continuous  electrodes.  The  ratio  of 
this  transverse  force  to  the  retarding  Lorentz  force  is 


F 

_JL  ^  _ 

F  1  +  0'  K 

X  e  H 


which  is  of  order  of  magnitude  unily,  for  3^  >  1;  this  must  be  balanced  by 
a  pressure  gradient  of  the  gas  normal  to  the  plane  of  the  electrodes  similar 
to  continuous  electrodes. 

(iv)  Voltage  Control,  Hall  Generator 

Since  the  maximum  power  density  is  developed  when  K  =  -^  ,  see 
(14.43),  and  the  efficiency  is  a  maximum  when  K  (see  Figure  14.7), 
the  Hall  generator  will  generally  be  operated  on  the  rising  part  of  the  power 
curve,  see  Figure  14.5.  This  means  that  if  the  generator  load  consists  of 
resistances  inparallel,  as  the  load  is  reduced  by  removing  resistances, 
the  total  load  resistance  will  increase,  increasing  the  load  voltage.  If 
the  initial  point  of  operation  was  K  then  the  new  value  of  the  loading 

1  b  2 

factor  may  still  be  less  than  —  ,  so  that  the  power  to  the  remaining  load  is 

now  greater  than  the  original  power  to  all  resistances.  Morever,  as  K  . 

H 

increases,  the  efficiency  decreases,  in  contrast  to  the  previous  two 
geometries.  This  may  be  seen  as  follows:  The  load  voltage  is  given  by 
K  3 '  U  B  L,  where  L  is  the  generator  length;  the  load  current  is  cr  h  b  P  UB* 

b  ®  r  2  2l  -]  ® 

(1-K)  (1  +  3  3  )  +  3  ;  thus,  the  load  resistance  R  =  V/I,  is  given 

^  L.  6  1  6  J  l_j 

by 


“  *^i  1-K 


(14.59) 
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where  the  internal  resistance  of  the  Hall  generator  is  defined  as: 


R. 

1 


abh  (1  +  Pg  Pj) 


(14.60) 


Solving  (14.59)  for  K,  one  obtains: 

H 


K  = 
H 


(14.61) 


For  K  ,  R.  >  ;  thus,  the  external  resistance  must  be  less  than 

H  2  1  L  - 

the  internal  resistance  to  obtain  good  local  efficiencies.  This  is  also  in 

contrast  to  flie  previous  cases.  Thus,  as  R^  increases,  K  will  also  in- 

L  H 

crease,  and  ^  will  decrease.  To  remedy  this  situation,  the  load  resis- 
tances  can  be  connected  in  series  and  resistances  can  be  removed  from  the 
circuit  by  shunting  them,  or  the  magnetic  field  can  be  reduced  to  keep  the 
load  voltage  constant,  if  the  load  resistances  are  in  parallel.  The  analysis 
of  the  voltage  control  is  simplified  if  the  ion  slip  is  neglected.  Then,  to 
keep  the  axial  electric  field  constant,  the  magnetic  field  must  be  varied 
so  that 


(14.62) 


where  K^,  are  the  initial  or  design  values  of  the  loading  factor  and 
magnetic  field,  respectively.  From  (14.59,  14.60),  the  load  resistance 
is  given  by: 


L  (1  +  e/)  K 

p  =  _ 6  " 

L  abh  (1  -  K) 

H 


(14.63) 
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and  the  design  load  is 


L 

abh 


(1  +  e  K 
o  o 

(1-Ko) 


(14.64) 


The  ratio  of  off-design  load  resistance  to  load  resistance  is  obtained  by 
dividing  (14. 63)  by  (14.  64).  Thus,  the  off-design  load  factor  K  is  given  by: 


(14.65) 


where  use  has  been  made  of 
given  by; 


e  =  8  ^K. 
e  o  o 


The  off-design  efficiency  is 
A 


Kh(1-K4 

1 


1-K. 


8 


+K 
2  H 


1 


8 

o  o 


+  1 


(14.66) 


For  example,  consider  a  Hall  generator  which  has  been  designed  for  8^  = 
5,  K  =  0.2,  V  .  Suppose  the  load  resistance  is  increased  by  38.5%, 
and  the  magnetic  field  is  reduced  to  keep  the  load  voltage  constant.  From 

(14.  65)  the  off-design  value  of  K  becomes  0.4;  from  ( 14.  62),  the  magnetic 

H 

field  and  8^  are  reduced  by  0.707,  and  the  efficiency  drops  to  0.5,  from 
( 14. 66).  If  the  magnetic  field  had  been  kept  constant  while  the  load  resis¬ 
tance  was  increased  by  the  same  amount,  the  load  factor  would  have  in¬ 
creased  to  0.53,  and  the  efficiency  would  have  been  reduced  to  0.436. 


(v)  Annular  Hall  Generator 

Since  there  is  no  electric  field  in  the  y-direction,  it  is  possible  to 
"wrap"  this  geometry  into  a  cylindrical  annulus;  the  axial  coordinate  is 
then  x;  the  y  coordinate  becomes  the  circumferential  coordinate  ©,  and 
the  z-coordinate  becomes  the  radial  coordinate,  see  Figure  14. 10.  In 
this  device,  the  flow  is  introduced  in  the  axial  direction  and  the  magnetic 
field  is  in  the  radial  direction.  In  practice,  it  is  not  possible  to  make 
radial  magnetic  fields  which  have  large  axial  lengths;  it  is  necessary  to 
periodically  reverse  the  direction  of  the  radial  magnetic  field,  see  Reference 
2  for  one  method  of  accomplishing  this.  However,  the  power  generation  or 
efficiency  does  not  depend  on  the  direction  of  the  radial  magnetic  field;  in 
fact,  it  may  be  shown  that  the  power  density  depends  only  on  the  average  of 
the  square  of  the  magnetic  field. 

On  the  other  hand,  the  circumferential  Lorentz  force  on  the  gas,  which 
is  caused  by  the  axial  current  j^,  does  depend  on  the  direction  of  the  radial 
magnetic  field.  If  the  radial  magnetic  field  were  to  be  in  the  same  direc¬ 
tion,  this  circumferential  force  would  cause  the  gas  to  rotate;  that  is,  the 
flow  would  become  helical.  But  if  the  direction  of  the  radial  magnetic 
field  were  alternated,  then  the  direction  of  the  circumferential  force  would 
also  alternate  in  the  x  direction.  This  would  tend  to  minimize  the  rotation 
of  the  gas. 

For  this  geometry,  the  local  performance  is  given  by  the  previous 
formula,  if  the  annulus  is  narrow. 

14. 3d  Series-Segmented  and  Helical  Flow  Generators 

In  the  preceding  sections  it  has  been  shown  that  the  segmented  electrode 
has  the  combination  of  the  hipest  power  density  and  highest  efficiency  of  the  three 
configurations  which  have  been  discussed;  however,  it  suffers  from  the  require¬ 
ment  of  multiple  loads.  This  problem  can  be  relieved  by  connecting  the  electrodes 
which  are  at  the  same  potential  in  series,  see  Figure  14. 11.  The  lines  of  constant 
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1— d 


SEGMENTED  ELECTRODES 


Figure  14.11.  Connection  Diagram  For  the  Series  -  Segmented 
Electrode  Geometry,  tan  0=  a‘‘  = 


potential  are  normal  to  the  electric  field  shown  In  Figure  14.  6.  The  require¬ 
ment  for  connection  of  ( 
from  (14.32)  becomes: 


ment  for  connection  of  electrodes  is  that  =^<P,orEh  =-E  d,  which 

y  X  y  X 


d  ^  K 

h  (1-K)  P' 
e 


(14.67) 


The  connections  shown  in  Figure  14. 11  are  for  a  single  load;  if  additional 
electrodes  are  inserted  between  those  shown  and  connected  in  series,  two  loads 
may  be  used,  and  so  forth. 

(i)  Maximum  Power  Density 

Now  it  is  natural  to  inquire  whether  there  is  some  other  linear  geometry 
which  can  have  a  higher  power  density  than  the  three  geometries  which  Have 
been  considered  for  a  given  efficiency.  To  determine  this,  let 


E  =  K,  UB 
X  1 

E*  =  K„  UB 

y  2 


(14.68) 


so  that  E  =  UB  +  E*  =  (K„  +  1)  UB.  The  current  can  be  obtained  in  terms 
y  y  2 

of  K  ,  K  from  (14. 18),  and  the  power  becomes: 

1  M 


=  T)  j  B  = 
L  •'y 


UB 


(1  +  P  Pj  (1  +  P" 

el  e 


(14.69) 


For  V  constant,  differentiation  of  (14. 69)  with  respect  to  K  shows  that  the 

1j  2 

power  density  is  a  maximum  when 


(14.70) 
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Now,  the  efficiency  is  also  given  by  {-E  j  -  E  j  )/j  UB,  or; 

XX  y  y  y 


(14.71) 


For  V  constant,  differentiation  of  (14.71)  yields; 

J_j 


d  K  d  K 

2K.  —  =  -  e;  (1-^3^)  -2K,-l+r?^ 


1  dK^ 


(14.72) 


Finally,  substitution  of  from  (14.71)  and(d  K^/d  K2)from  (14.70)  yields; 


Substitution  of  this  value  of  into  (14. 71)  yields  the  following  value  of  K^; 

K,  =  S'  (Tj  -  1)  =  S'  K  (14.74) 

1  e  '  L  e  2 

Use  of  (14. 68)  then  shows  that 


E  =  e'  E*  (14.75) 

X  e  y 

so  that  the  axial  current  is  zero,  which  is  precisely  the  condition  in  the 
segmented  electrode  geometry.  Thus,  for  a  linear  generator,  the  segmented 
electrode  geometry  yields  the  highest  power  density  for  any  given  efficiency. 
This  condition  may  also  be  interpreted  as  follows;  the  highest  power  density 
and  efficiency  is  obtained  when  the  Lorentz  force  exactly  opposes  the  gas 
velocity. 
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(ii)  Load  Variation  -  Series  Segmented 

The  series-segmented  geometry  is  designed  for  just  one  loading,  since 
the  electrode  pitch  d  is  selected  on  the  basis  of  the  loading  factor  K.  Thus, 
once  K  is  selected. 


(14.76) 


remains  constant  regardless  of  the  loading,  where 


(14.77) 


and  is  the  local  efficiency  for  the  design  load,  from  (14.  67).  For  an 
off-design  load,  E  =  K  UB,  so  that  E  =  -  a  K  UB,  and  the  load  voltage  is 

y  ^ 

V  =  K  UBh  (1  +  L/d)  ~  K  UBtt  L  (14.78) 


where  L  is  the  generator  length.  The  load  current  is  the  sum  of  the  current 
to  an  electrode  plus  the  axial  current; 

I  =(-iyd  (14.79) 

From  (14.18,  14.76), 


.  ^  _ g  UB 

y  (1  +  3  0)  (1  +  fi' 2) 

el  e 

This  may  also  be  shown  by  using  I  =  P/V,  where  V  is  given  by  (14.78); 
P  =  +  E yjy)h  Lb;  and  use  of  (14.76). 


K  -  1  -B 


'  2 


l-T) 
^ _ O 

T? 


K 


(14.80) 
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(14.81) 


a  UB  3' 

(i  +  3^ej)(i  +  3^  ) 

so  that  the  load  current  is: 

j.  • _ g  UBdb 

(l  +  3^3j)  (1  +  3^  ^ 

(14.82) 


(-K  .  7)^) 


The  value  of  K  which  corresponds  to  <^en  circuit  is  obtained  by  setting  the 
bracket  in  (14.82)  to  zero: 


K 

o.c. 


(14. 83) 


q’  2  ^ 

which  can  be  greater  than  unity.  In  fact,  for  3^  1,  the  open  circuit 

value  of  K  =  ^  /(l-'H  ).  Thus,  if  the  design  value  of  V  is  greater  than 
oc  o  o  o 

1/2,  K  is  greater  than  luiity.  The  power  generated  is  given  by  the  product 
oc 

of  (14. 78)  and  (14. 82);  the  ratio  of  the  off-design  power  to  the  design  power 
is  given  by: 
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For  3'  ^  »  1  this  reduces  to 
e 


This  has  a  maximum  when  K  is  equal  to  ^  /ll-V  ),  that  is,  half  the  open- 

2  o  o 

circuit  voltage;  this  maximum  is 


P 


4  r)  (i-T?  ) 
o  o 


(C) 


Since  ^  always  greater  than  4,  the  maximum  power  is 

always  greater  than  the  design  power.  Now,  if  t?  ^  ,  the  operating  point 

will  always  be  less  than  half  the  open-circuit  voltage.  Thus,  the  design 
point  will  be  on  the  rising  part  of  Ihe  power  curve,  see  Figure  14. 5.  In 
this  respect,  the  series-segmented  electrode  generator  resembles  Ihe  Hall 
generator.  If,  however,  the  design  value  of  V  is  less  than  ,  the  design 

O  it 

point  lies  on  the  falling  side  of  the  power  curve. 

(.lA.TS) 

The  load  resistance  is  given  by  the  ratio  of  the  load  voltage  to  Ihe  cur- 
rent;  the  design  value  corresponding  +© 


R, 


© 


abh 


3  (1  +  3  B  )  =■ 


e  r 


(Tbd 


^'0*Peftr) 


(d) 


while  the  off-design  load  is  related  to  the  load  factor  K  by: 


(e) 
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For  3’  ^ 
e 


»1,  the  reciprocal  of  (e)  becomes: 


+  1 


V 

o 


(f) 


which  shows  that  as  the  load  resistance  increases,  the  value  of  K  increases. 


The  efficiency  is  the  ratio  of  the  power  to  j  UB,  or: 


(14.84) 


For  large  values  of  3^  this  becomes: 


1  -  K 


(14.85) 


Thus,  as  the  load  resistance  increases,  the  load  voltage  will  increase,  and 
the  efficiency  will  decrease.  This  behavior  is  also  similar  to  that  of  the 
Hall  generator. 

(iii)  Helical  Flow  Generator 

It  was  shown  in  the  preceding  section  that  the  series-segmented  electrode 
geometry  has  the  advantage  of  the  high  efficiency  and  power  density  of  the 
segmented  electrode  geometry,  but  has  only  two  terminals  for  a  single  load. 
In  comparison  to  the  annular  Hall  geometry,  it  is  mechanically  complex 
because  of  the  necessity  for  having  multiple  electrodes  and  connections. 

Consider  the  electrode  arrangement  shown  in  Figure  14. 11,  where  the 
electrodes  are  along  constant  potential  lines.  If  the  electrodes  were  re¬ 
arranged,  as  shown  in  Figure  14. 12,  the  power  density  and  performance 
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would  be  the  same  as  for  Figure  14. 11;  that  is,  the  total  power  and  load 

voltage  are  identical.  Note  that  the  ionized  gas  must  pass  through  the  line 

of  electrodes  C-D.  Now,  this  geometry  is  symmetrical  about  the  line  A-A. 

Thus,  the  performance  is  not  changed  if  it  is  wrapped  into  an  annulus  with 

points  C  and  C  meeting,  and  points  D  and  D'  also  meeting.  Since  the 

equipotential  lines  now  meet,  the  electrodes  at  ^  ™3,y  be  dispensed 

with;  only  the  electrodes  at  and  <P  ^  are  now  needed  for  cormection  to  the 

*  o  4 

external  load.  This  new  geometry  is  shown  in  Figure  14. 13;  it  is  essen¬ 
tially  the  same  as  the  annular  Hall  generator,  except  that  now  the  gas  must 
be  admitted  with  a  circumferential  velocity  V0  ,  such  that  Vq/u^  =  dA,  or: 


u 

X 


tan  6 


n 

o 


(l-V  )  6* 
o  e 


(14.86) 


which  can  be  accomplished  by  inlet  turning  vanes.  The  axial  electric  field 
is: 


E 

X 


u  B  t]  (1  +  tan^6) 
X  o 

tan  6 


B' 

e 


(1  +  tan^B)  u  B 

j _  X 

1  +  b'  tan  0 
e 


(14.87) 


while  the  axial  electric  current  density  is  given  by: 

Cil-V  )  u  B  tan  6 
o  X _ 

X  1  +  60^ 

e  1 


(14.88) 


The  efficiency  has  been  taken  as  in  (14. 86  -  14.  88),  where  the  efficiency 
is  now  defined  as: 


V 


X  •’x 


^X  X 


u  •  (j  X  B) 


(u  j  -  V  j  )  B 
X  y  y  X 


(14.89) 
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Figure  14.  13.  Helical  Flow  Generator 
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The  power  density  is  obviously  3^,  which  is  obtained  from  (14. 87,  14. 88), 
while  (14. 89)  may  be  used  to  obtain  j^.  The  load  voltage  is  -  E^L  and  the 
load  current  is  the  product  of  j  and  the  annular  cross-sectional  area. 

X 

The  load  resistance  is  obviously  equal  to  the  ratio  of  load  voltage  to  load 
current. 

Althou^  the  flow  is  helical,  the  Lorentz  force  is  parallel  to  the  stream¬ 
lines,  so  that  the  Lorentz  force  does  not  change  the  ratio  of  V0/u^,  that  is, 
the  helix  angle  9  will  remain  constant.  The  helix  angle  is  determined  by 
the  selected  efficiency  which  in  turn  is  determined  by  the  load  resistance. 
If  the  load  resistance  is  changed,  the  current  and  voltage  will  change  as 
described  in  the  previous  section;  but  also,  the  direction  of  the  Lorentz 
force  will  be  different  than  the  flow  helix  angle  9.  This  can  be  remedied 
by  changing  the  turning  angle  of  the  inlet  stator  vanes  so  Uiat  tan  6  = 

Xj 

(1-^,),  where  9  is  the  new  helix  angle.  Also,  voltage  control 

J_(  6  Xj  Xj 

may  be  obtained  by  varying  the  magnetic  field  which  in  turn  changes 
and  6  again. 

The  helical  flow  geometry  has  two  possible  disadvantages:  first,  the 
helical  flow  will  cause  a  radial  pressure  gradient,  which  if  sufficiently 
large  will  alter  the  performance.  Also,  to  keep  the  helix  angle  small, 
large  values  of  are  necessary.  Second,  the  radial  magnetic  field  must 
be  in  the  same  direction;  switching  the  direction  of  the  magnetic  field  re¬ 
verses  the  required  helix  angle.  If  it  is  necessary  to  reverse  the  direction 
of  the  radial  magnetic  field  in  the  downstream  direction  as  in  the  annular 
Hall  geometry,  then  additional  turning  vanes  are  required  at  each  point 
where  the  magnetic  field  is  reversed.  It  will  not  be  necessary  to  reverse 
the  magnetic  field  if  the  generator  length  L  is  sufficiently  short,  which  re¬ 
quires  a  large  magnetic  interaction  parameter.  It  therefore  appears  that 
the  helical  flow  generator  is  most  applicable  to  large  values  of  both  the 
Hall  parameter  and  magnetic  interaction  parameter. 
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14. 3e  Magnetically  Induced  Ionization 

Whenever  an  electric  current  flows  in  a  gas,  the  "temperature"  of  the 
charge  carriers,  refered  to  the  mass-average  velocity,  will  be  higher  than  the 
uncharged  particles.  This  temperature  depends  on  the  loss  of  energy  of  the 
charge  carrier  when  it  suffers  a  collision  with  an  uncharged  particle;  in  Chapter 
V  it  was  shown  that  this  effect  is  most  pronounced  in  li^t  charge  carriers, 
that  is,  electrons.  As  the  electrons  enter  the  geijerator,  their  temperature  is 
changed  in  accordance  with  the  electron  energy  equation  as  follows; 
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where  is  the  gas  temperature,  n^  is  the  local  ionization  or  deionization  rate, 
and  Tj  is  the  ionization  temperature  equal  to  e  E^A,  where  is  the  ionization 
potential. 

From  (14.90),  the  relaxation  length  to  change  the  electron  temperature 


X 


> 


(14.91) 


where  is  the  electron  mean-free  path,  and,^  C 
of  the  heavy  particles.  Since  the  electron  velocity 
is  about  the  same  as  the  gas,  the  factor  at  c. 

the  gas  Mach  number;  thus,  (14.91)  becomes; 


\is  the  mean  thermal  speed 
'M 

in  the  streamwise  direction 
'sis  approximately  equal  to 
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(14.92) 
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If  TH-  6~  1,  then  (14.92)  shows  that  the  relaxation  length  is  of  the  order  of 
2  3 

10  -  10  mean-free  paths,  which  in  an  MHD  generator  is  much  shorter  than 
the  magnetic  interaction  length  given  by  (14. 26).  If  the  asymptotic  temperature 
is  considerably  higher  than  the  gas  temperature,  then  considerable  additional 
ionization  may  occur,  which  makes  n^  positive,  and  decreases  the  rate  of  in¬ 
crease  of  the  electron  temperature  somewhat. 

The  asymptotic  electron  temperature  is  obtained  by  setting  the  first  and 
second  term  on  the  right  side  of  (14.90)  equal  to  zero.  This  tacitly  assumes 

«  I 

that  the  relaxation  length  (14.92)  is  much  less  than  the  length  over  which  the  gas 

temperature  or  density  changes  appreciably,  that  is,  the  magnetic  interaction 

length.  Equation  (14. 26).  The  relaxation  length  for  ionization  is  of  the  order  of 

-4  -2 

n  /n  of  the  interaction  length;  since  n  /  is  of  the  order  of  10  to  10  ,  the 

e  e  n 

assumption  appears  to  be  valid.  Note  also  that  j  •  E*  is  always  positive,  so 
that  the  electron  temperature  is  always  greater  than  the  gas  temperature. 

To  obtain  the  asymptotic  electron  temperature,  the  relation  between 
and  E*  must  be  specified.  Since  this  relation  depends  on  the  generator  geometry, 
this  will  be  analyzed  for  the  three  basic  geometry  shapes.  Note  that  the  electron 
conduction  current  is  given  by; 


Thus,  the  electron  heating  depends  only  on  the  magnitude  of  the  electric  field, 
(i)  Continuous  Electrodes 

In  this  geometry,  the  only  component  of  the  electric  field  is  E*  = 
(K-1)  UB.  The  component  of  the  electron  current  in  this  direction  is 


bO 


jg^  =  a  E*/(l  +  Pg  )•  Equatiiig  the  sum  of  the  first  two  terms  on  the  ri^t 
side  of  (14.90)  to  zero,  substitution  of  E*  and  3  ,  and  solving  for  the 

electron  temperature,  one  obtains: 


T  y  (1-K)^  e 

~  =  1  +  - (14.95) 

g  3  6(1  +  3^2^ 

Where‘S  is  the  Mach  number  of  the  flow.  Now,  it  is  more  convenient  to 
express  the  electron  temperature  in  terms  of  the  stagnation  temperature 
of  the  flow,  given  by: 
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The  ratio  of  the  electron  temperature  to  the  local  stagnation  temperature 
is  therefore; 
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(14.97) 


In  order  for  the  electron  temperature  of  the  gas  to  exceed  the  local  stagna¬ 
tion  temperature,  it  is  necessary  for  the  coefficient  of^^  in  the  numerator 

2 

of  (14.97)  to  exceed  the  coefficient  of  9?^  in  the  denominator;  that  is  ^ 
de-f\ned  \>y  r 


(14.98) 


bl 


must  exceed  unity.  To  make  0  large,  it  is  necessary  for  0^  to  be  large 
and  0  small.  In  addition,  the  generator  should  be  close  to  short-circuited 
so  that  K  ~  0,  Under  these  conditions,  the  maximum  value  of  $  is  5/3. 

Thus,  if  the  Mach  of  the  flow  is  large,  T^/T^— >  5/3. 

Although  this  electron  temperature  may  cause  additional  ionization, 
it  is  not  practical  to  operate  an  MHD  generator  at  short  circuit.  Note 
also  that  6  must  be  small  for  appreciable  electron  heating. 

From  Chapter  V,  values  of  S  close  to  unity  may  be  possible  in  monatomic 

gases  at  hi^  electron  densities.  However,  in  combustion  gases  in  which 
2  3 

6  is  10  -  10  ,  no  appreciable  electron  heating  will  occur, 

(ii)  Segmented  Electrodes 

In  the  segmented  electrode  generator,  the  total  Hall  current  is  zero 
so  that  the  electric  field  from  (14. 32)  is  given  by: 

E*  =  -  (1-K)  UB 
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Substitution  of  (14.99)  into  (14,94)  yields: 
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(14. 100) 


In  the  numerator  of  the  bracket,  0^  may  be  neglected  in  comparison  to  0^; 
substitution  of  (14.  100)  into  (14.90)  yields: 


T  y  (i-K)^  0^7rL 

=  1  +  - ^ ^  (14.101) 
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for  the  asymptotic  electron  temperature;  the  ratio  of  the  electron  temper¬ 
ature  to  the  local  stagnation  temperature  is  then; 


y  (i-k)^b 

’^o  1+^  {7-1)711^ 


(14.102) 


Comparison  of  (14. 102)  with  the  corresponding  expression  for  continuous 

2  2 

electrodes,  (14.97),  shows  that  the  factor  )  which  is  always 

less  than  unity  is  replaced  by  3  ^/(l  +  3  3  ),  which  can  be  much  greater 

^  ®  ®  ^  2 
than  unity.  In  fact,  as  the  magnetic  field,  i.e.,  3  is  increased,  3  /(i  +  3  3  ) 

1/2  ^  eel 

has  a  maximum  when  3  =  (M  /M  )  ,  which  can  be  of  flie  order  of  10  -  20. 

e  el 

The  maximum  value  of  the  factor  is  therefore; 

2 


3 


B 


max 


1  _®. 
4  M, 


(14.103) 


max 


which  is  of  the  order  of  25  -  250.  Thus,  considerable  electron  heating  is 
possible  in  a  segmented-electrode  generator.  If  the  ionization  follows  the 
electron  temperature,  extremely  high  electrical  conductivities  may  be  ob¬ 
tained.  For  example,  calculations  have  been  made  for  a  mixture  of  argon 
and  1%  molar  potassium  for  a  stagnation  temperature  of  1500°K  as  a  func¬ 
tion  of  the  local  stagnation  pressure  and  Mach  number  of  the  resulting 
ciurrent  density,  see  Figure  14. 14.  The  corresponding  values  of  3^  = 
are  shown  in  Figure  14. 15.  The  parameter  is  (1-K)  UB,  where  U  is 
related  to  the  local  Mach  number  by; 
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^ykT/M  = 


(14. 104) 
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and  the  local  pressure  is  given  by 


(14. 105) 


when  the  temperature  ratio  is  given  by  (14.96).  In  these  calculations,  the 
effect  of  ion  slip  was  neglected;  the  perfect  gas  law  was  used,  and  the 
effect  of  the  ionization  energy  on  the  energy  equation  was  neglected.  The 
electrical  conductivity  is  obtained  by  dividing  the  current  density  by  the 
transverse  electric  field  (1-K)  UB.  The  resulting  electrical  conductivity 
is  several  orders  of  magnitude  greater  than  that  of  the  equilibrium  elec¬ 
trical  conductivity  at  T^.  Note  also  that  the  electrical  conductivity  increases 

with  increasing  Mach  number,  which  merely  indicates  that  the  coefficient 
2 

in  the  numerator  of  (14. 102)  is  greater  than  (y-l)/2.  Calculations 
at  other  stagnation  temperatures  indicate  that  the  electrical  conductivity 
is  approximately  proportional  to  the  stagnation  temperature.  Thus,  mag¬ 
netically  induced  ionization  may  be  of  interest  in  closed  cycle  applications 
of  MHD  generators.  As  of  the  present  time,  however,  magnetically  in¬ 
duced  ionization  has  not  been  demonstrated  experimentally.  Experiments 

8  9 

using  an  electric  field,  however,  have  verified  the  basic  theory.  ’ 

(iii)  Hall  Generator 

For  the  Hall  generator; 
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(14. 106) 


For  the  Hall  generator,  the  components  of  the  electric  field  are  given  by: 
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Use  of  (14. 106,  14,94,  14. 90)  yields  the  following  expression  for  the 
electron  temperature: 
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As  the  magnetic  field  is  increased,  the  first  factor  3^  /(I  +  8^  )  rapidly 

reaches  a  maximum  of  unity.  The  factor  in  the  bracket  reaches  a  maximum 
2 

when  6^  =  thus,  this  value  of  3^  yields  the  maximum  electron  temp¬ 

erature.  Note  that  this  maximum  is  the  same  as  that  for  the  segmented 

2 

electrode  geometry.  The  maximum  value  coefficient  of  >"%  /3  then 

12  2  12  2 
1  +  —  K  3^  ,  which  compares  to  —  ( 1-K)  3^  for  the  segmented  electrode 

generator.  Since  in  the  Hall  generator  K  is  comparable  to  (1-K)  in  the 


segmented  electrode  generator,  the  maximum  electron  temperature  is 


approximately  the  same  in  both  generators. 
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14.4  POLYTROPIC  EFFICIENCY  OF  MHD  GENERATORS 

The  local  efficiency  which  was  defined  in  the  previous  chapter  is  not  useful 
for  thermodynamic  cycles;  the  reason  being  that  the  previous  definition  was 
based  on  local  conditions  instead  of  stagnation  conditions.  The  local  efficiency 
was  defined  as 
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(14.108) 


Since  the  power  density  is  -E  •  j  and  the  ohmic  heating  per  unit  volume  is  E*  •  j, 

rs^  ^  rsf 

this  definition  is  the  same  as 


^  _  _ Power _ 

L  Power  +  Dissipation 

The  denominator  of  (14.  108)  is  also  equal  to  the  ideal  power  oulput  corresponding 
to  zero  dissipation  for  the  same  pressure  drop  due  to  Lorentz  forces,  that  is, 
for  the  same  current  density.  Thus,  the  ideal  power  is  given  by: 


P_  .  =  (-E  •  j)_ .  ,  =  ~  (E*  -  V  X  B)  .  j 

Ideal  ~  Ideal  ~ 


(14. 109) 


For  the  same  Lorentz  force  but  zero  ohmic  dissipation,  E*  must  be  zero,  which 

yields  the  denominator  of  (14. 108).  The  local  efficiency  can  also  be  represented 

on  a  Mollier  diagram  along  the  static  line  A-B,  see  Figure  14. 16.  Consider  two 

axial  stations  1  and  2  along  an  MHD  generator  which  are  a  distance  x  apart. 

As  the  gas  passes  between  these  two  stations,  the  total  enthalpy  per  unit  mass  is 

changed  by  ^H,  which  is  also  equal  to  the  change  in  static  enthalpy  if  the  change 

in  gas  velocity  is  negligible,  while  the  static  pressure  change  is  p  -  p  .  The 

1 

ideal  change  in  static  enthalpy  between  these  two  pressures  is  ^  h..  The  local 
efficiency  may  therefore  also  be  defined  as: 

(14.110) 
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Figure  14.  Ki.  Relation  between  Local  and  Polvtropic  Efficiencies 
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However,  the  thermodjmamicist  is  usually  interested  in  the  stagnation 


behavior  of  the  gas,  as  along  the  line  The  reason  for  this  is  that  the  gas 

usually  starts  at  a  small  Mach  number  in  a  combustor  or  heat  exchanger,  and 
after  being  accelerated  in  a  nozzle  and  passing  through  the  generator,  is  again 
brought  almost  to  rest  to  recover  the  pressure  or  to  pass  the  gas  through  a  heat 
exchanger  or  some  other  low  speed  device.  Thus,  one  is  interested  in  the  stag¬ 
nation  conditions  of  the  gas  in  the  generator.  The  local  efficiency  calculated  on 
the  basis  of  stagnation  conditions  is  called  the  poly  tropic  efficiency  some¬ 
times  called  the  "stage"  or  "local  adiabatic"  efficiency.  From  Figure  14.  16, 

V  =  Now,  because  of  the  spreading  of  the  constant  pressure  lines 

P  1 

with  increasing  enthalpy,  ^H.  ^  ^h.;  hence, 

A  slightly  different  definition  of  the  polytropic  efficiency  is; 

d  H 

1 

where  d  is  the  increment  in  total  enthalpy  due  to  electric  power  generation^ 

-  E  •  j,  and  d  H.  is  the  isentropic  or  ideal  change  in  total  enthalpy  between  the 
same  increment  in  pressure.  The  difference  in  d  and  the  actual  change  in 
total  enthalpy  d  H  is  due  to  heat  transfer  from  the  gas.  In  this  definition,  the 
heat  transfer  is  charged  as  a  loss  against  the  generator,  although  this  energy 
may  be  utilized  elsewhere  in  the  cycle.  This  definition,  of  course,  is  not  an 
adiabatic  efficiency.  For  the  ideal  or  isentropic  generator,  there  is  zero  heat 
transfer,  the  first  law  of  thermodynamics  yields; 

dQ.  =0  =  d  H.  -  — ^  (14.112) 

*  ‘  's 

where  p  and  are  the  stagnation  pressure  and  density,  respectively.  Sub- 
O  S 

stitution  of  (14.  112)  into  (14,  111)  and  use  of  the  ideal  gas  law  for  the  stagnation 
density  yields; 
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(14.113) 
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The  ratio  d  H  /d  H  is  the  fraction  of  the  total  enthalpy  which  is  converted  into 
e 

electrical  energy;  the  remainder  is  heat  transfer  from  the  generator;  let  this 

fraction  be  a  =  d  H  /d  H.  Also  d  H  =  C  T_,  and  C  /R  =  y/(  V-l)  where  7  is 
e  p  S  p 

the  ratio  of  specific  heats.  Then  (14. 113)  may  be  written  as; 
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Then  (14. 114) 


H 

H 


may  be  integrated  to  obtain; 


(14.115) 


where  the  exponent  of  the  stagnation  pressure  ratio  is  an  average  defined  by  the 
integration  of  (14. 114)  as  follows; 
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(14. 116) 


Now,  in  the  complete  MHD  generator,  the  stagnation  enthalpy  progresses  from 

I  to  3,  see  Figure  14. 16,  that  is,  from  H  to  H„.  The  ideal  change  in  enthalpy 

o  3 
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between  the  same  stagnationpressures  is  given  by  -  Hg,.  The  overall  generator 
efficiency  is  defined  as: 
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(14.117) 


since  only  the  fraction  a  of  (H  -  H  )  is  converted  into  electrical  energy.  The 

O  u 

relation  between  H  ,  and  H  is  obtained  from  (14. 116)  with  V  =  a  =  1: 

O  O  P 


(14. 118) 


while  H  is  obtained  from  (14. 11^  directly.  The  generator  efficiency  therefore 

<3 

becomes; 


(14. 119) 


For  small  pressure  ratios  across  the  generator,  that  is,  for  p  /p  ~  1,  the 
generator  efficiency  given  by  (14. 119)  becomes  just  But  for  large  pressure 
ratios,  the  pressure  ratio  terms  in  the  (14. 119)  become  small,  and  the  generator 
efficiency  approaches  a;  if  the  pressure  ratio  is  sufficiently  large,  the  generator 
efficiency  becomes  independent  of  the  polytropic  efficiency.  This  effect  is  well 
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known  in  the  design  of  steam  turbines  and  is  the  reason  for  designing  them  for 

very  large  pressure  ratios.  The  reason  for  the  high  efficiency  is  as  follows; 

the  dissipation  in  each  stage,  or  at  each  section  of  the  generator,  appears  as 

heat  energy  which  is  passed  onto  the  next  stage  or  section;  the  main  effect  of 

the  dissipation  being  a  loss  of  total  pressure.  But,  if  the  pressure  ratio  of  the 

machine  is  large  enough,  the  loss  of  pressure  is  small,  due  to  the  flattenii^  of 

the  constant  pressure  lines  at  low  pressure  on  the  Mollier  diagram,  as  shown 

in  Figure  14. 17.  For  moderate  pressure  ratios,  ^  is  equal  to  a  (H^  -  H^)/ 

(H  -  H„.).  For  large  pressure  ratios,  V  is  shown  as  a  (H  -  H„)/(H 
o  2*  or-  g  o  3  o  3' 

From  Figure  14. 17,  it  is  seen  that  the  generator  efficiency  is  larger  for  the 
larger  pressure  ratio.  Note  that  large  pressure  ratios  imply  large  temperature 
ratios,  since  (14. 119)  can  also  be  written  as; 


o 


W\'‘V>  -VVie  use  of 


(14. 120) 


However,  in  a  MHD  generator,  in  which  combustion  gases  are  used  at  pres¬ 
sures  greater  than  one  atmosphere,  the  lowest  total  temperature  at  which  the 
gas  still  has  sufficient  electrical  conductivity  is  about  4000°R;  the  highest  possible 
flame  temperature  is  about  SOOO'^R.  For  these  conditions,  it  is  important  that 
the  polytropic  efficiency  be  as  high  as  possible  because  the  generator  efficiency, 
as  calculated  from  (14. 119)  or  (14. 120),  will  be  only  slightly  larger  than 
The  generator  efficiency  will  be  improved  appreciably  over  the  polytropic  effi¬ 
ciency  if  the  pressure  ratio  is  very  large,  which  requires  some  method  for 
improving  the  electrical  conductivity  at  low  temperatures. 

14.4a  Polytropic  Efficiency  in  Linear  MHD  Generators 

In  this  section,  the  general  expression  is  derived  for  the  polytropic  efficiency 
in  linear  MHD  generators  under  steady  flow  conditions  at  a  particular  generator 
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Figure  14.17.  Effect  ( 
Genera 


station  x,  where  the  flow  velocity  in  the  x-direction  is  u  (x) ,  and  the  transverse 
velocity  components  are  negligible.  As  before,  the  local  magnetic  field  is 
transverse  to  the  generator  in  the  z-direction.  It  is  assumed  that  the  magnetic 
Reynolds  number  is  small  so  that  perturbations  of  the  imposed  magnetic  field 
may  be  neglected.  If  current  is  allowed  to  flow,  then  there  will  be  a  Lorentz 
force  in  the  direction  of  the  flow,  and  power  generation  per  unit  volume  equal 
to  -E  •  j.  In  addition,  there  is  a  friction  pressure  drop  (dp/dx)  which  is 
negative,  given  as  follows; 

/'dp\ 

W)f  ‘  -  “ 


where  T  is  the  average  shear  stress  at  the  wall;  C  is  the  perimeter  length  of 

w 

the  cross-section,  and  A  is  the  cross-sectional  area;  each  of  these  will  generally 
be  a  function  of  the  generator  station  x.  In  general,  a  friction  factor  c^  is  de¬ 
fined  as  follows; 
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In  addition  to  the  skin  friction,  there  may  also  be  heat  transfer  from  the 
gases  in  the  MHD  generator  to  the  walls;  this  is  given  by; 


C/A 


where  q  is  the  average  heat  transfer  rate  per  unit  area  over  the  perimeter  of 

w 

the  cross-section.  The  heat  transfer  rate  is  generally  related  to  the  film  coef¬ 
ficient  h  as  follows; 
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where  T  is  the  local  bulk  gas  stagnation  temperature  and  T  is  the  wall  temp- 
S  W 

erature;  the  film  coefficient  is  usually  made  non-dimensional  as  follows: 


®t  Pu  C 


where  is  the  local  Stanton  number.  Now,  if  the  Prandtl  number  of  the  gas, 

is  close  to  unity,  as  is  the  case  for  most  gases,  if  the  degree  of  ioniza¬ 
tion  is  less  than  a  few  per  cent,  then  as  a  general  rule  the  Stanton  number  is 
equal  to  half  the  friction  factor,  *  or: 


~  1 
2  ®f 


The  one-dimensional  momentum  and  energy  equations  therefore  become: 
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where  b  is  the  fraction  of  the  pressure  drop  due  to  the  Lorentz  force: 


(14.123) 


*The  presence  of  a  magnetic  field  may  alter  this  relation,  if  the  magnetic  inter¬ 
action  parameter  is  much  larger  than  unity. 
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and  a  is  die  fraction  of  the  enthalpy  drop  which  is  due  to  electrical  power  genera¬ 
tion  as  was  previously  defined.  From  (14. 122), 


dH  E  •  j 
_ e  ^  ^ 

dx  Pu 


(14. 124) 


and  from  (14. 113),  the  polytropic  efficiency  is; 


p  dH  /dx 

T)  =  — ^  - 

p  RT  dp  /dx 


E*  j 

4^ 


Pu  RTg  d  (-tn  Pg)/dx 


(14.125) 


To  evaluate  d  (-tn  p  )/dx,  use  is  made  of  the  following  isentropic  relation; 


y 

y-i 


or 


•tnpg  =  ■^np+:^  ^ 


(a) 


Differentiation  with  respect  to  x  yields; 


dx  p 


y 

y-i  T 


+ 


y 

y-i 


(b) 


where  (  )'=  (d/Ax).  Next,  substitute 

H'  -  C  T'  +  uu'  (c) 

p 


7b 


in  the  energy  equation,  and  eliminate  c  uu'  between  the  energy  equation  and  the 
momentum  equation.  The  result  is: 


El  _2L  II  ^  1  PiLSi 

p  ”  y-1  T  p  L  b 


(d) 


also,  from  the  energy  equation, 


'T't  _ 

S  apu  C 

P 


Substitution  of  (e,d)  into  (b),  and  then  into  (14. 125)  yields: 


(e) 


V 

P 


u 


a 


•  (3  X  B) 

/w  ^ 

E  •  j 


a 

b 


(f) 


since  lt(3  X  B)=  u  •  (j  X  B);  hence,  the  first  factor  in  the  denominator 
is  just  the  reciprocal  of  the  previously  defined  "local"  electrical  efficiency 
Thus,  the  polytropic  efficiency  becomes; 


n 

P 


a  T, 


V  bT 


(14.  126) 


where  T  /T  is  the  ratio  of  the  local  stagnation  temperature  to  static  temperature, 
which  for  a  perfect  gas  is  given  by: 


(14. 127) 
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where  the  local  Mach  nvimber.  To  proceed,  expressions  are  required  for 
a,  b.  Several  possible  cases  will  be  examined  in  the  next  three  sections: 
a  =  b  =  1,  corresponding  to  zero  friction  and  zero  heat  transfer;  a  =  1,  b  ^  1, 
corresponding  to  zero  heat  transfer  but  finite  friction,  and  a  1,  b  1,  corres¬ 
ponding  to  both  heat  transfer  and  friction. 

14.4b  Polytropic  Efficiency  for  Zero  Friction  and  Zero  Heat  Transfer 
For  this  case  a  =  b  =  1,  and  use  of  (14. 127)  in  (14.  126)  yields; 


r) 

P 


(14,  128) 


Thus,  for  the  polytropic  efficiency  to  be  close  to  the  local  efficiency  either 
V  must  be  close  to  unity  or  else  (7-  must  be  small.  Thus,  the  local 

J_ 

efficiency  is  identically  equal  to  the  polytropic  efficiency  of  a  subsonic,  zero 
heat  transfer,  zero  friction,  flow.  For  combustion  gases  where  7  1.2,'>^~  1, 

and  7)^  =  0, 8,  the  polytropic  efficiency  is  only  2%  less  than  But  for 
gases  in  which  7  =  5/3,^^)^=  2,  and  =  0.8,  the  polytropic  efficiency  is  almost 
25*^  less  than  the  local  efficiency  V 

14.4c  Polytropic  Efficiency  With  Friction  But  Zero  Heat  Transfer 

With  friction  but  zero  pressure  drop,  a  =  1,  and  from  the  definition  of  b. 


1  ^  ‘Tc/A 

b  “  (jxB) 

rs/  ^ 


(14.  129) 


To  proceed  further,  an  expression  is  needed  for  (j  x  B)  ,  which  requires 

/W  /W  X 

that  the  geometry  be  specified.  We  will  consider  all  three  basic  geometries. 

Actually,  since  the  local  efficiency  is  equal  to  K  for  both  the  continuous  and  segmented 
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electrode  geometries,  they  may  be  considered  as  one  case  with  an  effective 
transverse  electrical  conductivity  which  is  evaluated  as  follows; 


Continuous  Electrodes:  o  „ 

eff 


_ 0 _ 

(1  +  p  P  +  P  ^ 

el  e 


Segmented  Electrodes: 


_ a 

^eff  1  +  P  P 

e  I 


Since  the  local  efficiency  is  different  for  the  Hall  generator,  this  will  be  con¬ 
sidered  separately. 

(1)  Continuous  and  Segmented  Electrode  Generators 

2 

For  these  two  geometries,  (j  x  B)  =  (K-1)  a  uB  ,  so  that  the  ex- 

X  cli 

pression  for  b  becomes; 


1 

b 


1  + 


'f 


2  I  (1-K) 


(14.130) 


^farawe'Ver 

where  I  is  the  magnetic  interaction  based  on  a  length  equal  to  A/C  as 


follows; 


1  = 


DuC 


(14.131) 


Finally,  a  parameter  Z  may  be  defined  as: 


^  =¥1 


(14.132) 
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so  that  b  is  given  by; 


1 

b 


=  1 


+ 


Z 

1-K 


(14. 133) 


Substitution  of  (14. 133)  into  the  expression  for  the  polytropic  efficiency 
yields: 


T? 

P 


1 

1-K  +  Z  1  ^ 
K  (1-K)  “  T 


(14, 134) 


Now,  for  any  given  value  of  Z  and  T  /T,  there  is  a  value  of  K  which 
minimizes  the  bracket  in  (14. 134)  and  which  therefore  yields  the  highest 
polytropic  efficiency.  This  value  of  K  =  K^  is  given  by: 


K  =  1  +  Z  - 
m 


Z  +  Z 


2 


(14. 135) 


or  alternatively, 


Z 


2  K  -  1 
m 


(14. 136) 


The  maximum  polytropic  efficiency  is  therefore: 


2  K 

m 


pm 


1-(1-K  )(7-l)'^^ 
m 


(14.  137) 
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For  subsonic  flows  where  the  Mach  number  is  negligible,  ~  2 

substitution  of  K  in  terms  of  V  into  (14.  136)  then  yields; 
m  pm 


{i-n  y 

pm 

1  ^  ^pm 


(14. 138) 


which  is  identical  to  the  expression  for  the^ipaximum  local  efficiency  of  a 

‘j'vtrt  by  ^  2 

Hall  generator,  with  Z  replacing  I/S'  .  This  indicates  that  small  values 
A  ® 

of  Z  are  required  to  achieve  high  values  of  the  polytropic  efficiency.  Now, 

I  based  on  A/C  will  seldom  be  more  than  0.1,  while  the  friction  factor  f 
is  generally  inversely  proportional  to  the  Reynolds  number  raised  to  the 
power  of  1/2  for  laminar  flows  and  0.  2  for  turbulent  flows.  *  The  result 
is  that  very  large  flow  Reynolds  numbers  are  required  to  obtain  high  values 
of  the  polytropic  efficiency;  that  is,  Z  must  be  small. 

For  Z  much  less  than  unity,  and  finite  Mach  number,  the  above  ex¬ 
pressions  may  be  simplified  as  follows: 


Z  1; 


(14. 139a) 


ri 

pm 


(14.  139b) 


which  clearly  shows  that  the  polytropic  efficiency  decreases  with  both 
increasing  Z  and  Mach  number. 

In  the  jireceding  discussion,  it  was  assumed  that  the  Mach  number 

was  held  constant.  However,  in  an  MHD  generator,  the  mass  flow  rate 

is  usually  siieeifiecl.  It  is  therefore  of  interest  to  determine  the  manner 

in  which  V  changes  with  Mach  number  for  a  constant  mass  flow  rate, 
pm 

Mt  has  been  assumed  here  that  the  magnetic  field  does  not  affect  the  friction  factor. 
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For  simplicity,  consider  a  square  cross-section  which  has  a  characteristic 
dimension  h,  so  that 


A/C  =  7  h,  and  h  =  V m/ P  u 
4 

where  m  is  the  total  flow  rate.  Also,  the  skin  friction  factor  c^  is  inversely 
proportional  to  the  Reynolds  number  to  a  power  n; 


c 


f 


GM 


n 


.  n/2 
m 


(Pu) 


-n/2 


where  n  =  0. 5  for  laminar  flow  and  0. 2  for  turbulent  flow,  and  G  is  a 
constant.  Also, 


Du 


s 


Thus, 


(rp,/a,)' 

eff 


n 

4 


(14. 140) 


The  terms  which  depend  explicitly  on  Mach  number  increase  with 
increasing  Mach  number,  as  is  shown  in  Figure  14. 18. 
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Figure  14.  18 


If  the  effective  electrical  conductivity  also  decreases,  then  the  maxinaum 
polytropic  efficiency  given  by  (14. 140)  will  decrease  with  increasing  Mach 
number.  On  the  other  hand,  if  die  electrical  conductivity  increases  with 
increasing  Mach  number,  as  for  example  is  the  case  with  magnetically 
induced  ionization,  then  the  effect  of  increasii^  Mach  number  on  the  poly¬ 
tropic  efficiency  is  mitigated. 

(li)  Hall  Generator 

M 

From  (14.44)  and  (14. 123)  the  egression  for  b  becomes: 

1 

Z  -  I  *  I  - 1 - 1  z  (14.141) 

SO  that  the  expression  for  the  polsrtropic  efficiency,  (14. 126)  becomes: 


(1+Z)  +  K  +  Z 

3  ^  J 

®  (14.142) 

This  expression  is  essentially  the  same  as  that  given  by  (14. 128)  in  the 

absence  of  friction,  except  that  in  computing  the  local  efficiency  of  the 

Hall  generator,  3’  should  be  divided  by  the  factor  1  +  Z  +  3  Z 
©  © 

14. 4d  Polytropic  Efficiency  With  Friction  and  Heat  Transfer 
For  this  case,  the  factor  a  defined  as 

-1  ,  V 


=  1  1 
n  ^  ■  T  K(l-K) 
p 
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is  now  larger  than  unity.  Use  of  for  a  crossed 

field  qert^rAf-ot'  Q^=  "'*‘w  ),  ^t=l^/puc  C  - ZS^  ^ 

ana?  T  ^  a- 

_ 

'^S  /  (14. 143) 

where  Z  is  defined  by  (14. 132).  Then,  use  of  (14. 133,  14. 143)  in  (14. 126)  yields; 


n 

p 


+  1 


z 

^  K(l-K) 


(14. 144) 


The  last  term  in  the  above  expression  represents  the  effect  of  skin  friction  and 

is  identical  to  14. 134,  except  that  T  has  replaced  T  in  the  last  term.  Thus, 

W  D 

wall  cooling,  that  is,  T  T  ,  leads  to  a  reduction  in  the  pressure  drop  due  to 
friction.  This,  of  course,  is  the  well  known  effect  in  one-dimensional  compres¬ 
sible  flow,  namely,  that  cooling  increases  the  stagnation  pressure. 
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14.  5  END  LOSSES  IN  LINEAR  MHD  GENERATORS 


In  addition  to  the  friction  losses  described  in  Section  14.4,  there  may  be 
electrical  losses  at  the  inlet  and  exit,  caused  by  eddy  currents  which  flow  from 
the  positive  electrode  to  the  negative  electrode  through  the  conducting  gases  in 
the  region  of  the  channel  which  is  .external  to  the  main  magnetic  field,  that  is, 
upstream  and  downstream  of  the  electrode  region,  see  Figure  14. 19.  At  present, 
only  the  simplest  analyses  of  this  problem  have  been  made  because  of  its  com¬ 
plexity.  The  assumptions  generally  made  are: 

1.  The  average  electron  collision  time  is  constant. 

2.  The  presence  of  viscous  boundary  layers  is  neglected.  The  basis  for 
this  assumption  is  that  the  viscous  Reynolds  number  will  be  large  in 
MHD  generators,  so  that  the  viscous  boundary  layer  will  be  quite  thin. 

3.  The  gas  velocity  is  a  constant,  taken  as  U  in  the  downstream  (x)  direc¬ 
tion.  This  assiunption  really  contains  two  assumptions: 

a)  The  flow  is  essentially  incompressible,  that  is,  the  Mach  number 
is  small. 

b)  The  induced  Lorentz  forces  do  not  appreciably  alter  the  velocity 

profile.  This  assumption  is  valid  if  the  magnetic  interaction 

2 

parameter  based  on  the  channel  width,  aB  h/PU,  is  small  com- 
10 

pared  to  unity. 

4.  The  magnetic  field  in  the  generator  is  equal  to  the  imposed  magnetic 
field.  This  assumption  is  valid  if  the  magnetic  Reynolds  number 

a  Uh  is  small,  which  is  generally  true  in  MHD  generators. 

5.  The  imposed  magnetic  field  is  in  the  transverse  direction  only  of 
the  magnitude  is  a  function  only  +h«t  lorttjiVw.dirv:.\  ccofUimiiin  the 
region  near  the  ends  of  the  magnetic  field,  that  is,  in  the  "shaded" 
region.  This  assumption  violates  Maxwells'  equations,  but  the  error 
is  not  serious  if  the  pole  face  separation  distance  is  sufficiently  small; 
the  results  are  qualitatively  correct. 
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Negative 

electrode 
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Electric  Current  Lines  Due  to  Electric  Fields  in  the  Entrance  Region  of  the  Channel 


6.  The  insulating  sides  are  perfect  electrical  Insulators  and  the  electrodes 
and  leads  have  zero  electrical  resistance. 

7.  With  the  above  assumptions,  the  problem  is  two-dimensional  in  the 
plane  normal  to  the  magnetic  field. 

8.  The  flow,  magnetic  field  and  electric  field  are  steady.  The  problem  is 
then  to  calculate  the  electrical  current  distribution  in  the  generator, 
the  power  generaced,  the  pressure  drop,  and  the  overall  efficiency 

of  the  generator. 

14. 5a  Basic  Equations 

The  basic  equations  are. 

k 

V  xE  =  0  (14.145) 

V  •  j  =  0  (14.146) 

Since  the  flow  is  steady,  E  =  -  V  ^ .  Thus,  from  (14. 145),  the  electric  potential 
is  continuous,  while  from  (14. 146),  the  electric  current  is  continuous. 

The  boundary  condition  at  the  insulator  is  that  the  normal  component  of 
the  electric  current  must  be  zero.  In  the  geometry  of  Figure  14. 19,  this  means 
that 


j^fX,  h)  =  i^(X,  0)  =  0;  X<  0 


(14. 147) 


At  the  electrodes,  the  boxmdary  conditions  depend  upon  the  type  of  electrodes 
which  are  used.  If  the  electrodes  are  continuous,  then  they  are  equipotential. 

On  the  other  hand,  if  the  electrodes  are  segmented,  then  one  may  specify  either 
(a)  the  current  distribution  along  the  electrodes,  (b)  the  transverse  voltage  dif¬ 
ference  (X,  0)  -  (X,  h),  or  (c)  the  ratio  of  the  voltage  to  the  current,  which 

is  the  same  as  specifying  the  distribution  of  external  load  resistance.  In  the 
ensuing  sections,  we  will  consider  both  continuous  electrodes  and  segmented 
electrodes  for  generators  which  use  the  Faraday  current. 
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In  the  generator,  there  usually  exists  an  insulated  inlet  section,  the  elec¬ 
trode  region,  and  an  insulated  outlet  section;  in  principle  all  three  regions  should 
be  analyzed  simultaneously.  However,  it  has  been  found  that  if  the  generator 
aspect  ratio  a,  defined  as  the  ratio  of  the  electrode  region  to  the  generator 
width  h,  is  greater  than  about  0.7,  then  the  inlet  and  exit  regions  can  be  analyzed 
separately,  and  the  results  combined  to  yield  the  generator  efficiency;  that  is, 
the  exit  region  has  no  effect  over  the  inlet  region  and  vice  versa.  In  addition, 
if  the  velocity  profile  is  constant,  then  the  current  distribution  at  the  inlet  and 
exit  are  similar;  hence,  only  the  inlet  region  is  analyzed  in  this  section.  How¬ 
ever,  in  order  to  calculate  the  generator  efficiency,  both  the  inlet  and  exit 
losses  must  be  considered.  Although  in  general,  the  MHD  generator  will  have  a 
diverging  cross-section  so  that  the  inlet  and  exit  losses  will  be  different,  for 
simplicity  we  will  assume  that  the  inlet  and  exit  losses  are  identical;  that  is, 
we  will  consider  the  generator  to  have  a  constant  cross-sectional  area.  Thus, 
it  is  assumed  that  the  insulating  region  extends  an  infinite  distance  upstream  of 
the  electrode  region  (X  0)  while  the  electrodes  extend  infinitely  downstream. 

Finally,  the  physical  dimensions  are  made  non-dimensional  as  follows: 

X  =  XTT/h 

y  =  Y^T/h 

so  that  the  non-dimensional  channel  width  is  ^ . 

14.5b  Continuous  Electrodes,  Scalar  Electrical  Conductivity 

With  the  electrical  conductivity  scalar,  the  current  equation  becomes: 

j  =  a(E  +  vxB)  (14.148) 

rsj  fs* 

so  that  (14. 146)  becomes: 

y  2 

=  0  (14.149) 
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since  v  =  i  U  and  B„  =  B  (X).  At  the  electrode  at  Y  =  0,  the  potential  may  be 

rsj  Cj  La 

taken  as  +0,,,,  while  at  the  electrode  at  Y  =  h,  the  potential  is  taken  as  -  <P 

w  w 

so  that  the  potential  difference  between  the  two  electrodes  is  2  0  .  Along  the 

w 

insulator  at  X  <  0,  the  normal  component  of  tiie  current  is  zero  so  that  (14. 148) 
becomes 


d0 
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(X,  0)  =  -5^  (X,  h)  =  -UB^  (X) 


(X  <  0) 


(14. 150) 


We  will  consider  three  different  variations  of  the  magnetic  field  in  the 
region  upstream  and  downstream  of  the  electrode  region; 

(i)  B^  (X)  =  0;  X  <  0 


(ii)  B^  (X)  =  B^  exp  (  X/X);  X  <  0 


(iii)  B  (X)  =  B^;  -X^  <  X  <  0 
L  00 


B  (X)  =0  X<  -X^ 

L  O 


Case  (i)  corresponds  to  the  magnetic  field  ending  at  the  end  of  the  electrode 
region.  Case  (ii)  corresponds  to  an  exponential  decrease  in  the  magnetic  field 
beyond  the  electrode  region  with  an  e-folding  length  X.  Case  (iii)  corresponds 
to  an  extension  of  the  magnetic  field  at  constant  magnitude  B^  for  a  distance  X^ 
upstream  (or  downstream)  of  the  electrode  region. 

This  is  then  a  problem  of  solving  Laplace's  equation  in  an  infinite  strip. 
However,  it  is  complicated  by  the  boundary  conditions:  for  x  ^  0,  the  value  of 
the  potential  is  specified  at  boundaries,  while  for  x  0,  the  normal  gradient 
is  specified.  For  this  reason,  solutions  cannot  be  obtained  analytically  by  means 
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of  separation  of  variables.  However,  this  problem  may  be  solved  by  means  of 
various  techniques  which  utilize  complex  variables;  the  simplest  method  is  to 
conformal  map  the  strip  so  that  the  boundary  on  which  the  normal  gradient  is 
specified  is  perpendicular  to  the  boundary  on  which  the  potential  is  specified. 

The  simplest  mapping  which  accomplishes  this  transformation  is  shown  in  Figure 
14.20,  and  is  given  by; 

Z  =  'C-nsinZ'  (14.151) 

where  Z  =  X  +  i  y  is  the  original  complex  plane  and  Z '  =  X'+  iy*  is  the  trans¬ 
formed  plane.  In  the  Z'  plane,  Laplace's  equation  must  be  satisfied  with  respect 
to  x',  y';  the  potential  along  the  electrodes  remains  the  same,  but  the  normal 
gradient  must  be  appropr  i  ately  transformed.  The  three  cases  are  considered 
next; 

(i)  Zero  Extension  to  Magnetic  Field 

For  zero  magnetic  field  upstream  of  the  electrodes,  from  (14. 150) 

(  3  <p/  3y)  =0  along  the  boundaries;  hence,  in  the  transformed  plane 

(  3  (p/  cyi)  =  0.  The  solution  to  Laplace's  equation  in  the  Z'  plane  is  then 

obviously; 


W 


(14. 152) 


Next,  the  total  electrode  current  to  some  point  X^  can  be  calculated  as 
follows; 


where  (14. 148)  has  been  used.  The  first  term  in  the  bracket  can  be  made 
dimensionless,  and  the  second  term  integrated  as  follows; 
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0 

a.  PHYSICAL  PLANE 

I 

y 

~w 

11^^77777777////  ^ 

“  '2  0 

b.  TRANSFORMED  PLANE  Z  *  Ln  sin  z' 

Figure  14.20.  Conformal  Mapping  of  End  Losse 
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J  (XJ  = 
y  1 


-1 

-V  v'“- 


x)  dx  -  aUBX, 


Now  (  9<P/a  y)  d  X  =  -  (  a<p/  ^x')  dy'  along  the  tflectrode,  and  from  (14. 152), 
S<p/dx' =2  0^/n.  Thus, 


Jy(X^)  .  2(7(»^yyn  -oOBXj 


(14. 153) 


Along  th^  electrode  Z  =  x  and  Z*  =  ^/2  +  i  y',  so  that  (14. 151)  becomes 


X  =  -f-n  cosh  y* 


For  large  values  of  the  argument,  cosh  y*  ~  ^  exp  y'  so  that 


4  ~  (2  ~  yj 


Substitution  of  y^  in  terms  of  X^  into  (14. 153)  then  yields; 


2  ri" 

JIXj)  =  -B  |-jr~  *  ■^”2'  -OUBX 


1 


(14. 154) 


Now,  the  definition  of  the  load  factor  may  be  taken  as; 


K  =  2<P„/UBh 

w 

so  that  (14. 154)  becomes 

Jy(Xi)  =  -  o  UB  (1-K)  +  K  o  UBh  -tn  2  (14.155) 
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The  first  term  in  (14. 155)  is  the  "ideal"  current,  which  would  exist  in 
the  absence  of  end  losses;  the  second  term  represents  the  current  loss, 
due  to  the  shunting  effect  of  the  conducting  fluid  upstream  of  the  magnetic 
field.  Note  that  this  term  is  proportional  to  the  voltage  difference  between 
the  electrodes,  that  is,  the  load  voltage.  The  generated  power  up  to  point 
is  given  by  the  product  of  the  current  J  (X^)  and  the  load  voltage, 
KUBh.  Thus,  the  second  term  in  (14. 155)  also  represents  the  power  loss. 

If  the  exit  of  the  generator  has  the  same  dimensions  as  the  inlet,  a 
similar  loss  exists  at  the  exit,  thus,  the  total  lost  power  P  per  unit 
height  is; 

P  =  2  •  KUBh  •  K  UB  h-tn  2  (14.156) 

Xj 


while  the  ideal  power  P.  for  a  generator  of  length  L,  width  h,  per  imit 
height  is; 


Pj  =  OK  (1-K)  U^B^h  L  (14.157) 

The  actual  power  P  is  the  difference  of  (14. 157)  and  (14.  156); 

P^  =  oU^B^Kh  [(1-K)  L  -  2TT"^  h  ^n  2]  (14.158) 


To  calculate  the  efficiency,  the  flow  work  W  through  the  generator 
must  be  calculated; 


W 


V  j  X  B  d  Y  d  X 


0  0 


(14. 159) 


The  integration  is  taken  over  X  to  obtain  the  average  flow  work  across  the 
generator,  while  the  integration  over  Y  is  taken  only  where  the  magnetic 
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field  Is  non-zero.  With  v  =  i  U,  and  B  in  the  Z  direction,  the  only  compo- 
nent  of  current  in  (14. 159)  is  thus. 


W  = 


4/  +  UB^  BdYdX 

0  0 


(14. 160) 


Integration  with  respect  to  Y  and  then  X  and  use  of  the  definition  of  K 
yields; 


W 


-  a  (l-K)  U^B^  Lh 


(14. 161) 


which  is  identical  to  the  e^ipression  which  one  obtains  if  the  end  regions 
are  ignored.  The  efficiency  is  then  given  by 


V 


P 

a 

W 


(14.162a) 


from  (14. 158,  14. 161),  and  where  a*  is  a  modification  of  the  aspect  ratio. 


a*  =  ^  t  "^"2  (14.162b) 

Note  that  as  the  generator  becomes  very  long,  that  is,  as  a* — >  0,  then 
the  ejqaression  for  the  efficiency  becomes  identical  to  the  "local"  efficiency 
K.  Also,  K  =  0  still  corresponds  to  short  circuit,  but  open-circuit  condi¬ 
tions  correspond  to  the  bracket  in  (14. 16/2a.' equal  to  zero,  e.g. ,  when 


K 

o.c. 


1 

1  +  a* 


(14. 163) 


Thus,  the  larger  the  value  of  a*,  the  smaller  the  open-circuit  voltage. 
However,  when  the  voltage  is  equal  to  the  open-circuit  voltage,  the  effi¬ 
ciency  is  equal  to  zero,  rather  than  equal  to  K.  The  efficiency  is  therefore 
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greater  than  zero  only  for  0^K'^(l+a*)  A  typical  efficiency  cuirve 
is  shown  in  Figure  14. 21  for  an  aspect  ratio  LA  =  1.  Note  that  the  maximum 
efficiency  is  only  0.29,  corresponding  to  a  loading  factor  of  0,45.  As  the 
aspect  ratio  LA  increases,  a*  decreases,  and  both  the  maximum  efficiency 
and  loading  factor  at  maximum  efficiency  increase  .  For  any  given  aspect 
ratio,  the  maximum  efficiency  is  obtained  by  equating  to  zero  the  derivative 
of  (14. 162a).  The  maximum  efficiency  is  then  given  by; 


This  maximum  efficiency  is  shown  in  Figure  14. 22.  It  is  seen  that  large 
aspect  ratios  are  required  to  obtain  efficiencies  greater  than  0.6. 

For  a  given  aspect  ratio,  the  efficiency  can  be  improved  by  the  insertion 
of  insulating  vanes  at  the  inlet  and  exit  parallel  to  the  flow  but  external  to 
the  magnetic  field  region.  For  example,  a  single  inlet  vane  and  exit  vane 
placed  in  midstream  will  double  the  apparent  aspect  ratio;  two  equally 
spaced  vanes  will  triple  the  apparent  aspect  ratio,  etc. 

Another  meAod  of  improving  the  efficiency  is  to  extend  the  magnetic 
field.  Actually,  most  magnetic  fields  decrease  exponentially  away  from  the 
pole  faces.  Hence,  this  type  of  "shading"  is  considered  next. 

(ii)  Exponential  "Shading"  of  Magnetic  Field 

To  solve  this  case,  it  is  convenient  to  introduce  the  complex  potential 


f  (z)  =  <P  +  i  fp  (14.  165) 

such  that  ‘P,  obey  the  Cauchy- Riemann  conditions;  furthermore,  the 
equation  for  0  is  also 

V  2  , 

^  =  0  (14.166) 


9b 


GENERATOR  EFFICIENCY 


.8 


LOAD  FACTOR,  K 


Figure  14.  Zl,  Generator  Efficiency,  For  Channel  Length  Equal 
to  Channel  Width 
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max 


Figure  14.  2Z.  Maximum  Efficiency  Considering  End  Losses  - 
Zero  Extension  to  the  Magnetic  Field 
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Thus,  J/)  is  a  "stream  function"  for  the  potential.  Also,  the  contribution 
to  the  current  caused  by  the  potential  gradient  is  given  by  -V  qc: 


l(p 


X 


gTT 

h 


Re 


d  ^ 
dz 


j 


<P 


y 


gTT 

h 


Im 


d$ 

dz 


or; 


X  y 


gTT  d  $ 
h  dz 


(14. 167) 


Thus,  between  any  two  points  z  and  z  ,  the  total  current  between  these 

1  M 

two  points  is; 


=  -  g 


(14. 168) 


In  particular,  if  z  and  z  are  points  on  the  same  equipotential ,  then; 
1  ^ 


=  -ig  [V^(z^)-i^(z2)] 


or 


J 


(P 


X 


y 


-Im  (z^)  -  ^(Zg)] 
a  Re  [iP(z^)  -  (zj 


(14.  169) 
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Now,  for  this  case, 


B  =  B  e"  ;  x  <  0 
o 


(14, 170) 


so  that  along  the  insulator,  for  to  be  zero: 


h  yx 

—  =  -  -  UB  e 

d  y  n  Q 


(14. 171) 


where 


y  =  h/TT  X 


In  the  z'  plane,  along  y'  =  0, 


5  <P  _  ^  Q  Ip 

S  y'  3  y  (J2'  9x' 

Use  of  (14. 151)  and  (14. 171)  in  (14. 172)  and 
x'  yields: 


(14.  172) 


integration  with  respect  to 


ip 


.  >  , 

~yjr~  sin  x' 


(14. 173) 


Now,  Ip  must  obey  Laplace's  equation  in  the  transformed  plane;  thus,  the 
solution  may  be  written  as: 


Ip 


2<P  y' 
W 


T 

n=0 


b  cos  ;inx'  e 
n 


-2ny' 


(14. 174) 
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The  first  term  above  satisfies  the  boundary  condition  at  the  electrodes;  it 
is  in  fact  the  solution  obtained  in  case  (i).  The  second  term  represents 
the  Fourier  series  solution  to  Laplaces  equation.  Only  the  cosine  series 
is  used  because  the  sine  series  does  not  satisfy  the  boundary  conditions 
along  the  electrodes.  The  cosine  series  thus  represents  the  additional 
current  due  to  the  shading  of  the  magnetic  field;  this  will  be  designated  as 
Now,  from  (14.169),  the  additional  current  to  the  electrode  between 
(  ^/2,  0)  and  some  point  (^/2,  3;^)  where  y^  1  is; 


J  (y^)  = 


^  *  Yl)  -  ^  ( 2  ’ 


(14. 175) 


From  (14. 174),  the  first  term  in  the  above  bracket  is  just  the  first  term 
in  the  Fourier  series,  b^.  This  is  obtained  by  solving  for  the  Fourier 
coefficients  along  y'  =  0  from  (14. 173,  14. 174); 


h  UB 


,/,+  o  .  y  .  T-  .  „  .  -2nv' 

>P  =  ^  y —  sm  x'  =  2,  b^  cos  2  n  x'  e 


from  which 


2  h  UB 


lim  >1^  (-  .  y^)=  b^  = 

y'— 


^/2 


sin  x' d  x' 


\.*/  e  f  e  T*  *  S 
From  (14.  173), 


h  UB 
_ o 

try 


(14.  176) 


(14. 177) 
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Substitution  of  (14.  176,  14.  177)  into  (14.  175)  yields  the  additional  electrode 
current; 


■)  -  - 


o  UB  h 
_ o_ 

n  y 


(14.  178) 


The  additional  power  at  the  inlet  is  given  by  2  ‘P  J^.  If  the  configuration 

4.  t  W  ^ 

at  the  exit  is  identical  to  that  at  the  inlet,  an  equal  amount  of  additional 

power  will  be  generated  at  the  exit  due  to  the  shading  of  the  magnetic  field. 

a 

To  calculate  the  efficiency,  the  flow  work  is  calculated  in  manner 
analogous  to  (14. 160);^^ 


W 


2  2  2 
o  U  B  h 
o 


^  2 


-  ( 1-K)  L/h 


^  n=2,4,6 


A 

n 


1 

y 


t 

\ 


(14. 179) 


where: 


A 

n 


_ r-(y+i) _ 


The  efficiency  is  obtained  by  dividing  the  total  power  by  the  flow  work  given 
by  (14.  179);  this  is  shown  in  Figure  14.21  for  two  different  e-folding  lengths 
of  the  magnetic  field,  h  and  h.  It  is  seen  that  considerable  improvement 
in  the  efficiency  is  obtained  by  slight  shading  of  the  magnetic  field. 
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(iii)  Constant  Magnitude  Extension  of  the  Magnetic  Field 

12 

The  magnetic  field  may  also  be  extended  at  constant  magnitude  ; 
the  efficiency  for  an  extension  equal  to  ten  channel  widths  is  shown  as  the 
dashed  line  in  Figure  14.21.  It  is  apparent  that  the  effectiveness  of  this 
type  extension  is  poor  in  comparison  to  the  exponentially  shaded  field.  This 
is  due  to  the  sudden  termination  of  the  magnetic  field  at  a  point  where  the 
channel  is  an  insulator.  Thus,  almost  the  full  open-circuit  voltage  is 
developed  there.  This  leads  to  a  maximum  eddy  current  loss  of  power  at 
this  point,  since  this  loss  is  essentially  proportional  to  the  voltage  difference 
across  the  channel. 

14.  5c  Effect  of  Tensor  Conductivity 

In  general,  the  electric  current  vector  can  be  separated  out  into  two  parts: 
one  due  to  the  stationary  potential  gradient  V  and  the  other  due  to  the  induced 
field  V  X  B  as  follows: 

}  =  +  ir,  (14.180) 

This  procedure  has  already  been  utilized  in  the  preceding  section,  see  (14. 167). 
Since  the  values  of  v  and  B  are  assumed,  the  induced  current  is  obtained  easily, 
but  to  obtain  the  current  due  to  the  potential  gradient  requires  the  solution  to  the 
potential  problem.  The  boundary  conditions  for  the  potential  problem  are  usually 
either  a  specified  potential  or  normal  current.  When  the  electrical  conductivity 
is  scalar,  the  normal  current  boundary  condition  becomes  a  normal  gradient 
of  the  potential.  When  the  potential  or  the  normal  potential  gradients  are  specified 
on  the  boundaries,  the  problem  is  well  posed. 

However,  when  the  electrical  conductivity  is  tensor,  the  specification  of  the 
normal  current  does  not  lead  to  a  normal  gradient  of  the  potential.  For  example, 
consider  the  velocity  in  x  direction  only,  and  the  magnetic  field  in  the  z  direction 
only.  Then  the  current  components  become: 
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1  =  e  .  .  _  o UB 

^B  1  R  2  ’  ^B  "  “  ,  f,  2 

X  1+P  y  1+P 

e  e 


a<p  b(p 

■Sx  **6^? 


(14. 181) 


^<P 

“ay  "  e  Sx 


Thus,  if  j„  is  specified  as  zero  on  a  boundary,  then  (  ^  <P/  ^  Y)  =  -  3  (  ^  (p/  9x). 

_i  ® 

Thus,  the  equipotential  lines  have  a  slope  of  with  respect  to  the  x-axis. 

This  type  boundary  condition  is  neither  Dirlchelet  nor  Neuman,  and  therefore 

requires  special  treatment.  To  solve  this  type  of  problem,  one  may  conformal 

map  the  original  region  such  that  the  equipotential  lines  are  straight  and  equally 

spaced  in  the  transformed  plane.  The  boundaries  along  which  the  potential  lines 

are  sloped  are  then  mapped  at  the  appropriate  angle  to  the  equipotential  lines. 

Using  this  technique,  the  potential  map  for  the  continuous  electrodes  have  been 
0 

calculated  when  the  electrical  conductivity  upstream  of  the  electrode  region 
is  zero.  The  resulting  equipotential  and  current  flow  lines  are  shown  in  Figure 
14. 23,  for  =  1*  The  potential  and  current  due  to  the  induced  electric 

field  UB  have  been  omitted.  It  is  seen  that  there  is  a  current  concentration  at 
the  upper  left-hand  comer  of  the  electrode.  There  is  no  effect  of  the  distortion 
at  the  inlet  on  the  efficiency;  however,  the  total  power  is  reduced  somewhat  by 
the  end  effect. 


Using  this  same  technique,  the  potential  field  and  current  has  been  obtained 

when  the  magnetic  field  and  electrical  conductivity  are  extended  infinitely  up- 

13 

stream  of  the  electrodes  and  when  the  electrodes  are  skewed  with  respect  to 
one  another. 
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Figure  14.  23.  Potential  Mesh  (Solid  Lines)  and  Current 
Flow  (Dashed  Lines)  For  to  T  -  And 
0^  =  0  For  X  <  0.  The  Curr^nf  and  Potential 
Lines  Due  to  UB  Have  Been  Omitted. 


For  values  of  'T  in  excess  of  unity,  the  mapping  technique  has  also 
been  used  to  determine  the  current  flow  in  the  vicinity  of  segmented  electrodes, 
when  the  electrode  length  in  the  direction  of  flow  is  equal  to  the  length  of  insula¬ 
tion  between  electrodes.  The  calculated  current  distributions  are  shown  in 
Figure  14. 24.  It  is  seen  that  there  is  a  concentration  of  current  in  the  down¬ 
stream  portion  of  the  electrode;  and  also,  the  distortion  of  the  current  lines 
extends  into  the  stream  only  a  distance  equal  to  the  electrode  length. 

If  the  electrode  pitch  distance  were  infinitely  small,  the  effect  transvei^se 

conductivity  would  be  equal  to  the  scalar  electron  conductivity,  in  the  absence 

of  ion  slip;  the  effect  of  finite  size  electrodes  and  insulators  decreases  the 

15 

effective  transverse  conductivity.  The  theoretical  reduction  is  shown  in 

Figure  14. 25  for  ^  =  1»  where  it  is  seen  that  increasing  pitch  distances 

decrease  the  effective  transverse  conductivity.  In  addition,  increasing  values 

also  decreases  the  effective  transverse  conductivity  as  shown  in  Figure 

14. 26.  Also,  increasing  the  electrode  width  as  compared  to  the  electrode 
16 

pitch  is  deleterious.  Finally,  it  has  been  shown  that  the  viscous  boundary 

layer  along  the  electrode  surface  has  no  effect  on  these  results,  and  also, 

15 

finite  size  electrodes  do  not  effect  the  efficiency,  neglecting  end  losses. 

The  effect  of  end  losses  on  efficiency  with  segmented  electrodes  has 

17 

also  been  investigated  theoretically.  The  maximum  efficiency  is  shown  in 

Figure  14. 27  for  various  aspect  ratios,  where  it  has  been  assumed  that  the 

electrode  pitch  distance  is  infinitely  small.  In  the  vicinity  of  the  inlet  or  exit, 

one  may  specify  either  (i)  the  current  to  each  electrode,  (ii)  the  transverse 

potential  difference,  or  (iii)  a  relation  between  the  electrode  current  and 

potential  difference.  To  investigate  the  effect  of  these  different  conditions, 

calculations  were  made  for  w  T  =  o  for  (i)  constant  electrode  current  and  (ii) 

6  0 

constant  interelectrode  potential.  It  can  be  seen  that  the  efficiency  for  case  (i) 

is  greater  than  (ii).  Also,  it  can  be  shown  that  the  results  for  case  (iii)  will 

be  between  (i)  and  (ii);  hence,  only  the  "constant  current"  case  was  calculated 

for  finite  w  t  ,  it  is  seen  that  increasii^  w  t  increases  the  efficiency 
©  ©  ©  © 

slightly  as  does  extending  the  magnetic  field.  It  therefore  appears  that  the 
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CURRENT  DISTRIBUTION 


ure  14.24.  Current  Flow  Lines,  (a)  Segmented  Electrodes  with  ojT  = 

(b)  Segmented  Electrodes  With  cor  =  3;  Note  That  Most  of  the 
Current  Flows  Into  the  Right  Hand  Corner  of  the  Electrode. 

(c)  Infinitely  Thin  Electrode;  This  Case  is  Independent  of  («jr. 
The  Dashed  Lines  in  the  Figure  Represent  Current  Dividing 
Lines  Between  Adjacent  Electrodes. 
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OJT  =  I 


Figure  14.25.  Effect  of  Electrode  Pitch  on  Effective  Transverse 
Electrical  Conductivity  for  u:  e  "^e  '  ^ 
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Figure  14.Zfa.  Effect  of  On  Effective  Transverse  Electrical 

Conductivity  For  Segmented  Electrodes 


Fitjure  14.  Z7.  Maximum  Generator  Efficiency 


generator  efficiency  for  segmented  electrodes  and  finite  values  of  is 

greater  than  for  continuous  electrodes. 

The  influence  of  tensor  conductivity  on  end  losses  has  also  been  calculated 
for  the  boundary  condition  at  the  electrode  that  =  0-  Although  this  is  not  as 
realistic  as  those  used  above,  the  resulting  field  distribution  and  conclusions 
are  essentially  the  same. 
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14.6  COMPRESSIBLE  FLOW  IN  FARADAY  CURRENT  MHD  GENERATORS  ’ 

In  any  magnetohydrodynamic  generator  which  uses  an  ionized  gas,  the 

effect  of  compressibility  must  be  included  if  there  is  a  large  pressure  drop 

through  the  generator.  Thus,  the  results  of  the  previous  sections  may  be 
ou 

regarded  only  as  "local"  analysis  over  a  short  section  of  the  generator,  in 
A 

which  the  change  in  properties  is  small.  As  the  temperature  and  pressure 

change  in  the  generator,  so  will  the  degree  of  ionization  or  composition  and 

hence  the  gas  is  not  perfect,  and  the  analysis  of  the  flow  may  require  numerical 
21  22 

integration.  ’  However,  if  the  degree  of  ionization  is  small  (<  0. 1%),  then 
the  total  amount  of  ionization  energy  is  small  compared  to  the  gas  energy,  and 
the  compressibility  factor  is  also  close  to  unity.  Under  these  conditions,  one 
may  assume  that  the  gas  is  perfect. 

The  effects  of  friction  and  heat  transfer  may  also  be  neglected  if  the  Z 
factor  (14.132)  is  small;  which  implies  "large"  generator  dimensions.  Finally, 
if  the  interaction  parameter  (14.131)  is  small  compared  to  unity,  which  is 
generally  correct,  then  the  properties  of  the  ionized  gas  change  "slowly"  in 
flow  direction.  With  this  assumption,  one  may  use  the  assumption  of  quasi-one 
dimensional  flow;  that  is,  the  transverse  flow  velocities  in  the  channel  are 
small;  and  the  axial  flow  velocity,  magnetic  field,  electric  field,  and  electrical 
conductivity  do  not  vary  across  the  cross  section  of  the  generator.  The  assump¬ 
tion  concerning  the  magnetic  field  is  valid  when  the  magnetic  Reynolds  number 
is  small.  If  x  is  the  direction  of  flow  and  u  is  the  mass  averaged  gas  velocity 
in  that  direction,  then  the  appropriate  equations  are  as  follows: 

Continuity 

puA  =  constant  =•  m  (14.182) 


Momentum 


pu 


(jxB) 


(14.183) 
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Energy 


pu 


=  E 


•  j 


(14. 184) 


State 

p  =  pRT  =  nkT  (14.185) 

Caloric 

dh  =  C  dT  (14.186) 

P 

where  A(x)  is  the  cross-sectional  area,  is  the  specific  heat  at  constant 
volume  assumed  constant,  and  R  is  the  gas  constant, 

R  =  k/<m>  (14.187) 

when  <  m>  is  the  average  mass  of  the  particles: 

n 

<m>  =  E  — ^  m  (14.188) 

g  n  s 

In  addition  to  these  equations,  it  is  necessary  to  specify  Ohm's  law,  and  the 
manner  in  which  the  conductivity,  or  the  factors  which  enter  the  expression  for 
the  conductivity,  vary  with  temperature,  pressure,  etc.  These  equations 
(14. 182)  -  (14.186)  have  been  used  to  analyze  the  compressible  flow  in  Faraday 
type  generators;  that  is,  for  continuous  or  segmented  electrodes.  At  present, 
compressible  analyses  have  been  not  performed  for  the  Hall  type  geometry. 

The  Hall  geometry  generally  requires  an  additional  equation  for  conservation  of 
current  in  the  downstream  direction, 

j^A  =  constant  (14.189) 
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In  general,  it  is  also  necessary  to  specify  the  manner  in  which  the  transverse 
magnetic  field  varies  in  the  downstream  direction.  For  convenience  it  will  be 
assumed  that  the  magnetic  field  is  constant:  methods  for  including  the  variation 
of  the  magnetic  field  will  be  indicated.  Also,  the  electric  field  must  be  speci¬ 
fied;  in  general  the  electric  field  will  be  taken  as  a  fraction  of  the  transverse 
induced  electric  field  uB.  Even  with  the  magnetic  field  and  electric  field 
specified,  equations  (14.182)  -  (14.186)  are  not  determinate:  there  is  one  more 
unknown  than  there  are  equations.  It  is  therefore  necessary  to  specify  one  of 
the  variables  to  solve  these  equations.*  In  the  following  sections,  solutions 
will  be  presented  for  various  variables  such  as  velocily,  Mach  number,  tem¬ 
perature,  pressure,  and  area  specified  constant.  The  paths  which  are  taken 
by  keeping  these  variables  constant  are  shown  in  Figure  14.28. 


14.6a  Constant  Velocity  Generator 

With  the  axial  velocity  held  constant,  the  electrical  energy  is  extracted 
from  only  the  thermal  energy  of  the  gas.  Actually,  this  may  be  regarded  as  a 
two-stage  process,  since  the  Faraday  induction  removes  electrical  energy  only 
from  the  kinetic  energy  of  the  flow.  However,  the  gas  is  then  accelerated  to 
restore  the  velocity,  but  in  the  process  of  acceleration  the  static  temperature 
decreases.  The  momentum  and  energy  equations  become,  respectively; 


dx 


(j  X  B)  .  i 


(a) 


p  uC 

P 


dT 

dx 


E  •  j 


(b) 


♦There  is  a  school  of  thought  which  asserts  that  even  when  the  magnetic 
Reynolds  number  is  small,  one  must  include  Maxwell's  equation  jjj.  =  u 
However,  it  is  physically  more  plausible  that  the  current  in  the  gas 
causes  an  x-component  of  the  magnetic  field  and  that  (SB  /ax)  may  be 
neglected.  ^ 


ax 
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Figure  14.28.  Paths  on  a  Mollier  Diagram  For  Various  Linear 

MHD  Generators.  Point  0  =  Stagnation  Conditions 
Upstream  of  Generator;  1  =  Generator  Inlet;  2  = 
Generator  Exhaust;  3  =  Stagnation  Conditions 
Corresponding  to  Generator  Exhaust. 
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The  first  integral  of  equations  (a),  (b)  is  obtained  by  division  of  (a)  by  (b)  and 
integration  as  follows: 


V  -  1  df,n  p  _  1  _  Jj^ 

y  d-tnT  K 


(c) 


where  use  has  been  made  (14. 108),  and  the  fact  that  in  a  Faraday  type  gener¬ 
ator,  the  local  efficiency  is  equal  to  the  loading  factor,  whether  the  electrodes 
are  continuous  or  segmented.  Since  the  generator  is  expected  to  have  a  certain 
overall  efficiency,  it  is  desirable  to  have  the  jMjlytropic  efficiency  approximately 
constant;  hence,  the  loading  factor  K  will  be  chosen  constant.  This  also  insures 
that  every  section  of  the  generator  will  be  generating  electrical  power.  With 
K  and  y  constant,  equation  (c)  may  be  integrated  immediately. 


T. 


y 


(14.190) 


Equation  (14. 190)  may  be  compared  to  the  usual  isentropic  expression: 

xzA. 


Since  r).  is  always  less  than  unity,  (14.190)  shows  that  for  the  same  pressure 

J.J 

ratio,  the  temperature  decrease  is  always  less  than  isentropic.  This  is  a 
direct  result  of  the  ohmic  heating  of  the  gas. 

The  entropy  change  may  be  calculated  from: 

A  S  =  C  -In  ^  -  R  ^n  (e) 

P  Tj  p^ 
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Substitution  of  (14. 190)  in  (e)  yields: 


AS  _  1  - K 

C  K  T 

P 


(f) 


which  shows  the  rapid  increase  in  entropy  of  the  gas  with  decreasing  K.  This 
is  also  due  to  the  increased  ohmic  heating  associated  with  the  increased  elec¬ 
trical  current  density.  Equation  (f)  may  also  be  rearranged  as  follows: 


T 

=  ®^P 


KAS 
(1  -  K)  C 


(g) 


To  obtain  the  variation  of  the  gas  properties  in  the  downstream  direction, 
it  is  necessary  to  integrate  either  equation  (a)  or  (b),  making  use  of  the  gener¬ 
alized  Ohm's  law  for  a  Faraday  type  generator: 


E  =  KuB 

j  =  a^jjCE  -uB)=  (1  -  K) 


(14.191) 


where  the  magnetic  field  B  is  taken  in  the  Z -direction,  the  electric  field  E  is  in 
the  y-dlrection  and  is  given  in  Section  14.4b.  For  example,  substitution  of 
(14. 191)  into  (a)  yields: 


(14.192) 


Now  is  generally  not  a  constant,  but  depends  upon  temperature,  pressure, 
and  possibly  the  recombination  rate.  If  is  purely  a  "state"  variable,  then 
it  may  be  expressed  in  terms  of  pressure  alone  by  the  use  of  (14. 185,  14. 190); 
and  (14. 192)  may  be  integrated  directly.  Also,  if  B  varies  in  a  prescribed 
manner  with  x,  this  variation  may  also  be  included  in  (14. 192),  and  the  integration 
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performed.  A  different  method  of  handling  the  variation  of  conductivity  and 
magnetic  field  is  to  use  a  transformed  variable  defined  as  follows: 


'Vf 


dx 


(14. 193) 


where  the  subscript  1  refers  to  conditions  at  the  generator  inlet.  Then  (14. 192) 
becomes: 


which  may  be  integrated  immediately  to: 

Pj  -  P  (? )  =  (1  -  K)  (a^j^  B^)^  u^  (14. 194) 


Now  with  the  use  of  (14. 185,  14. 190),  p,  T,  and  consequently  may  be 
obtained  as  a  function  of  The  variation  of  ^  with  x  is  then  obtained  by  inte¬ 
gration  of  (14. 194)  as  follows: 


X 

(B/Bj)^  dx  (14. 195) 

0 


The  variations  of  B  and  may  therefore  be  regarded  as  a  "scaling"  of  the 
length  parameter  x;  that  is,  the  value  of  ^  which  is  required  to  obtain  a  given 
pressure  drop  is  equal  to  the  length  x  required  to  obtain  the  same  pressure  drop 
in  a  generator  in  which  the  magnetic  field  and  effective  conductivity  remain 
constant.  Decreases  in  either  the  magnetic  field  or  conductivity  increase  the 
length  X  required  to  obtain  the  same  pressure  drop. 
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In  spite  of  the  variable  conductivity  or  magnetic  fields,  (14. 19  )  may  be 
used  to  define  an  "interaction  length"  L,  which  is  the  length  in  units  of  ^  for 
which  the  exhaust  pressure  is  much  smaller  than  the  inlet  pressure  p^,  and 
hence  may  be  neglected.  Thus, 

Pi 

L  =  - - -  (14.196) 

(1  -K)  )jU 

«.  }  A 

The  interaction  length  is  therefore  a  measure  of  the  length  of  the  MHD  generator 
which  is  required  to  obtain  a  large  pressure  drop;  that  is,  an  appreciable  tem¬ 
perature  drop  of  the  gas .  Note  that  L  is  different  than  the  magnetic  interaction 
parameter  I  which  was  defined  in  order  to  compare  the  magnitude  of  the  Lorentz 
force  to  the  friction  pressure  drop. 

(i)  Optimum  Inlet  Mach  Number,  Constant  Velocity  Generator 

Now  if  L  as  calculated  from  (14.196)  is  large  compared  to  the  trans¬ 
verse  dimension  of  the  generator,  then  it  may  be  desirable  to  change  the 
inlet  conditions  to  decrease  the  generator  length:  the  most  convenient 
method  for  accomplishing  this  being  to  change  the  inlet  Mach  number, 
e.g. ,  the  inlet  static  temperature.  The  interaction  length  L  is  a  minimum 
when  the  derivative  of  u)j  is  set  to  zero.  For  example,  we  may 

consider  a  segmented  electrode  generator  with  thermal  ionization. 
Neglecting  electron-ion  collisions,  and  letting  the  collision  cross  sections 
be  constant,  the  scalar  electrical  conductivity  is: 


.,^^2  ^3/4  exp(-E./2kT) 

V 8m  <Q>p^^^  (n  /n)^^^ 

©  s 


(14. 197) 


where  it  is  assumed  that  the  degree  of  ionization  of  the  seed  material  S 
is  small.  The  gas  velocity  u  is  given  by: 


u  =  V2C  (T  -  TJ 
p  o  1 


(14. 198) 
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The  inlet  static  pressure  is  given  by; 


(14.199) 


Substitution  of  (14.197-199)  into  (14.196),  differentiation  with  respect  to 
T^,  and  equating  the  resulting  expression  to  zero  yields  the  following 
expression  for  for  which  L  is  a  minimum: 


o 


The  results  of  applying  (14.200)  is  shown  in  table  14.3,  for  two  typical 
gases:  combustion  gases,  for  which  y  ~  1. 20,  and  inert  gases,  for  which 
y  =  5/3.  These  results  indicate  that  the  minimum  length  is  obtained  for 
combustion  gases  when  the  flow  is  supersonic,  while  for  inert  gases,  sub¬ 
sonic.  It  is  to  be  expected  that  the  optimum  Mach  number  is  less  for  inert 
gases  than  for  combustion  gases  because  the  temperature  of  the  former 
decreases  more  rapidly  than  the  latter  as  the  gas  velocity  increases. 

For  magnetically  induced  ionization,  as  the  Mach  number  is  increased, 
the  conductivity  increases  up  until  ion  slip  begins  to  decrease  the  electron 
temperature,  hence  (14.200)  may  not  be  used  in  this  case.  Also,  (14.200) 
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is  not  a  useful  criteria  if  the  pressure  drop  due  to  Lorentz  forces  is 

small  compared  to  p^^;  or  if  the  interaction  length  is  already  sufficiently 

small.  In  this  case,  one  may  wish  to  optimize  the  power  density 
2  2 

K(1  -  K)  au  B  .  Then,  for  thermal  ionization,  constant  collision  cross 
section,  and  small  ionization,  the  optimum  inlet  static  temperature  is 
given  by: 


T 


o 


+ 

2  - 


2(62-1) 


(14.201) 


where 


e 


2 


_ 1— 

2(y  -  1) 


2kT 

o 


The  results  of  (14.201)  are  also  shown  in  Table  14.3.  Note  that  the 
optimum  inlet  Mach  number  for  combustion  gases  is  now  subsonic,  while 
the  optimum  Mach  number  for  the  inert  gas  is  increased  only  slightly. 


TABLE  14.3.  COMPARISON  OF  OPTIMUM  INLET 
MACH  NUMBERS 


E.  =  3. 87  ev 
1 


o 

Scsii 

Length 

Power  Density 

Gas: 

y 

'^1 

Combustion 

TTi 

1.36 

0.81 

hiert 

1.67 

0.55 

0.60 
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(ii)  Effect  of  Variable  Conductivity 

To  illustrate  the  direct  integration  technique  with  variable  conductivity, 
consider  a  thermally  ionized  gas  whose  electrical  conductivity  is  given  by 
(14. 197),  which  may  be  represented  as  follows: 


(14.202) 


nct 

It  is  convenient  to  use  an  expression  as  complicated  as  that  given  by 

y\ 

(14.202);  to  simplify  the  analysis,  let: 


o 


(14.203) 


where  ui,  v  are  exponents  which  are  obtained  by  matching  with  (14.202), 
that  is: 


V  - 


1 

2 


LC 


3 

4 


+ 


^n  (T/T^) 


(14.204) 


Although  cc  varies  with  temperature,  the  variation  is  small,  and  the  tem¬ 
perature  T  in  (14.204)  may  be  chosen  at  some  average  value.  With  the 
variation  of  conductivity  given  by  (14.203),  either  the  momentum  equation 

(14. 192)  or  the  energy  equation  may  be  integrated;  to  illustrate  this, 

ZA- 

eonsider  the  energy  equation  (for  constant  velocity): 


P  u  C 

P 


dx 


-  K  (1  -  K)  a  uB 


2 


(14.205) 
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law,  (l'+>‘85)^ 

Use  of  the  perfect  gas,  (14.190),  and  (14.203),  (14.205)  ifeductis +0; 
A 


-  1 


(14.206) 


where 


yp. 


= 


(y  -  1)  K(1  -  K)  B 


Equation  (14.206)  then  integrates  immediately  to: 


ax 

X* 


-1/a 


(14.207) 


and  with  the  use  of  (14.190)  the  corresponding  pressure  ratio  is  obtained. 

To  illustrate  the  use  of  the  scaling  law,  we  may  integrate  (14. 195) 
with  constant  magnetic  field.  From  (14. 194),  (14.190),  and  (14.203),  the 
conductivity  as  a  function  of  ^  is  given  by: 


-2. 

a. 


1  - 


^  (1  -  K)  uB  e 


w  (y-  1) 
y 


K  - 


(14.208) 
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integration  of  (14.195)  then  yields: 


(1  -  K)  uB  i 


+  1  -ooK 


1  - 


Oj  (1 


K)  uB 


(14.209) 


which  is  the  required  scaling  relationship.  Successive  substitution  of 

(14.209)  into  (14. 194)  and  (14. 190)  then  yields  the  same  relation  (14. 207) 
as  obtained  by  direct  integration. 

Other  examples  of  direct  integration  are  given  in  reference  20,  which 

includes  the  effects  of  electron -ion  collisions,  and  the  Hall  reduction  in 

conductivity.  The  resulting  expressions  are  generally  more  complex  than 

(14.207).  The  constant  velocity  case  with  thermal  ionization  and  constant 

24 

collision  cross  section  was  Hcst  integrated  by  Brocher, 

(iii)  Generator  Efficiency 

The  overall  generator  adiabatic  efficiency  is  defined  as 


V  = 


T  -  T 
o  3 

T  -  T„' 
o  3 


(14.210) 


where  it  is  assumed  that  the  specific  heat  of  the  gas  is  constant,  see 
Figure  14.28.  Now, 


o  2  1  o 
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path  0-1  occurs  in  the  nozzle,  and  is  close  to  isentropic.  The 
path  2-3  occurs  in  the  diffusor;  this  is  actually  not  isentropic  since  all 
diffusors  have  an  entropy  loss.  However,  for  the  present  purposes  this 
path  will  be  considered  isentropicv  the  effect  of  diffusor  losses  is  given  in 
re-f.  37  .  Thus,  for  paths  0-1  and  2-3,  the  temperatures  are  related 

to  the  pressures  by  the  isentropic  expressions  so  that  (a)  becomes: 


(b) 


If  K  is  constant  for  the  entire  generator,  the  relation  given  by  (14. 190) 
may  be  used  as  follows: 


(c) 


Substitution  of  (c)  into  (b)  yields: 


(d) 


Next,  it  is  necessary  to  express  the  static  pressure  ratio  across  the 

generator,  p„/p, ,  in  terms  of  the  desired  stagnation  pressure  ratio  p  /p  . 

AX.  «j  o 

We  may  write 
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The  first  and  last  pressure  ratios  are  given  by  the  following  isentroplc 
relations: 


Substitution  of  (f,  g,  h)  into  (e)  yields: 


Unfortunately  this  is  an  implicity  relation  for  p„/p-  in  terms  of  the  desired 
stagnation  pressure  ratio.  Thus  for  specified  values  of  and  K, 

(14.211)  must  be  solved  numerically  to  obtain  (p  /p  ).  The  alternative, 

A  X 

of  course,  is  to  choose  (p  /p  )  and  solve  for  p  /p  . 

M  X  O 
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Finally,  T  '/T  is  given  by  the  isentropic  relation: 

o  O 


(14.212) 


overdU 


With  the  lise  of  (14.210,  14.212)  in  (14.209),  the  efficiency  becomes: 

A 


y-  1 


(14.213) 


where  (P-/p„)  is  related  to  p  /p  by  (14.211).  When  K  becomes  very  close 

1  O  O 

to  unity,  ?)— >1  as  expected.  However,  with  K  less  than  unity,  the  effi¬ 
ciency  may  still  become  close  to  imily  If  the  pressure  ratio  p  /p  is 

O  O 

sufficiently  large.  This  is  the  familiar,  internal  reheat  of  turbomachinery  - 

in  the  case  of  MHD  generator,  the  ohmic  heating  is  passed  downstream  and 

is  available  for  conversion  into  electrical  power.  This  effect  may  also  be 

attributed  to  the  "spreading"  of  the  pressure  lines  on  the  temperature- 

entropy  diagram:  At  very  low  pressures,  a  large  entropy  change  causes 

only  a  small  change  in  temperature.  The  various  efficiencies  are  shown 

in  Figures  14.29,  14.30.  In  Figure  14.29  the  value  of  << 

and  the  local  efficiency  77  and  p>olytropic  efficiency  rj  are  both  equal  to  the 

L  p 

loading  factor  K.  But  as  the  pressure  ratio  p  /p  is  increased,  the  generator 

o  O 

efficiency  becomes  greater  than  K. 

The  situation  is  not  nearly  as  favorable  for  large  values  of  the  Mach 

number;  in  Figure  14.30  the  various  efficiencies  are  shown  for 
2 

1/2  (y“l)7y[j  =  1.  The  polytropic  efficiency  is  now  much  less  than 
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EFFICENCY 


LOADING  FACTOR.K 


Figure  14.29.  Frictionless  Adiabatic  Efficiency,  Cros sed-Field 
MHD  Generator,  Constant  Velocity,  «  1. 
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Figure  14,30.  Frictionless  Adiabatic  Efficiency,  Crossed-Field 
MHD  Generator,  Constant  Velocity, 

^  =  1.  =  0.7 


129 


the  loading  factor;  and  the  overall  generator  efficiency  is  only  slightly 
better  than  the  polytropic  efficiency. 

The  variation  of  the  cross-sectional  area  is  obtained  from  the  con¬ 
tinuity  equation;  thus 


A(xl  ^  ^i  ^  ^1  .  TM 

Aj  p(x)  p(x) 

_  1 


(14.214) 


Generally,  the  smaller  the  inlet  Mach  numer  and  the  larger  the  value  of 

19  ^ 

■y,  the  smaller  the  area  increase.  For  continuous  electrodes,  each  of 
which  is  maintained  at  the  same  potential,  the  interelectrode  spacing  h 
must  remain  constant,  since  u  and  K  are  both  constant.  This  means  that 
all  of  the  area  increase  must  be  in  the  direction  of  the  magnetic  field. 

Another  disadvantage  of  the  constant  velocity  generator  is  that  the  local 
Mach  number  at  the  exit  is  higher  than  at  the  inlet,  due  to  the  decrease  in 
static  temperature.  This  increases  the  problem  of  pressure  recovery  in 
the  exit  diffusor,  especially  if  the  exit  Mach  number  is  near  or  greater 
than  unity.  For  this  reason  it  is  desirable  to  consider  other  flows  in  which 
the  exit  Mach  number  is  equal  to,  or  less  than,  the  inlet  Mach  number. 


14.6b  Constant  Mach  Number  Generator 

If  the  Mach  number  of  the  gas  flow  in  the  generator  is  kept  constant,  then 
as  electrical  power  is  extracted  from  the  flow,  both  the  gas  velocity  and  static 
temperature  of  the  gas  will  decrease.  For  this  case,  the  momentum  and  energy 
equations  become: 

p  uu'  +  p’  =  -  (1  -  K)  Ogjj  uB^  (a) 
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2„2 


pu  (C  T'  +  uu*)  =  -ctK  (1  -  K)  u  B 


(b) 


where  '  indicates  d/d^  .  The  local  Mach  number  is  defined  as: 
^  = 


yRT 


(c) 


Since  the  Mach  number  is  taken  constant,  (c)  may  be  differentiated  with  respect 
to  X,  which  yields: 


(d) 


Substitution  of  (d)  into  (a)  and  (b);  then  division  of  (a)  by  (b),  and  rearrangement 
yields: 


J9L 

P 


XL  Z. 

T  K 


1 

-  + 

y  - 1 


Ml 

2 


(1  -  K) 


which  can  be  immediately  integrated  to  obtain  the  following  first  integral  of 
(a,  b): 


-E. 

P' 


1+  ^  (1  -K){y-1)^^ 


(14.215) 


Since  the  poly  tropic  efficiency  is  given  by: 

_ K _ 

1  +  ^  (1  -  K)  (y-1)^^  ’ 
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then  (14.215)  can  be  rewritten  simply  as: 


P' 


(14.216) 


Since  both  the  Mach  number  and  loading  factor  K  have  been  assumed  constant 

for  this  generator,  then  T]  is  also  constant. 

P 

To  obtain  the  variation  of  temperature  or  pressure  in  the  direction  of 

flow,  (14.216)  is  substituted  into  (a)  and  integrated,  which  yields: 

A 


(T-1)  % 

y- 1 

P  > 


2 


1 

2 


(1  -  K)  Uj  t 


(12.217) 


Equation  (14.217)  gives  the  length  in  imits  of  however,  (a)  can  be  integrated 
directly  for  the  case  of  variable  conductivity  of  the  form  given  in  (14.203),  since 
the  pressure  can  be  expressed  in  terms  of  the  temperature  by  means  of  (14.216). 

Because  the  Mach  number  is  constant,  the  stagnation  temperature  ratio  is 
equal  to  the  static  temperature  ratio.  Thus,  the  overall  adiabatic  efficiency 
(neglecting  friction)  is  given  by  the  following  expression: 


T]  (T-i)/r 


(14.218) 
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1  2 

which  is  identical  to  that  for  the  constant  velocity  generator  when  r  (7  -  l)/7L 

A 

<<  1.  hi  fact,  for  moderate  pressure  ratios  the  performance  of  the  constant 
Mach  number  generator  is  identical  to  the  constant  velocity  generator ,  even  for  lovv 
supersonic  Mach  numbers 

14.6c.  Constant  Temperature  Generator 

In  the  constant  temperature  linear  MHD  generator ,  all  of  the  electrical 
energy  is  extracted  from  the  kinetic  energy  of  the  flow.  The  maximum  energy 
is  extracted  when  the  exhaust  velocity  is  very  small;  the  fraction  of  energy 
removed  from  the  gas  is  then: 


T  -T, 
o  1 

T 


o 


1 


-1 

2 


'TfL 


2 

1 

2 

1 


(a) 


Thus,  to  obtain  an  appreciable  energy  extraction  from  the  flow,  the  flow  Mach 
number  at  the  inlet  should  be  greater  than  unity,  depending  on  the  value  of  y. 
The  constant  temperature  generator  has,  however,  the  advantage  that  the  exit 
Mach  number  is  always  less  than  the  inlet  Mach  number,  thus  reducing  the 
dlffusor  losses,  hi  fact,  the  theory  indicates  a  throatless  transition  from 
supersonic  to  subsonic  flow  may  be  possible,  although  this  has  not  yet  been 
proven. 

The  quasi  one -dimensional  equations  for  constant  temperature  flow  are: 

Momentum 


RT 


uu'  +  p' 


-K)uB 


(a) 


Energy 


RT 


uu* 


K(1  -  K)  uB 


(b) 
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Subtraction  of  (b)  from  (a)  yields: 


P'  =  -agjf(l-K)W,  (c) 

then  division  of  (b)  by  (c)  yields: 

j!U_  du/dx  ^  _K...  , 

RT  dp/dx  1  -  K  '  ' 

which  may  then  be  integrated  to  obtain  the  first  integral  by  noting  that  (du/ dx)/ 
(dp/dx)  =  (du/dp).  The  result  is 


£_ 

Pi 


e 


1-K 

2KRT^ 


(14.219) 


To  determine  the  variation  of  velocity  with  the  energy  equation  (b)  may  be 
Integrated  by  use  of  p  =  p/RT^,  and  substitution  of  (14.219).  The  resulting 
expression  is 


[/K-l)r7??jV2K]  f 


1  -K 
2K 


1-K 

2K 


=  (1-K) 


‘^eff 


Uj  B^4/p^ 


(14.220) 


If  the  electrical  conductivity  varies  in  a  manner  given  by  (14.203),  the 

integration  may  be  easily  performed  to  obtain  the  variation  of  the  velocity  with 
25 

the  physical  leng^. 
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The  adiabatic  efficiency  is  given  by 


T?  = 


C  (T 
P  o 


1  -  (u//u^^ 

(14.221) 


Now  the  stagnation  pressure  ratio  is  given  by: 


(e) 


The  first  and  last  factors  in  (e)  are  given  by  the  isentropic  relations,  so  that  (e) 
becomes 


y 

P„  Pi 

Use  of  (14.219)  for  p_/p,  and  ^lft[  ^  =  (u  ^/u  ^)^  ^  in  (f),  the  stagnation 

M  X  A  A  X  X 

pressure  ratio  is  given  by: 


y 


(14.222) 
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Thus,  with  the  Inlet  Mach  number  and  stagnation  pressure  ratio  chosen,  (14.222) 
is  used  to  obtain  the  velocity  ratio  u  /u  ,  which  is  then  substituted  into  (14. 221) 

Ca  X 

to  obtain  the  efficiency.  Figure  14.31  shows  the  results  of  such  a  calculation, 
for 

y-1 

(P3/P0)  =  0-5  and  I  =  1, 

so  that  the  isentropic  velocity  ratio,  u  /u  for  K  =  1,  is  zero  according  to 
(14.222).  Note  that  the  efficiencies  are  quite  close  to  K. 

The  variation  of  efficiency  for  K  =  0. 9  with  stagnation  pressure  ratio  is 

shown  in  Figure  14.32.  It  may  be  seen  that  large  pressure  decrease  the 

A 

efficiency  markedly.  This  is  due  to  the  high  inlet  Mach  number  associated 
with  large  pressure  ratios,  which  causes  a  large  decrease  in  the  polytropic 
efficiency.  Thus,  the  constant  temperature  magnetohydrodynamic  generator 
is  not  suitable  for  very  large  pressure  ratios. 

14. 6d  Constant  Pressure  Generator 

As  can  be  seen  from  the  temperature -entropy  diagram,  Figiure  14.28, 
if  the  local  static  pressure  is  kept  constant,  the  local  static  temperature  will 
increase.  Thus,  only  a  part  of  the  initial  kinetic  energy  is  converted  into 
electrical  energy;  the  remainder  is  converted  back  into  thermal  energy  of  the 
gas .  This  may  be  of  some  advantage  if  the  gas  is  thermally  ionized  since  the 
electrical  conductivity  will  increase  in  the  downstream  direction.  Also,  by 
converting  the  kinetic  energy  back  into  thermal  energy,  the  generator  acts  as 
a  diffusor. 

For  this  generator,  the  momentum  and  energy  equations  become 
respectively; 

p  uu'  =  -  (1  -  K)  uB^  (a) 

P  (CpT'  +  uu')  =  -  (1  -K)  uB^  (b) 
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K 

Figure  14.  31.  Variation  of  Generator  Efficiency  with 
Loading  Factor,  Constant  Temperature 
Generator,  For  (p  =  0.5, 

1/2  (y-1)  uf  =  1. 
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Figure  14.  32.  Variation  of  Constant  Tenaperature 
Generator  Efficiency  with  Pressure 
Ratio  for  K  =  0.  9. 
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Subtraction  of  (a)  from  (b)  yields: 


pC  T» 
P 


(C) 


and  division  of  (a)  by  (c)  yields: 

u  du/dx  _  1 

C  dT/dx  “  “  1  -  K 
P 


(d) 


which  may  be  integrated  immediately  to  obtain  the  following  first  integral: 


1.2  2,  ^ -1  ^  ^  , 
-  (u  -  )  =  -  (1  -K)  Cp  (T  -  T^) 


(14.223) 


To  obtain  the  variation  of  velocity  with  distance,  use  the  perfect  gas  law 
to  eliminate  the  density  from  equaticm  (c),  and  then  eliminate  temperature  with 
the  use  of  (14.223),  as  follows: 


d(u/Uj) 


(y-1) 

2ypj 


^eff^a-K)^  u^dC 


(e) 


where: 


2 


m 


1  + 


_ 2 _ 


Integration  of  (e)  then  yields: 


_1. 

m 


2_ 

m 


-  tanh 


-(y-1) 

2yp^ 


(1-K)^u^4 
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(14.224) 


The  physical  length  can  also  be  obtained  for  variable  electrical  conductivity  by 

expressing  (e)  in  terms  of  temperature  instead  of  velocity,  since  the  conductivity 

25 

will  also  depend  only  on  temperature  because  the  pressure  is  constant.  The 
adiabatic  efficiency  (neglecting  friction)  is  given  by: 


V 


T  -  T„ 
o  3 

^  T  -  T„' 
o  3 


Now, 


(f) 


!3  _  ^  ^ 

Po  '  P2  Po 


y-1 


1  + 


y-i 


te) 


since  p^  =  p^.  Rearrangement  of  (g)  yields'^  2  i*'  terms  of  the  stagnation 
pressure  tatio; 


% 


& 


y 


(1+ 


(h) 


2 


The  temperature  ratio  can  be  expressed  in  terms  of  the  Mach  number  as 
follows: 


T 

o 


T  T 

1  o 


2 


(i) 
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From  (14.223), 


1  +  ll-KL(y-l)oy|  2 
2 


2 

^2 


(j) 


Equation  (j)  may  be  solved  for  T  /T  as  follows: 

A  J. 


T,  1. 

Tj  ■  1  ^  Il-K)(y-1)^  2  ® 

2  2 

Substitution  of  (k)  into  (i),  and  substitution  of  from  (h)  yields: 

^  1  +  (1-K) 

I  y 

KMl-K) 

(14.225) 

Substitution  of  (14.225)  into  f  yields  the  desired  expression  for  efficiency,  in 
terms  of  ,  y,  K,  and  the  stagnation  pressure  ratio  p  /p  . 

d.  o  O 

14. 6e  Constant  Cross-Sectional  Area 

For  the  purposes  of  performing  experiments  on  magnetohydrodynamic 
elecfrical  power  generation,  the  most  convenient  channel  configuration  is  a 
straight,  rectangular  duct  of  constant  width  and  height,  and  hence  constant  area. 
The  performance  will  then  depend  on  whether  the  electrodes  are  continuous  or 
segmented.  If  continuous,  the  voltage  difference  between  the  two  electrodes 
will  be  constant,  and  hence  the  electric  field  will  be  a  constant.  This  causes 
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the  loading  factor  K  to  vary  along  the  length  of  the  generator,  so  that  the  loco.  1 
conversion  efficiency  varies.  In  fact,  it  is  possible  to  have  one  part  of  the 
duct  acting  as  a  generator  while  some  other  part  is  acting  as  an  accelerator. 

On  the  other  hand,  if  the  electrodes  are  segmented,  then  the  voltage 
difference  between  pairs  of  electrodes  need  not  be  constant,  and  can  be  made 
proportional  to  the  local  velocity  by  setting  the  ratio  of  the  external  load  resis¬ 
tance  per  unit  length  to  the  internal  resistance  of  the  generator  per  unit  length 
equal  to  some  constant. 

For  either  configuration,  the  flow  Mach  number  tends  toward  !•  If 
the  flow  was  initially  subsonic,  the  temperature  decreases  while  the  velocity 
increases.  Thus,  additional  thermal  energy  is  converted  into  kinetic  energy 
which  is  not  recoverable  in  the  generator.  If  the  flow  was  initially  supersonic, 
the  velocity  will  decrease  and  temperature  will  increase,  so  that  some  thermal 
energy  cannot  be  converted.  For  these  reasons,  the  constant  cross-sectional 
area  generator  does  not  appear  to  be  practical.  In  fact,  in  multiple  stage  gas 
or  steam  turbines  the  flow  cross-sectional  area  is  always  increased  in  the 
downstream  direction  in  order  to  limit  the  axial  velocity.  It  is  found  that  for 
all  of  the  previous  cases  in  this  chapter,  the  cross-sectional  area  increases  in 
the  downstream  direction  also.  However,  because  the  constant  cross  section 
generator  is  simpler  to  design  and  construct  than  the  others ,  the  theory  is 
included  here. 

The  basic  equations  are: 


Conservation: 


u 


PfUi 


Momentum 


pu 


du  ^  dp 
dx  dx 


%ff  B  (E  -  uB) 


(a) 


(b) 
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Energy 


S  [V"  i 

where  E  is  the  constant  transverse  electric  field,  taken  equal  to  its  value  at 
the  inlet.  It  is  convenient  to  reference  this  value  of  the  electric  field  to  the 
inlet  velocity  ^oUows,; 

E  =  u^  B 


Thus  is  a  constant  for  the  generator.  This  definition  is  different  than  the 
one  used  in  previous  sections ,  where  the  electric  field  was  referenced  to  the 
local  velocity.  It  is  also  convenient  to  let  C^T  =  (p/p)f  and  to 

make  the  equations  non-dimensional  by  letting 


U  =  u/u^ 
P  = 

X  =  x/h 


I 


Then  equations  (b,c)  become: 


> 


dU  ^  dP 
dX  dX 


I  (K^  -  U) 


(d) 


(14.226) 


d_ 

dX 


Kj  I(K^  -  U) 


(14.227) 
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26,27 

Equations  (14.226,  14.227)  have  been  integrated  by  first  eliminating 
I(K  -  U)  between  them: 


This  is  integrated  immediately  as  follows: 

PU  +  =  K  U  +  PK.  +  constant  (e) 

f  1.  A  J.  1 

2  2 

At  the  inlet,  U  =  1,  and  P  =  )  ;  which  allows  the  constant 

in  (e)  to  be  evaluated.  Rearrangement  of  (e)  then  yields  the  following  first 
integral: 


To  obtain  the  variation  of  velocity  with  distance,  the  left  side  of  (14.227)  is 
expanded,  and  (14.226)  and  (14.228)  are  used  to  elinoinate  (dP/dX)  and  P, 
respectively.  Rearrangement  yields: 
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where: 


6  = 


1  + 


1 


(g) 


If  the  effective  conductivity  is  assumed  to  be  constant,  then  the  interaction 
parameter  I  is  also  constant,  and  (f)  can  be  integrated  to  become: 


+  liyL-Jl  _  y(y-i)  KeKn 


^  2  (V  - 1)^ 

^  j‘^1  2y  (>  +  !)  + Kj  6X7-1);^  f,n 


U 


K 


1  -■ 


y  .  u 


y  -1  K, 


JL 


(>  -l)Kj  J 


[2  K  1 

r  1 

LT:7«-  2  O'*'*] 

- 1 

I 

>  -1  K, 


1  - 


(y  -i)K 


1 


=  Kj  IX 


(14.229) 


The  length  required  to  reach  O/l  „  =  1  has  been  calculated  in  reference  27 
and  is  shown  in  Figure  14.33.  It  is  seen  that  the  higher  the  initial  load  factor 

or  the  lower  the  inlet  Mach  number,  the  longer  the  generator  length. 

nu-iTibt*' 

Also,  high  loading  factors  or  low  inlet  Mach  lead  to  the  largest  electrical 

A 

energy  extraction  from  the  flow,  as  is  shown  in  Figure  14.34.  On  the  other 

hand,  the  largest  average  power  density  is  obtained  when  K  %  0.5,  for  any 

27  ^ 

subsonic  inlet  Mach  number. 

For  the  second  case,  in  which  segmented  electrodes  are  used  so  that  the 

local  transverse  electric  field  is  a  fixed  fraction  of  the  local  value  of  uB,  the 

20 

corresponding  first  integral  is  given  by: 


AH/H 


Figure  14.  33.  Length  to  Reach  Mach  !•  0  at  Outlet 
(From  Ref.  27) 


Figure  14.  34.  Energy  Extraction  as  a  Function  of  Inlet 
Mach  Number  (From  Ref.  27) 


14b 


(14.230) 


-■/tKl 


K  -  1 


_2JL 

y  -1 


-  K 


+ 


+ 


1  -  K 


2y 

y  -1 


K(y-i) 


and  the  variation  with  length  is  given  by: 


(1-K)ae«  B  ^ 


PlUi 


2y  -(y-iL  _u 

2y  -K{y-1)  ^ 


+ 


y 

2y  -K(y  -1) 


+ 


1  -  K 


_22_ 

y-1 


-  K 


2y-K(y-i)-| 

y-K(y-i) 


(14.231) 
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14.7  ALTERNATINQ  CURRENT  POWER  GENERATION 


For  many  applications  it  is  desirable  to  have  the  output  power  of  the 
MHD  generator  as  alternating  current  rather  than  direct  current;  in  addition, 
if  ac  induction  generators  could  be  made,  the  electrode  problem  would  be 
eliminated.  However,  there  are  two  problems  associated  with  induction 
generators:  first,  it  appears  that  they  would  have  a  poor  power  factor;  and 
second,  if  the  power  factor  could  be  increased,  the  power  density  in  the  gen¬ 
erator  would  be  reduced.  *  ’ 


The  first  effect,  of  poor  power  factor,  can  be  estimated  from  the  following 

simple  argument:  The  simplest  method  of  generating  ac  power  is  to  alternate 

the  imposed  magnetic  field  at  some  frequency  f.  Now  the  magnetic  field  energy 
2  , 

per  unit  volume  is  B  thus,  the  rate  of  circulation  of  magnetic  field 

energy  per  unit  volume,  P  ,  is: 

O 


P 


B 


(a) 


while  the  generated  electrical  power  per  unit  volume  is  approximately: 


P  ^  K(l-K) 


(b) 


Thus ,  the  ratio  of  the  circulating  magnetic  field  power  to  the  generated  power 
is 


K(l-K) 


(c) 


Next,  a  characteristic  length  L  may  be  defined  as  the  distance  the  gas  will 
move  during  one  period: 


L  = 


u 

f 


(d) 
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so  that  (c)  becomes; 


^  _ 1 _  ^  1 

P  K(l-K)^^uaL  K(1-K)R^ 

where  R  is  the  magnetic  Re3molds  number.  This  ratio  may  be  estimated  with 
in  2 

some  typical  values;  With  a  =  100  Mho/Meter,  u  =  10  M/sec,  and  f  =  60 

cycles/sec,  R^^^  --  2.  The  maximum  value  of  K(1  -K)  is  0.25.  Thus,  the  ratio 

of  circulating  magnetic  power  to  generated  power  is  about  2,  which  means  that 

the  power  factor  is  poor.  In  addition,  the  circuit  used  for  circulating  the 

current  for  the  magnetic  field  must  have  a  low  ohmic  resistance;  otherwise 

this  circuit  will  consume  all  of  the  generated  power. 

14.7a  Induction  MHD  Generator 

The  electrodeless  induction -type  MHD  generator  operates  by  means  of  a 
traveling,  alternating  magnetic  field,  which  induces  currents  in  the  gas.  hi 
turn,  this  induced  current  creates  an  induced  magnetic  field  which  cuts 
the  field  windings  of  the  magnet.  If  the  gas  velocity  is  greater  than  the  velocity 
of  the  traveling  magnetic  field,  the  induced  currents  will  be  in  such  a  direction 
that  the  Lorentz  force  retards  the  flow;  and  in  addition,  there  will  be  power 
output  from  the  field  windings. 

A  very  simplified  analysis  of  the  induction  generator  to  illus¬ 
trate  its  performance,  see  Figure  14.35.  more  detailed  analyses 

have  been  made.  It  is  assumed  that  the  generator  channel  is  infinite  in  the 
y-direction.  The  field  coils  are  also  in  the  y-direction,  and  are  so  excited  that 
the  z -component  of  the  magnetic  field  is  approximately  sinusoidal  and  moving 

in  the  downstream  direction  x  with  velocity  u  .  The  gas  is  assumed  to  have  a 

B 

constant  scalar  conductivity,  to  be  incompressible,  and  moving  in  the  down¬ 
stream  direction  with  constant  velocity  u.  Also,  it  is  assumed  that  the  induced 
current  loops  close  at  y  =  +  <».  In  practice,  this  may  be  accomplished  by  the 
use  of  an  annular  geometry. 
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In  the  absc'iicc  of  a  conducting  fluid  in  the  generator,  the  imposed  mag¬ 


netic  field  B  is  given  by: 
o  ” 


B 

A 


O 


Re  B* 
o 


(x  -  U^t) 


(a) 


where  X  is  the  wavelength  of  the  magnetic  field.  In  the  presence  of  a  conducting 
fluid  there  will  also  be  an  induced  magnetic  field  ,  caused  by  the  currents  in 
the  fluid  which  may  be  represented  by:  ^ 


Re  B! 


1 


i  2  n 

T" 


(X  -  Ugt) 


(b) 


such  that 


VxB. 


“o  Ji 


(c) 


WriAtiooel  iVir 

where  j.  is  the  induced  current  in  the  fluid.  It  is  assumed  that  Ihe^magnetic 
field  is  only  in  the  x -direction,  and  the  current  is  only  in  the  y -direction.  Then 
(c)  becomes: 


-^B 

1 


3x 


-  p 


o 


(cl) 


Use  of  (b)  for  B 

T 


yields: 


i  2:: 
A 


B 


z. 

1 


■  ^  •'y. 
'  1 


(e) 
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The  induced  current  may  also  be  obtained  from  Ohms  Law.  First,  the  electric 
field  in  the  y-direction  is  obtained  from  Maxwell's  equation; 


9x  at 


(f) 


Let 


E 

y 


i2  n 

X 


(X  -  Ugt) 


(g) 


B  =  B  +  B 
z  z  z. 

O  1 


B* 


i2n 

X 

e 


(X  -  Ugt) 


(h) 


Substitution  of  (g),  (h)  into  (f)  yields; 


(i) 


Ohms  law  then  yields; 


j  =  a  (E 
Yi  5 


uB  ) 
z 


<"b  -  ®z 


=  -  (1  -  S)  CT  u  (B  +  B  ) 
z.  z 


(j) 


where  S  is  the  slip,  S  -  Equating  the  two  expressions  for  the  induced 

current  (e)  and  (j),  one  may  solve  for  the  induced  magnetic  field; 


B 


B 


R 


m 


(k) 


1  SZ 


z. 

1 


1 


where  R  is  a  magnetic  Reynolds  number  defined  by; 
m 

R  =u  auX(l-  S)/2  rr 
m 


The  total  magnetic  field  is  given  by: 


B  =  B  +  B 
z  z.  z 
1  o 


B 

z 

_ o 


1  +  iR 

m 


(1) 


(m) 


and,  from  (i): 


E 

y 


B 

B  z 


1  +  iR 

m 


(n) 


Now  the  current  in  the  field  windings  may  be  related  to  B  by  means  of  Stokes 

z 

theorem  Va>'  iVte  cifeo-  4he  Itine  • 


B  (x)  -  B 
z  z 

L  o 


(x  +  Ax) 


b 


2u  J  Ax 
o  y 


(o) 


where  J  is  the  field  current  per  unit  length  in  the  x-direction  on  each  side  of 

y 

the  generator.  Equation  (o)  may  be  written  as: 


2  u  J 
o  y 

b 


(P) 


Thus, 


J 
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i  2  rrb 
2fx  X 
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(r) 
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Now  the  complex  power  per  unit  length  and  width  delivered  to  the  upper  or 
lower  field  windings  is  given  by  P*  =  ^  J^,  where  is  the  complex  con¬ 

jugate  of  E^.  Thus,  the  total  power  per  unit  volume  is; 


p*  =  — 

2 


E  J 

y  y 


n  B  u„o 
o  B 

M  ^ 


R  -1 
m 


1  +  R 


m 


(14.232) 


where  equations  (n,  r)  have  been  used.  The  real  power  is  given  by  the  real 

part  of  (14.232)  while  the  reactive  power  is  given  by  the  imaginery  part  of 

(14.232).  Equation  (14.232)  thus  clearly  illustrates  that  the  reactive  power  is 

larger  than  the  actual  power  when  R^  <1.  To  increase  the  power  factor, 

must  be  increased,  but  then  the  power  density  will  decrease  for  R^  >  1. 

Equation  (14.232)  may  also  be  compared  to  the  direct  current  generator  by 

substitution  of  (1)  for  R  in  the  numerator; 

m 


Re  P  =  - 


S(1  -  S) 


1  +  R 


m 


Thus,  the  slip  velocity  S  replaces  the  load  factor  K.  The  factor  of  ^  is  caused 
b}’  the  sinusoidal  variation  of  power  output,  as  is  usually  associated  with  alter¬ 
nating  current  power.  Equation  (14.232)  clearly  shows  that  increasing  the 
magnetic  Reynolds  number  will  decrease  the  power  density.  However,  it  is 
unlikely  that  magnetic  Reynolds  numbers  in  excess  of  unity  can  be  achieved,  so 
that  the  main  problem  is  the  power  factor.  From  (14.232),  it  is  seen  that  the 

phase  angle  between  the  voltage  and  current  is  given  by  0  =  tan  ^  R  -1;  more 

m 

28 

detailed  analyses  have  verified  this  poor  power  factor. 
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14.8  MHD  POWER  CYCLE  ANALYSIS 


In  the  previous  sections ,  methods  were  presented  for  calculating  the 
performance  and  efficiency  of  various  magnetohydrod3mamic  electric  power 
generators.  These  methods  allow  for  the  calculation  of  the  decrease  in  total 
enthalpy  in  the  generator,  the  pressure  drop,  and  the  equivalent  turbine 
efficiency.  In  the  present  section,  it  is  assumed  that  this  performance  can  be 
calculated,  and  it  is  now  desired  to  determine  the  overall  efficiency  of  a  gen¬ 
erator  system. 

There  are  three  general  types  of  MHD  generator  systems.  In  the  simplest 
system,  liquid  or  solid^fuel  and  liquid  oxidizer  are  burned,  seed  added,  and  the 
resulting  gas  mixture  passed  through  the  generator  with  no  further  recovery  of 
the  remaining  thermal  energy  or  seed  material  in  the  exhaust  gases .  Such  a 
system  is  feasible  only  for  small  or  short  duration  generators.  Since  the 
amount  of  energy  required  to  pump  the  liquid  (or  solid)  fuel  and  oxidizer  is 
small,  all  of  the  generated  electrical  power,  is  available  for  the  load  except  that 
required  to  energize  the  magnet,  if  it  is  an  electromagnet.  This  system  will 
be  discussed  further  in  section  14.8a. 

For  the  generation  of  bulk  power,  fossil  fuel  must  be  used  for  economic 
reasons,  and  the  oxidizer  must  be  air,  although  some  enrichment  with  addi¬ 
tional  oxygen  may  be  possible .  In  order  for  the  combustion  temperature  to  be 
high  enough  for  thermal  ionization,  preheating  of  the  incoming  air  is  necessary. 
The  exhaust  temperature  of  the  gas  from  the  generator  is  determined  by  the 
decrease  in  the  ionization  of  the  gas;  but  because  combustion  gases  must  be 
thermally  ionized,  the  exhaust  temperature  is  still  very  high.  For  economic 
operation  further  extraction  of  this  energy  is  necessary,  which  may  be  accom¬ 
plished  by  using  the  exhaust  gases  to  generate  steam  power.  These  systems 
are  described  further  in  Section  14.8b. 

Finally,  the  gas  which  passes  through  the  generator  may  be  in  a  com¬ 
pletely  closed  cycle;  that  is,  the  gas  (or  vapor)  is  heated  in  a  heat  exchanger, 
expanded  through  a  nozzle,  passed  through  the  generator,  cooled,  and 
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recompressed  (or  pumped)  back  through  the  heat  exchanger.  For  such  a  cycle, 
the  maximum  temperature  of  the  gas  is  set  by  the  allowable  working  tempera¬ 
ture  of  the  materials  in  the  heat  exchanger.  At  the  present  state  of  development 
of  refractory  materials ,  it  is  highly  unlikely  that  sufficiently  high  temperatures 
can  be  sustained  over  long  periods  of  time  to  permit  thermal  ionization;  thus, 
some  form  of  non-thermal  ionization  is  necessary. 


At  present,  the  main  use  for  such  closed  cycles  appears  to  be  for  vehicle 
power.  The  most  interesting  potential  application  for  such  cycles  is  for  the 
generation  of  electrical  power  in  space  where  the  heat  rejction  must  be  by 
radiation.  This  is  described  in  Section  14. 8. c. 

14.8a  Simple  Open  Cycle 

In  the  simple  open  cycle  fuel  and  oxidizer,  without  preheating,  are  burned, 
and  "seed"  material  is  added.  The  combustion  products  are  accelerated  by  a 
nozzle  and  passed  into  the  generator  and  exhausted  into  the  atmosphere.  Since 
seeded  combustion  gases  are  not  sufficiently  conducting  below  about  4000°F  at 
a  pressure  of  one  atmosphere,  the  flame  temperature  must  exceed  4000^F  in 
order  to  extract  any  electrical  power.  Since  the  temperature  of  most  flames 
with  air  is  below  4000°F,  it  is  necessary  to  use  either  pure  oxygen  or  a 
chemical  oxidizer.  Even  if  pure  oxygen  is  used  the  flame  temperature  with 
ordinary  hydrocarbons  is  not  especially  high  due  to  dissociation  of  the  products 
•  of  combustion,  see  Table  14.4. 

The  overall  efficiency  for  this  type  open  cycle  generator  can  be  expressed 
in  terms  of  the  enthalpy  difference  between  the  flame  temperature  and  ambient: 

Net  work/unit  mass  of  gas 

’’t  '  H  -  h 

o  a 

where  is  the  total  enthalpy  of  the  products  of  combustion  at  the  combustor 
pressure,  and  h  is  the  sensible  enthalpy  of  the  same  gas  if  cooled  to  ambient 

3. 
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pressure  and  temperature.  The  net  work  is  the  difference  between  the  gener¬ 
ator  inlet  total  enthalpy  and  the  exit  total  enthalpy,  less  any  heat  transfer,  and 
power  for  auxiliary  equipment.  Thus, 


H  -h„ 
o  2 


/2  -  Q  -  P 
2  aux 


H  -h 
o  a 


The  thermal  efficiency  is  therefore  limited  by  the  minimum  temperature 
at  which  there  is  sufficient  ionization,  and  by  the  exhaust  velocity.  If  the 
exhaust  velocity  is  reduced  by  adjusting  the  cross-sectional  area,  then  addi¬ 
tional  energy  can  be  extracted  from  the  gas  before  T  .  is  reached;  however, 
decreasing  the  gas  velocity  generally  increases  the  generator  length.  For  this 
reason,  it  is  desirable  to  use  the  higher  flame  temperatures. 


TABLE  14.4.  FLAME  TEMPERATURES  OF  FUELS  BURNED 
WITH  LIQUID  OXYGEN  AT  20  ATM.  PRESSURE 

Fuel  Flame  Temperature.  °F 


*Ethyl  Alcohol  5250 

*Kerosene  557  0 

♦Hydrogen  ~  5400 

♦Methane  5040 

Cyanogen  ~  8200 


♦From  Rocket  Propulsion  Elements  by  G.P.  Sutton,  2nd  Edition,  John  Wiley 
and  Sons,  Inc.,  1956,  N.  Y. 


Also  listed  in  Table  14.4  is  the  flame  temperature  with  Cyanogen,  C^N^.  The 
high  flame  temperature  results  not  so  much  from  a  larger  heat  of  combustion, 
but  from  the  fact  that  the  products  of  comljustion  do  not  dissociate  as  in  hydro¬ 
carbon  flames;  that  is: 
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22 

Detailed  calculations  of  the  performance  of  an  ojqt -cyanogen  generator 
indicate  that  over  36%  of  the  input  chemical  energy  can  be  converted  into  ter¬ 
minal  power.  This  is  to  be  compared  to  about  15%  for  hydrocarbon  fuels. 

With  an  initial  total  pressure  of  one  atmosphere,  K  =  0.5 ,  =  0. 9,  and 

segmented  electrodes,  only  10  cm.  of  length  is  required  to  obtain  the  overall 

2 

efficiency.  The  magnetic  field  was  2  webers/M  at  the  inlet  and  0.45  at  the 

exit.  Because  of  this  short  length,  it  is  necessary  to  have  a  small  cross- 

sectional  area  in  order  to  have  quasi-one  dimensional  flow.  For  an  inlet  area 
2  2 

of  1  cm  ,  the  exit  area  is  63  cm  ,  and  the  total  generated  power  was  10. 87 
kilowatt.  No  magnet  power  was  required  because  it  was  assumed  that  the  mag¬ 
netic  field  was  supplied  by  a  permanent  magnet.  For  larger  total  powers, 
larger  inlet  cross  sections  are  required,  and  the  length  must  be  increased,  but 
this  can  be  easily  accomplished  by  reducing  the  magnetic  field. 

In  these  calculations,  it  was  assumed  that  the  seed  material  was  10%  pure 
cesium.  Because  of  the  high  cost  of  cesium,  this  system  is  attractive  only  for 
small  units  and  short  duty  cycles.  A  cheaper  source  of  cesium  is  pollucite  ore 
which  contains  27%  cesium.  Certain  potassium  compounds  such  as  the  hydroxide 
or  carbonate  are  also  relatively  inexpensive.  The  halides  are  generally  not  as 
suitable  for  seed  materials  due  to  the  formation  of  negative  chlorine  atoms 
which  decrease  the  electron  density. 

(i)  M  agnet  Power 

A  fraction  of  the  output  power  may  be  used  to  energize  an  electro¬ 
magnetic,  instead  of  using  a  permanent  magnet.  To  determine  the  con¬ 
ditions  which  must  be  met  for  self -energizing  the  field  coil,  we  may 
examine  the  equation  for  the  field  current  to  the  electromagnet,  I: 


L  ^  +  m  =  huB 
dt 


(a) 


where  R  is  the  sum  of  the  internal  and  external  resistance.  If  it  is 
assumed  that  the  magnet  power  is  obtained  from  some  fraction  o  of  the 
total  length  of  generator  ^  ,  then 


2N-t 


(b) 


R 


ahal 


where  cr  and  a  are  the  gas  and  wire  electrical  conductivities,  respec- 

w 

tively,  N  is  the  number  of  turns  of  the  field  coil,  h  is  the  transverse 
distance  between  electrodes,  and  b  is  the  transverse  distance  in  the 
direction  of  the  magnetic  field,  see  Figure  14.36. 

The  magnetic  field  for  an  air  core  and  iron  core  magnets  are  given 
by  respectively: 

Air  Core 


B 


u  NI 
^o 


(C) 


Iron  Core 


B 


u  NI 


(d) 


where  is  the  average  length  of  the  flux  path  for  the  air  core  magnet, 
ajid  b'  is  the  iron  magnet  gap.  Combining  (a  -  d),  it  is  seen  that  the 
requirement  for  self-exitation  vvhiob  lO  JI/ itt.  ke.t.c,ytes- 


Air  Core 


'U  u  b  f  a  N 
o 


H  A  uh 
W  o  W 


<  1 


(e) 


Iron  Core 


^  ub  C  a  N 
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(f) 
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AIR  CORE  MAGNET 


Let  us  consider  the  air  core  magnet  first.  Equation  (e)  may  be  written  as 


1  1 

R  •  R  " 
m  m 


<  1 


(g) 


where  R^'  and  are  modified  magnetic  Reynolds  numbers  defined  by:* 


^  ’ 


(h) 


R  M  = 

m 


\  > 


(i) 


•fkeldL 

Now,  the  value  of  the  desired  magnetic  affects  both  (h  and  i),  because 

A 

increasing  the  magnetic  field  decreases  the  length  of  the  generator.  It 
is  easily  shown  that  the  left  side  of  (g)  is  a  minimum  with  respect  to  i 
when  R^*  =  R^".  Equating  (i),  (h)  and  solving  for  I,  one  obtains 


I  = 

w  b  tW 


2a  Nb 
{c^)  ^ 


R  '  >  2 
m 


0) 

(k) 


Substitution  of  (h,  j)  into(14  and  rearrangement  yields: 

Air  Core 


AwN 


W 


9  9 

ou  u  a 


(1) 


*It  is  assumed  that  ion  slip  is  negligible,  and  that  the  electrodes  are 
segmented. 
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The  left  side  of  (1)  is  just  the  total  cross-sectional  area  of  the  wire, 

bh  is  just  the  average  cross-sectional  area  of  the  channel,  while  ^  ^  is 

the  average  area  enclosed  by  the  magnetic  flux  path.  The  larger  the  cross- 

2 

secticmal  area  of  the  channel,  the  larger  oO  will  be.  Hence,  the  ratio 
(bh/«^^  is  approximately  a  constant,  and  is  independent  of  the  cross- 
sectional  area  of  the  generator. 

We  may  next  determine  the  effect  of  generator  size  on  the  required 

total  cross  section  of  wire,  for  a  given  gas  velocity,  and  fraction  of 

power  directed  to  the  magnet  a.  Since  the  right  side  of  (1)  is  a  constant, 

independent  of  cross-sectional  area,  the  total  cross  section  of  wire  is 

also  constant.  Thus,  the  dimensions  of  the  field  windings  are  independent 

of  the  generator  power  output,  e.g. ,  generator  cross-sectional  area. 

Thus,  air -core  ms^ets  favor  large  generators.  For  example,  let 

=  10®  Mhos/M;  u  =  10®M/sec,  a  =  10  (bh/<^^)  =  10  For 

hydrocarbon -ojQTgen  flames,  a  ~  30  Mho/M,  and  (1)  predicts  a  minimum 

2 

wire  cross  section  of  0.8  Meters  .  This  large  cross-sectional  area  of 

wire  makes  it  impractical  to  build  generators  with  small  duct  cross- 

sectional  areas.  For  example,  if  we  set  an  arbitrary  Limit  that  the 

cross-sectional  area  of  the  generator  duct  must  be  at  least  (0.01  A  N, 

2  ^ 

then  the  minimum  duct  cross-sectional  area  is  80  cm  . 


On  the  other  hand,  if  cyanogen -oxygen  is  used,  with  a  conductivity  of 

about  ten  times  larger,  the  minimum  wire  cross-sectional  area  is  about 
2 

800  cm  ,  which  makes  it  possible  to  have  a  minimum  duct  area  of  about 

„  2 
8  cm  . 

A  similar  analysis  can  be  made  for  iron  core  magnets.  Then  (1)  is 
replaced  by 


A^N 


2  2 

a/U  u  a  bh 


W 


(m) 
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The  minimum  wire  cross-sectional  area  is  now  much  smaller,  since, 

for  b  h,  b'  is  only  slightly  larger  than  b;  let  us  assume  that  b'  =  2b  =  2hi 

then  the  minimiun  wire  area  is  reduced  from  the  previous  example  by  a 

2  2 

factor  of  25,  which  becomes  320  cm  for  hydrocarbons,  and  32  cm  for 
cyanogen.  The  cross-sectional  area  of  the  iron,  on  the  other  hand,  is 
directly  proportional  to  the  generator  cross-sectional  area,  except  for 
very  small  sizes  in  which  case  sufficient  space  must  be  allowed  for  the 
field  coil  windings.  Thus,  the  weight  of  an  iron  core  magnet  is  approxi¬ 
mately  directly  proportional  to  the  generator  output,  except  for  very 
small  sizes.  This  is  in  contrast  to  the  cross  section  of  the  wire,  which 
is  almost  independent  of  generator  size. 

(ii)  Wall  Losses 

The  two  principal  wall  losses  are  skin  friction  and  heat  transfer. 
Skin  friction  increases  the  required  pressure  ratio  across  the  generator. 

If  the  fuel  and  oxidizer  are  both  liquid,  the  increased  pressure  ratio  does 
not  appreciably  affect  the  net  power  output  because  the  pumping  energy 
for  liquids  is  small.  On  the  other  hand,  if  the  increased  pressure  ratio 
is  obtained  by  increasing  the  inlet  pressure,  then  for  thermal  ionization, 
the  electrical  conductivity  will  decrease  and  the  generator  length  will 
increase. 

The  other  major  wall  loss  is  heat  transfer,  if  it  is  necessary  to 
cool  the  walls.  Heat  removal  reduces  the  enthalpy  available  for  conversion 
into  electricity  before  reaching  the  minimum  temperature  at  which  the 
gas  is  sufficiently  ionized. 

14. 8b  Open  Cycle  with  Recovery 

The  previous  open  cycle  is  clearly  unsuited  for  economic  operation  for 
three  reasons:  the  seed  material  is  not  recovered;  chemical  oxidizers  are 
expensive;  and  finally,  the  thermal  energy  of  the  exhaust  gases  is  not  recov¬ 
ered.  For  the  economic  generation  of  power,  for  example  bulk  power,  these 
three  problems  must  be  considered. 
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The  cheapest  available  oxidizer  is  obviously  air,  although  power  is 

required  to  compress  it  to  the  required  combustion  pressure,  hi  addition, 

comnaon  fossil  fuels  must  be  used  if  the  fuel  cost  is  to  be  minimized.  Now,  it 

is  impossible  to  burn  fossil  fuels,  such  as  natural  gas,  oil,  or  coal,  and  achieve 

sufficiently  high  combustion  temperatures  to  ionize  the  seed  material.  Hence, 

preheating  of  the  inlet  air  is  required  to  increase  the  combustion  temperature; 

29 

that  is,  a  regenerative  cycle  is  required,  see  Figures  14.37,  14.38. 

To  design  such  a  cycle  it  is  first  necessary  to  select  the  generator  exhaust 
velocity  u  ,  exhaust  stagnation  pressure  p  and  exhaust  temperature  T  .  The 

u  O  iu 

exhaust  pressure  should  be  high  enough  to  allow  the  gas  to  pass  through  the 
regenerator  and  any  other  heat -recovery  equipment.  The  exhaust  static  tem¬ 
perature  should  be  high  enough  so  that  there  is  sufficient  thermal  ionization  of 
the  seed  material.  From  the  Mollier  diagram  for  the  combustion  products  and 
seed  material,  aiid  the  diffusor  losses,  T  may  be  determined  next.  It  is  not 
necessary  to  select  the  allowable  temperature  drop  AT  =  T  -  T  through  the 

O  / 

regenerator.  The  smaller  the  value  of  AT,  the  larger  will  be  the  physical  size 
of  the  regenerator  and  vice  versa.  From  the  formulas  of  Section  14.6,  the 
generator  may  be  designed,  but  since  the  combustion  pressure  has  not  yet  been 
selected,  this  must  be  done  parametrically  for  various  combustion  pressures, 
p^'.  The  results  may  be  plotted  as  line  O'  -  3. 

Next,  for  an  air  preheat  temperature  T^,  the  flame  temperature,  including 
the  seed  material,  may  be  calculated  for  various  values  of  the  combustion  pres¬ 
sure  p^'.  Some  typical  calculations  are  shown  in  Figure  14.39. 

The  calculated  combustion  temperatures  may  be  plotted  as  T^"  as  a  func¬ 
tion  of  p^'.  The  intersection  of  the  O'  line  and  the  0"  line  gives  the  desired 

combustion  temperature  and  pressure.  The  compressor  pressure  p  is  then 

6 

selected,  allowing  for  the  required  pressure  drop  through  the  cold  side  of  the 

regenerator.  From  standard  compressor  formulas,  T  is  then  calculated  with 

6 

T  ,  p  selected  at  ambient.  Finally,  the  regenerator  hot-side  exhaust  tern- 

O  u 

perature  is  calculated  from; 

•  • 

”^6  ^7  ~  ^6^  "  "'o  ^3  "  ^4^ 

TVie  ai4erc.ee  between  \  <s  He  ivi«i  flow  ra+e  ^.f 

fvxe\  pluj  Seed  material 
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TEMPERATURE 


ENTROPY 


Fi'jurt-  14,  37,  Res^cnerative  Cycle,  Temperature  -  Entropy  Diagram 


1  6  “) 


DC -AC  Converter 


Open  System  MHD  Cycle 


Figure  14.38.  Regenerative  Cycle,  Flow  Schematic  (After  S.  Way) 


16b 


TO.COMBUSTION  TEMPERATURE,  “K 


Figure  14,  39.  Dependence  of  Connbustion  Temperature  on  Air  Preheat 
Temperature 

Carbon  -  9,7  atm.  pressure,  J.  Fabre  and  J.  Pericart, 

/ 

"Quelques  etudes  preliminaires  en  vue  de  la 
Generation  Magnetohydrodynami<^ue  de  I'energie 
Electrique",  presented  at  Societe  Hydrotechnique 
de  France,  17  March  1961 

Gases  -  1.0  atm.  ,  G.  Ribaud  and  N.  Nauson,  Applications 
des  constantes  et  donnees  thermodynamiques  des 
me  langes  gazeux  aux  temperatures  elevees. 
Publications  Scientifiques  et  Techniques  du 
Ministere  de  Pair,  No.  341  (1958) 
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Now,  because  the  specific  heat  of  combustion  gas  is  larger  than  that  for  air,  the 

regenerator  hot -side  temperature  T  will  normally  be  considerably  higher  than 
A  4 

the  regenerator  cold -side  inlet  temperature,  and  additional  energy  can  be 
recovered  from  the  exhaust  gases,  by  pro^cing  steam  and  using  a  steam  tur¬ 
bine.  It  is  possible  to  couple  the  turbine  directly  to  the  air  compressor,  and  to 
use  any  remaining  shaft  power  to  generate  electrical  power.  The  net  power 
output  is  then  the  direct  current  generated  in  the  MHD  generator,  plus  any 
remaining  steam  power,  less  the  magnet  power  and  auxiliary  power.  The  dc 
power  is  essentially  m  (h  -  h  ),  less  any  energy  lost  by  heat  transfer.  (It  is 

o  O  <j 

possible,  of  course,  to  use  the  MHD  generator  coolant  as  the  steam  cycle  feed 

water  heater.)  With  these  types  of  calculations,  overall  thermal  efficiencies 

29 

as  high  as  56%  have  been  calculated. 

Several  problems  in  the  cycle  are  evident.  First,  for  combustion  products 
and  about  1% molar  potassium  seeding,  T  is  about  4000°F,  and  the  regenerator 

A 

walls  will  then  be  close  to  this  temperature.  At  the  present  time,  economic 
materials  which  will  withstand  this  temperature  for  long  periods  of  time  are 

5 

not  known  (although  some  show  promise  for  future  development).  To  reduce 
the  regenerator  surface  temperature,  several  possibilities  have  been  explored. 
First,  T^  can  be  reduced,  which  will  help  lower  the  regenerator  surface  tem¬ 
perature.  However,  this  reduces  the  combustion  temperature  and  the  power 
generated  in  the  MHD  generator.  An  alternate  to  this  is  to  cool  the  gases  from 
temperature  Tg  by  use  of  a  steam  superheater,  before  passing  the  gas  into  the 

on 

regenerator.  This  also  reduces  the  combustion  temperature  T^. 

A  different  solution  is  to  reduce  T^,  as  explained  above,  but  to  obtain  the 

31 

same  combustion  temperature  by  enriching  the  inlet  air  with  oxygen.  Although 
this  alleviates  the  regenerator  materials  problems,  for  a  given  mass  flow  of 
combustion  products,  additional  fuel  is  required  to  burn  the  additional  oxygen, 
which  reduces  the  thermal  efficiency  somewhat;  in  addition,  additional  invest¬ 
ment  in  oxygen  separation  equipment,  and  additional  auxiliary  power  is  required. 

A  more  novel  remedy  which  has  been  suggested  is  to  use  molten  ash  as 
32 

the  regenerator,  which  is  similar  in  concept  to  the  "rotating  regenerator" 
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under  development  for  gas  turbines.  The  exhaust  gases  from  the  MHD  gener¬ 
ator  are  used  to  heat  ash  until  molten;  the  molten  ash  is  then  sprayed  through 
the  compressed  air  in  a  counter  flow  arrangement.  The  ash  solidifies  as 
granules,  and  is  mechanically  recirculated  to  the  exhaust  gas  heater. 

A  more  practical  method  of  reducing  the  regenerator  temperature  is  to 
use  pollucite  as  the  seed,  in  order  to  utilize  the  lower  ionization  potential  of 
cesium,  and  to  use  the  higher  field  strengths  that  may  become  available  through 
the  development  of  high  field  superconductors.  The  latter  also  has  the  advan¬ 
tage  that  the  magnetic  field  itself  consumes  no  power;  the  cmly  power  c(»sumed 
by  the  magnet  system  is  that  required  to  recondense  liquid  helium  to  keep  the 

windings  superconducting.  With  these  two  modifications,  it  has  been  estimated 
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that  a  reduction  in  temperature  of  as  much  as  700  C  can  6e  expected. 

Figure  14.39  also  shows  that  for  a  given  amount  of  preheat,  the  conibus- 

tion  temperature  varies  with  fuel.  Although  this  is  caused  in  part  by  the  larger 

specific  heat  of  the  combustion  products  containing  water  vapor,  it  is  important 

to  select  a  fuel  which  yields  the  largest  enthalpy  change  h„  -  h  ,  for  a  given 

3  o 

MHD  generator  exhaust  temperature  T-  and  preheat  temperature  T_.  hi  addi- 

it  7 

tion,  the  cost  of  the  fuel  should  be  considered.  On  these  two  bases,  coal  seems 
to  be  the  best  fuel  for  economic  operation  of  an  MHD  generator.  On  the  other 
hand,  the  high  ash  content  of  coal  causes  additional  difficulties.  Although  it  is 
possible  to  remove  as  much  as  90%  of  the  ash  in  a  "cyclone"  burner,  the 
remaining  ash  becomes  mixed  with  the  seed  material,  which  may  be  a  potassium 
compound  or  pollucite.  Since  economic  operation  requires  recovery  of  the  seed 
material,  this  means  that  the  10%  ash  must  be  recovered  as  well,  and  recircu¬ 
lated.  To  keep  the  "seed"  makeup  at  a  minimum,  scrubbers  and  electrostatic 
precipitators  are  required;  but  a  more  important  problem  is  the  large  weight  of 

ash  that  must  be  recirculated  through  the  system  in  order  to  make  use  of  the 

34 

seed  trapped  in  it.  Way  has  estimated  that  the  ash  recirculation  is  about  (me- 
third  the  weight  of  the  coal  being  burned!  Thus,  the  ];Aysical  size  of  the  scrubbers 
and  precipitators  becomes  enormous.  In  addition,  the  molten  ash  will  coat  out 
on  the  inside  of  the  generator  and  will  quite  likely  short  out  the  electrodes. 
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Also,  ash  will  coat  out  on  the  regenerator.  Thus,  the  problems  of  burning  coal 
are  quite  difficult. 

In  addition  to  the  problems  caused  by  ash,  there  are  two  other  important 

problems  associated  with  continuously  operated  generators.  First,  the  potas- 

46 

sium  attacks  insulating  materials  and  makes  them  conducting.  Second,  the 
electrode  materials  -  tungsten,  carbon,  or  silicon  carbide,  are  chemically 
eroded  by  combustion  gases.  This  last  problem  might  be  solved  by  using  con¬ 
sumable  electrodes;  that  is,  feed  the  electrodes  continuously  into  the  generator 

32 

in  a  manner  similar  to  that  used  in  electric  arc  furnaces . 

The  largest  problem  that  is  faced  by  large-scale,  continuous  generation 
of  power  by  fossil  fuel  -  MHD  generators  is  economics .  Although  this  system 
shows  promise  of  increasing  the  overall  thermal  efficiency  to  as  much  as  56% 
as  compared  to  40%  for  conventional  steam  plants,  the  additional  investment  in 
the  magnet,  generator  duct,  compressors,  regenerator,  scrubbers,  precipi¬ 
tators,  dc  to  ac  inverters  if  the  ac  is  required,  causes  the  plant  cost  of  the 

generated  electrical  power  to  be  essentially  the  same  as  that  for  present  day 
30 

steam  plants.  For  this  reason,  there  is  at  present  no  major  impetus  to 
develop  large  fossil  fuel  MHD  generator  plants  in  the  United  States.  However, 
in  other  countries  such  as  Japan,  England,  and  France,  where  the  fuel  cost  is 
much  higher,  and  capital  cost  is  lower  than  the  United  States,  there  may  be 
potential  economic  gain  in  the  use  of  MHD  power  plants.  For  this  reason, 
research  on  the  development  of  fuel-burning  MHD  generators  is  being  actively 
pursued  in  those  countries. 

If  the  cost  of  fossil  fuel  were  to  Increase  in  the  United  States,  then  nuclear 

power  would  become  economic  rather  than  fossil  fuel  -  MHD,  primarily  because 

of  large  successes  in  reducing  the  capital  costs  of  nuclear  plants  in  the  United 
30 

States.  On  the  other  hand,  the  capital  cost  for  nuclear  plants  in  other  coun¬ 
tries  has  remained  so  high  that  even  with  higher  fossil  fuel  costs,  nuclear  power 
is  not  economically  competitive.  This  situation  can  of  course  be  drastically 
altered  by  new  innovations  in  technology. 
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14.8c  Closed  Nuclear  -  MHD  Cycle  for  Space 

A  third  possible  application  for  MHD  power  generation  is  a  closed  cycle 

with  a  nuclear  heat  source.  This  appears  to  be  especially  attractive  for  the 

generation  of  electrical  power  in  space,  where  the  heat  rejection  must  be  by 

radiation,  because  the  heat  rejection  per  unit  area  varies  as  the  fourth  power 

a(i<i 

of  the  surface  temperature,  a  high  temperature  cycle  is  required  to  minimize 

A 

the  radiator  area.  At  present,  h^h  temperature  nuclear  reactors  are  under 
development;  but  the  corresponding  development  of  turbogenerators  is  difficult 
because  of  the  combination  of  high  operating  temperatures,  blade  stresses  in 
the  turbine,  and  the  corrosive  properties  of  some  high  temperature  working 
fluids ,  such  as  liquid  metals .  The  advantage  of  the  MHD  generator  is  that  the 
stress  level  is  reduced  drastically.  On  the  other  hand,  it  is  not  likely  fliat  the 
reactor  temperature  will  be  sufficiently  high  that  thermal  ionization  of  the 
working  fluid  in  the  MHD  generator  will  be  possible ;  thus ,  some  form  of  non- 
thermal  ionization  will  be  required.  Some  other  possibilities  will  be  discussed 
later  in  this  section. 

(i)  Optimum  Heat  Rejection  Temperature 

The  radiator  for  a  space  nuclear  closed  cycle  represents  a  major 
portion  of  the  weight  of  such  a  generator  system,  and  is  also  the  largest 
component  in  terms  of  physical  dimension,  which  is  cumbersome  for 
launch  and  also  decreases  the  reliability  of  the  system  due  to  possible 
puncture  by  meteoritic  dust.  For  this  reason,  it  is  desirable  to  minimize 
the  radiator  area.  For  a  given  output  power  and  maximum  cycle  temper¬ 
ature  T^,  the  radiator  area  depends  upon  the  reject  temperature.  If  the 
rej««t  temperature  is  set  too  close  to  the  heat  input  temperature  T^,  then 
the  cycle  overall  thermal  efficiency  is  too  low,  and  the  ratio  of  heat 
rejected  to  useful  power  is  large.  Since  the  radiator  area  is  proportional 
to  the  total  heat  which  must  be  rejected,  the  radiator  area  becomes  large. 
On  the  other  hand,  if  the  reject  temperature  is  set  too  low,  then  the  cycle 
efficiency  may  be  large,  but  since  the  heat  rejection  is  proportional  to  the 
fourth  power  of  the  rejection  temperature,  the  radiator  area  becomes 
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large  again.  Thus,  there  is  an  optimum  rejection  temperature  for  which 
the  radiator  area  is  a  minimum.  We  will  illustrate  this  optimum  by  con¬ 
sidering  a  simplified  Ranklne  (vapor)  cycle. 

Rankine  Cycle 

The  Rankine  cycle  is  extremely  attractive  for  space  applications 
because  the  heat  rejection  is  at  constant  temperature,  and  the  pumping 
power  is  small,  see  Figure  14.40.  The  top  temperature  in  the  cycle  is 
T^;  the  generator  operates  between  0  and  1;  the  working  fluid  is  then  con¬ 
densed  from  1  to  2  during  which  heat  per  unit  mass  of  working  fluid 
is  rejected  from  the  space  radiator.  The  liquid  is  then  pumped  from  2  to 
3,  and  heat  is  added  by^the  nuclear  reactor  from  3  to  0.  It  is  assumed 
that  there  is  no  temperature  drop  through  the  walls  of  the  radiator  so  that 
the  radiator  surface  is  also  at  temperature  T^^.  It  is  also  assumed  that 
there  is  no  radiation  Interaction  between  parts  of  radiator,  or  with  the 

spacecraft,  and  that  the  friction  pressure  drop  is  negligible. 

« 

Now  the  heat  rejected  is  equal  to  the  heat  radiated,  so  that 
•  4 

m  Q  .  =  a  e  AT,  (a) 

rej  1 

where  m  is  the  mass  flow  rate  of  working  fluid,  a  is  the  Stefan -Boltzmann 
constant,  e  is  the  total  emissivity  of  the  radiator  surface  (assiuned  con¬ 
stant),  and  A  is  the  total  radiator  area.  From  the  first  law  of  thermody¬ 
namics,  the  total  power  P  from  the  generator  is  given  by: 

P  =  m  (Q^  -  (b) 

and  the  overall  thermal  efficiency,  77,  is  given  by: 


V  = 


Q  -  Q  . 

^in  ^rej 


Q 


in 


(c) 


17Z 
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*  ‘ 


Figure  14.40.  Rankine  Cycle 


175 


From  (b)  and  (c),  m  may  be  eliminated  to  yield: 


m  Q 


rej 


=  P(1  -  JZl 
r) 


(d) 


Now  the  efficiency  for  a  Rankine  cycle  may  be  expressed  as  a  fraction  7J 

c 

of  the  Carnot  efficiency  for  a  Carnot  cycle  operating  between  the  same  inlet 
and  rejection  temperatures,  as  follows: 


Substitution  of  (d,  e)  into  (a)  yields: 


A* 


aeT  ^  A 
_ o 

P 


1  -..n 
n 


(14.233) 


where  A*  is  a  non-dimensional  area.  For  given  values  of  c,  T^,  and  the 
power  level  P,  the  radiator  area  is  then  proportional  to  A*,  which  is 
given  by  the  right  side  of  (14.233).  It  is  seen  that  as  the  overall  cycle 
efficiency  becomes  either  zero  or  T]^,  A*  becomes  infinitely  large.  To 
determine  the  minimum  value  of  A*,  (14.233)  is  differentiated  and  set 
equal  to  zero;  the  efficiency  which  then  minimizes  the  radiator  area  is 
given  by: 


5  -  J25  -  16  Tj 

rj=  - ^ ^  (14.234) 

and  the  optimum  reject  temperature  is  obtained  from  (e): 

T,  5  -  V  25  -  16  T? 

' - Ti, - ^ 


174 


while  the  optimum  value  of  A*  is  obtained  from  (14.233).  These  values 

are  shown  in  Figure  14.41  as  a  fimction  of  the  Carnot  efficiency.  The 

ratio  (T  /T  )  is  a  weak  function  of  77  ,  since  it  varies  from  0.8  for  77 
1  o  c  c 

=  0  to  0.75  for  77  =1.  On  the  other  hand,  A*  depends  strongly  on  77  . 

c  c 

For  typical  Rankine  cycles  for  space  applications,  77  ~  0.65,  77  ~  0.15, 

o 

and(T^/T^)Rj  0.77. 

It  is  also  possible  to  minimize  the  total  weight  of  the  entire  space 
power  system,  by  utilizing  the  specific  weight  of  the  radiator,  nuclear 
reactor,  and  generator;  such  a  stucfy  has  also  been  made  in  reference  35. 
However,  for  large  generator  system,  the  optimum  temperature  is  still 
between  0.75  and  0. 8.  At  the  present  state  of  development  of  refractory 
materials,  it  appears  that  the  upper  temperature  for  Rankine  cycles  may 
be  about  2500°R  with  a  corresponding  reject  temperature  of  1900°R. 

These  temperatures  require  the  use  of  a  liquid  metal,  such  as  potassium 
or  sodium  for  the  working  fluid,  buir  are  too  low  for 

thermal  ionization  even  with  the  addition  of  cesium,  so  that  some  form  of 
non-thermal  ionization  must  be  used.  Calculations  for  magnetically  in¬ 
duced  ionization,  however,  indicate  that  this  may  be  possible  in  potassium 
for  relatively  high  field  strengths,  but  low  Mach  numbers  and  total  temper¬ 
atures.  A  typical  calculation  is  shown  in  Figure  14.42.  Although  this 
effect  is  not  yet  confirmed  experimentally,  the  high  power  density  makes 
magnetohydrodynamic  generators  with  magnetically  induced  ionization 
extremely  attractive  for  space  power  applications. 

Brayton  (Gas)  Cycle 

The  Rankine  cycle  as  described  above  has  one  main  disadvantage:  to 
operate  at  high  temperatures ,  which  is  necessary  to  minimize  the  radiator 
size,  liquid  metals  must  be  used  as  the  working  fluid  which  are  quite  corro¬ 
sive.  On  the  other  hand,  an  inert  gas  such  as  helium,  which  could  be  used 
in  a  gas  cycle,  is  not  corrosive.  Although  the  gas  cycle  has  two  disadvan¬ 
tages  for  space  applications;  namely,  the  heat  rejection  is  at  a  continuously 
varying  temperature,  and  a  large  compressor  is  required;  nevertheless, 
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Figure  14.  42,  Theoretical  Calculation  of  Conductivity  for  Magnetically  Induced 

Ionization  in  Potassium  Vapor  (F.  Shair,  personal  communication) 


it  is  of  interest  to  calculate  the  radiator  area  for  such  a  cycle,  and  com¬ 
pare  it  to  the  Rankine  cycle.  The  comparison  can  be  made  on  two  different 
bases:  for  the  same  top  temperature  in  the  cycle,  the  radiator  areas 
may  be  compared;  or  for  the  same  radiator  area,  the  top  temperatures 
may  be  compared.  For  the  gas  cycles  certain  idealizations  will  be  made, 
similar  to  those  previously  made  in  the  Rankine  cycle;  namely  zero  tem¬ 
perature  drop  across  the  walls  of  the  radiator,  and  zero  friction  jpressxire 
drop  which  is  valid  for  high  pressure  operation.  The  cycle  is  shown  in 
Figure  14.43;  0-3'  -  4-5'  is  the  ideal  gas  cycle,  for  isentropic  compressicai 
(4-5')  and  generation  (0-3').  More  realistically,  the  compressor  and  gen¬ 
erator  are  not  isentropic  and  have  an  entropy  increase,  shown  as  4-5  and 
0-3,  respectively. 


Since  the  temperature  of  the  working  fluid  in  the  radiator  is  constantly 


changing,  we  will  write  the  radiator  heat  balance  for  a  small  amount  of 

heat  rejection  dQ  .  which  occurs  between  T  -  7  AT  and  T  +  ^  AT.  The 
rej  2  2 

differential  radiator  area  dA  required  is  therefore  given  by: 


m 


do  .  =  aeT^dA 

^ej 


(f) 


At  constant  pressure,  dQ  =  C  dT;  then  integration  of  (f)  yields: 

P 


(g) 


where  the  specific  heat  at  constant  pressure,  C^,  has  been  taken  constant. 
Now  the  total  heat  rejected  per  unit  mass  of  gas  is  given  by: 


Q 


rej 


Cp  (T3  -  T^) 


(h) 
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Figure  14.  43.  Gas  Cycle 
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Use  of  (h,  d)  in  (g)  yields: 


A*  = 


P 


(14.237) 


where  rj  is  the  overall  thermal  efficiency  of  the  gas  cycle,  given  by: 


= 


Q.  -Q  . 

in  rei 


Q 


in 


(i) 


If  we  let  the  compressor  and  generator  adiabatic  efficiencies  be  defined  as: 


T  I  _  T 
5  4 


wViere  I'o/Tj  =  T5  /  j 
then  it  is  easily  shown  that 


(j) 


(k) 


(1) 
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so  that  (14.237)  becomes: 


From  (14. 238)  it  is  seen  that  for  a  given  value  of  T^,  the  radiator  area 
depends  on  two  temperatures,  T  and  T  ,  in  contrast  to  the  Rankine  cycle, 

iS  4 

in  which  the  area  depended  only  on  the  condenser  temperature.  The  deter¬ 
mination  of  the  minimum  radiator  area  is  therefore  more  complicated. 

The  usual  procedure  used  to  determine  the  minimum  radiator  area  is  to 

simply  compute  the  area  for  various  values  of  (T  /T  )  and  (T  /T  )  and 

36  ^  °  ^  ^ 

locate  the  minimum  in  this  manner.  The  results  of  such  a  calculation 
are  shown  in  Table  14.5,  for  two  different  compressor  and  generator 
efficiencies.  Note  that  the  non-dimensional  radiator  area.  A*,  is  very 
sensitive  to  these  efficiencies,  and  is  much  larger  than  the  corresponding 
values  for  a  Rankine  cycle.  Note  also  that  the  values  of  (T  /T  )  are  0.75, 

O  O 

which  is  similar  to  the  result  for  an  ideal  Rankine  cycle,  while  the  heat 
rejection  starts  at  T  and  decreases  to  T  .  It  is  this  large  change  in  heat 

o  4 

rejection  temperature  which  is  responsible  for  the  large  radiator  area. 

The  low  values  of  T^/T^  are  required  to  obtain  net  power  output;  for  an 

ideal  gas  cycle,  then  the  minimization  yields  T  ~  T  ;  but  when  the  com- 

4 

ponent  efficiencies  are  less  than  unity,  it  is  necessary  that  T  <  T  in 

4  o 

order  to  obtain  net  power. 
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TABLE  14.5.  BRAYTON  CYCLE  FOR  SPACE  POWER 


(a) 

(b) 

77  =77  =0.8 
'c  'g 

7J  =77  =  0. 85 
'c  'g 

A** 

145 

84 

V 

0.15 

0.15 

T„/T 

3  o 

0.75 

0.75 

T/To 

0.3 

0.375 

T4/T3 

0.4 

0.5 

TJT. 

5  4 

1.57 

1.50 

0.69 

0.70 

2.55 

2.4 

COMPARISON  TO  RANKINE  CYCLE 
(T^  =  2460°  R,  A*=  16)  FOR  SAME  RADIATOR  AREA 


T 

4260°  R 

3730°  R 

0 

^3 

3190°  R 

2790°  R 

'^4 

1280°  R 

1395°  R 

T. 

2000°  R 

2090°  R 
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Column  (a)  was  calculated  by  the  minimization  process;  the  values 

in  column  (b)  were  obtained  by  minimizing  T  /T  while  keeping  the  cycle 

o  o 

efficiency  the  same.  Note  that  for  column  (a),  the  overall  cycle  efficiency 
is  15%,  which  is  close  to  the  Rankine  cycle,  which  implies  that  for  a 
given  power  level,  the  thermal  input  to  both  cycles  will  be  about  the  same. 


Although  the  radiator  area  for  the  gas  cycle  is  much  larger  than  the 
Rankine  cycle  for  the  same  cycle  maximum  temperature ,  because  the 
gas  cycle  can  use  inert  gases  and  avoid  the  corrosion  problem,  it  may  be 
possible  to  operate  the  gas  cycle  at  a  higher  temperature.  Now,  for  the 
same  radiator  area,  (14.237)  predicts  that  the  top  temperature  will 
vary  as  the  one-fourth  power  of  A*.  Such  a  comparison  is  made  at  the 

bottom  of  Table  14.5,  where  it  may  be  seen  that  for  n  =  n  =0.85,  the 

o  ^  ® 

top  temperature  is  3730  R.  Although  this  temperature  is  quite  high,  it 

may  be  possible  to  develop  nuclear  reactors  for  use  with  inert  gases  in 

this  range.  This  temperature  is  probably  too  high  for  the  use 

of  gas  turbines,  but  not  for  MHD  generators. 


A  major  problem  with  the  gas  cycle  is  that  a  gas  compressor  is 
required.  Since  the  gas  pressure  in  the  reactor  is  likely  to  be  high  in 
order  to  minimize  the  heat  exchanger  surfaces,  it  is  unlikely  that  a 
magnetohydrod3Tiamic  compressor  can  be  utilized,  and  an  electrically 
driven  rotating  compressor  is  required.  If,  for  the  example  in  column 
(a),  the  electric  motor  for  the  compressor  has  an  efficiency  of  74%,  then 
all  of  the  electrical  output  would  be  used  to  drive  the  compressor,  and 
there  would  be  zero  useful  power  output. 


Of  course,  a  magnetic  field  must  be  incorporated  into  a  nuclear- 
MHD  space  power  system,  which  raises  two  additional  problems,  namely, 
the  weight  of  the  magnet,  and  the  rejection  of  heat  from  the  magnet.  Now, 
most  space  applications  for  use  with  electrical  propulsion  will  require 
from  one  to  twenty  megawatts  of  power;  from  Section  14. 81,  the  weight 
of  a  copper  wound  air  core  magnet  is  prohibitively  large,  and  the  weight 
of  an  iron  core  magnet  is  only  slightly  less.  In  addition,  the  heat 
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generated  In  the  windings  must  be  removed  by  an  additional  radiator. 

For  example,  if  one  allows  10%  of  the  total  power  to  energize  the  field 
coils,  then  the  entire  power  plant  including  its  radiator  must  be  larger 
by  11%.  The  main  radiator  rejects  85%  of  the  thermal  energy  (if  the 
overall  thermal  efficiency  is  15%),  and  if  the  windings  are  operated  at 
the  same  temperature  as  the  main  radiator,  then  the  additional  radiator 
load  is  15%.  However,  since  the  electrical  resistivitjr  of  most  metals 
increases  with  temperature,  and  the  size  of  the  windings  increases  with 
resistivity  when  the  fraction  of  power  consumed  by  the  windings  is  held 
constant,  it  may  be  desired  to  maintain  the  windings  at  some  lower  tem¬ 
perature,  and  use  a  heat  pump  to  increase  the  rejection  temperature. 

This  reduces  the  size  of  the  radiator  for  the  windings,  but  additional 
power  is  required  for  the  heat  pump. 

The  problems  of  the  weight  of  the  windings  and  the  rejection  of  the 
ohmic  heating  in  the  windings  may  be  solved  by  the  use  of  hi^  field  super¬ 
conducting  materials.  The  hi^  electric  current  densities  of  supercon¬ 
ductors  reduce  the  weight  of  the  windings,  while  the  only  required  heat 
rejection  is  that  which  is  conducted  into  the  magnet  dewar,  which  tends 
to  boil  off  the  liquid  helium  in  which  superconductors  are  kept.  The 
liquid  helium  can  either  be  recooled  by  a  cryostat  which  requires  addi¬ 
tional  weight,  power,  and  radiators;  or  sufficient  liquid  helium  can  be 
stored  to  allow  boil -off  of  the  helium  during  the  required  mission  time. 

14. 8d  Other  Gas  Cycles 
Intercooled 

Since  the  radiator  area,  per  unit  heat  rejection,  decreases  with  in  creasing 
heat  rejection  temperature,  it  appears  advantageous  to  start  compression  when 
the  working  fluid  reaches  T_ ,  with  the  compression  divided  into  many  stages 
with  intercooUng  by  radiation  between  each  compression  stage.  Ideally,  this 
may  be  represented  as  an  infinite  number  of  infinitesimal  compressions,  and 
is  a  constant  temperature  line  on  the  temperature  entropy  diagram,  see 
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Figure  14.44 .  it  is  easily  shown  that  if  each  compression  stage  has  an  adiabatic 

compression  efficiency  77  ,  then  the  rejected  heat  from  the  intercooler  radiator 

c 

is  given  by: 
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From  the  definition  of  the  generator  adiabatic  efficiency, 
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the  pressure  ratio  can  be  obtained: 
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where  it  has  been  assumed  that  the  pressure  drops  are  negligible.  Also 


Q  =  C  (T  -  T_);  and  Q'  .  =  C  (T^  -  T_). 
in  p'  o  5  re]  p'  3  5' 


The  overall  efficiency  is 
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or,  with  the  use  of  (m,  n,  q): 


In  a  similar  manner,  the  radiator  area  can  be  expressed  as: 


In  order  to  effect  a  comparison  with  the  Brayton  cycle ,  calculations  have  been 

made  with  (14.239,  14.240),  assuming  that  t)  =  7)  =0.8,  and  77  =  0.15.  Then 

c  g 

(14.239)  gives  a  relation  between  (T  /T  )  and  (T  /T  ),  and  A*  may  be  calculated 

o  A  4  0 

for  various  values  of  T„/T  .  The  minimum  radiator  area,  for  these  assump- 

3  o 

tions  is  A*  =  220  as  shown  in  Table  14.6,  for  which  (T  /T  )  w  0.7,  and 

o  O 

(T  /T  )  ~  0.49.  Compared  to  the  Brayton  cycle  (a),  the  radiator  area  has  in- 
4  «3 

creased.  This  may  be  attributed  to  additional  heat  rejection  due  to  intercooling 

which  requires  that  T  ,  and  hence  T  ,  be  decreased  to  obtain  the  same  efficiency. 

r 

This  is  the  opposite  effect  of  the  usual  effect  of  intercooling  for  gound  based  gas 

A 

cycles,  where  the  lowest  temperature  in  the  cycle  corresponds  to  the  available 
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cooling  fluid;  intercooling  at  the  lowest  temperature  in  the  cycle  decreases  the 
net  heat  rejection. 

To  investigate  whether  complete  intercooling,  as  shown  in  Figure  14.44, 
decreases  the  radiator  area,  a  calculation  was  performed  for  f]  -  0.90. 

The  optimum  temperature  ratio  is  (T^/T^)  =  0.5,  but  the  radiator  area  is  still 
very  large  (A*  =  144)  as  shown  in  Table  14.6,  as  compared  to  a  simple  Brayton 
cycle  for  lower  component  efficiencies. 

Another  possible  modification  is  the  use  of  regeneration,  see  Figure  14.44. 

Again,  the  radiator  area  depends  on  ,  77  ,  (T  /T  ),  (T  /T  ),  and  the  temi)er- 

eg  u  O  4  0 

ature  drop  across  the  walls  of  the  regenerator.  Note  that  T  and  T  are  fixed 

5  6 

by  the  other  parameters.  Calculations  for  minimizing  the  radiator  area  show 

37 

that  a  value  of  A*  =  48  may  be  attainable,  which  further  reduces  the  required 
top  temperature  in  order  to  achieve  the  same  radiator  area  as  the  Rankine  cycle, 
see  Table  14.6.  However,  the  pressure  ratio  is  small,  p  /p  «  1.75,  which 

o  O 

simplifies  the  compressor  problem,  but  which  then  requires  that  the  pressure 
drop  in  the  regenerator  be  made  very  small.  Now,  the  pressure  drop  in  the 
regenerator  is  inversely  proportional  to  the  product  of  the  absolute  pressure  in 
the  regenerator,  the-number  of  tubes,  and  the  tube  cross-sectional  area.  Thus, 
to  keep  the  regenerator  small,  high  pressure  operation  of  the  cycle  is  required. 
On  the  other  hand,  in  order  to  achieve  non-thermal  ionization  in  the  MHD  gener¬ 
ator,  low  static  pressures  are  required,  which  then  requires  supersonic  flow  in 

the  generator  with  a  supersonic  dlffusor  downstream  of  the  generator.  Calcu- 

37 

lations  of  the  losses  associated  with  such  a  diffusor  indicate  that  there  is  an 
additional  decrease  of  the  generator  adiabatic  efficiency  of  about  10%  which  is 
undesirable  since  it  is  necessary  in  a  gas  cycle  that  the  generator  (and  com¬ 
pressor)  efficiency  be  as  high  as  possible. 

14. 8e  Two  Phase  Cycles 

Of  course,  it  is  not  necessary  for  the  working  fluid  o4a  MHD  generator  to 
be  an  ionized  gas;  it  could  for  example  be  a  liquid  metal,  or  a  mixture  of  gas  or 
vapor  and  liquid  metal.  Such  a  generator  has  the  advantage  of  the  very  high 
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TABLE  14.6.  RADIATOR  AREA  PARAMETER  FOR  VARIOUS  CYCLES 
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REGENERATIVE  BRATTON 


Figure  14.  44.  Simple  Vapor  and  Gas  Cycles  for  Space  Applications 
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electrical  conductivity  of  liquid  metals  at  all  temperatures.  The  main  problem 

is.the  conversion  of  the  thermal  energy  from  the  nuclear  reactor  into  kinetic 

energy  of  motion  of  the  liquid  metal;  normally  a  gas  or  vapor  is  used  for  this 

purpose.  This  can  be  accomplished  by  mixing  a  hot  liquid  metal  of  low  vapor 

k  I 

pressure  with  a  cooler  liquid  metal  of  high  vapor  pressure.  The  hot  liquid 
metal  causes  the  lower  temperature  liquid  metal  to  vaporize.  The  mixture  is 
passed  into  a  nozzle  which  allows  the  vapor  to  expand  and  accelerate  the  droplets 
of  liquid  metal.  The  mixture  is  then  separated,  without  appreciable  loss  of 
kinetic  energy  or  head  of  the  liquid  metal,  and  this  head  then  forces  the  liquid 

metal  through  the  generator  against  the  Lorentz  force.  One  such  possible 

38  '3 

cycle  is  shown  in  Figure  14.45.  The  lower  loop  for  the  liquid  metal  while  the 

A 

upper  loop  is  for  the  liquid  metal  vapor.  The  liquid  metal  enters  the  reactor  at 
9,  leaves  at  10,  and  is  mixed  with  the  high  vapor-pressure  liquid  metal  in  the 
"mixer,"  in  which  some  or  all  of  the  latter  liquid  metal  is  vaporized,  and  ther¬ 
modynamic  equilibrium  between  the  two  streams  is  reached  at  1.  The  entire 
mixture  then  passes  through  the  nozzle,  where  the  vapor  and  liquid  metal  droplets 
are  accelerated,  and  enters  the  separator  at  2.  The  separated  liquid  metal  uses 
the  acquired  kinetic  energy  to  flow  through  the  generator  and  reactor.  The 
separated  vapor  passes  through  a  regenerator  (heat  exchanger),  radiator  where 
it  is  condensed,  a  pump,  and  is  then  reheated  by  the  regenerator  before  entering 
the  mixer.  Calculations  of  the  efficiency  of  this  cycle  indicate  about  11%  is 

achievable  if  the  mass  flow  rate  through  tiie  lower  loop  is  100  times  the  mass 

38 

flow  rate  through  the  upper  loop;  this  means  that  the  Increase  in  temperature 
in  the  reactor,  T^^^  -  T^,  will  be  small.  The  high  mass  flow  ratio  of  liquid 
metal  is  required  to  keep  the  velocity  of  the  mixture  small  in  the  separator  since 
frictional  losses  increase  approximately  as  the  square  of  the  velocity.  The  high 
mass  flow  ratio  also  precludes  the  use  of  a  single  fluid,  since  a  single  fluid  if 
expanded  from  the  saturated  liquid  line  to  a  temperature  equal  to  about  3/4  of 
the  stagnation  temperature  will  have  a  large  mixture  velocity,  and  hence,  an 
excessive  friction  loss  in  the  separator. 
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Figure  14.  45.  Two-Fluid  MHD  Power  Conversion  Cycle 


Besides  the  problem  of  friction  losses  in  the  separator,  because  of  the 
small  temperature  change  in  the  reactor,  the  energy  extraction  per  unit  mass 
flow  in  the  generator  is  quite  small,  which  results  in  a  generator  of  extremely 
small  aspect  ratio,  see  Figure  14. 46,  and  hence  large  end  losses  in  the  absence 
of  inlet  and  outlet  vanes;  But  if  inlet  and  outlet  vanes  are  used,  these  will 
increase  the  friction  pressure  drop  further.  It  is  possible  that  there  are  other 
geometrical  arrangements  for  the  nozzle,  separator,  and  generator  which  will 
have  a  smaller  wetted  surface  area,  and  a  more  favorable  generator  aspect 
ratio. 

14. 8f  Heat  Rejection  to  Propellant 

One  possible  use  for  space  electrical  power  is  electrical  propulsion.  This 
raises  the  possibility  of  using  the  reject  heat  from  the  space  electrical  power 
system  to  preheat  the  propellant,  and  thereby  avoid  the  use  of  a  radiator.  For 
a  regenerative  cycle,  the  preheat  heat  balance  is  (see  Figure  14.44): 

m  Q  ,  =  P  0^^)  =  m  C  (T^  -  T.)  (r) 

rej  \  77  /  p  p  6  i 

P 

where  T  is  the  initial  temperature  of  the  propellant,  and  m  and  is  the  mass 
i  P  ^p 

flow  rate  and  specific  heat  of  the  propellant,  assumed  constant.  Also,  it  has 
been  assumed  for  simplicity  that  there  is  a  negligible  temperature  drop  between 
the  working  fluid  in  the  space  power  cycle  and  the  propellant  in  the  heat  exchanger. 
Now,  the  electrical  power  from  the  cycle  is  used  to  heat  the  propellant  further: 
this  may  be  in  a  plasma  jet,  resistance  heater  or  any  other  device;  we  assume 
that  the  efficiency  for  this  process  is  100%.  Thus: 

P  =  m  C  (T  -tJ  (s) 

p  p  o  6 
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Figure  14.46,  Nozzle,  Separator,  and  Generator  Configuration  for  500-kwe 
Conversion  System 


where  Tq  is  the  final  stagnation  temperature  of  the  propellant.  Also,  for  a 
regenerative  cycle, 


V  = 


1  - 


1  - 


(t) 


since  for  a  perfect  gas  cycle  (T  /T  )  =  (T  /T  ).  Division  of  (r)  by  (s),  and  use 

O  o  D  D 

of  (t)  yields: 


T  -T  = 
o  i 
P 


(u) 


Now,  to  maximize  the  propellant  stagnation  temperature,  must  be  made  as 
small  as  possible,  and  (u)  becomes  ~  T^.  Thus,  the  stagnation  temperature 
of  the  propellant  is  no  higher  than  if  the  propellant  had  been  directly  heated  by 
the  reactor!  For  fuel  rod  reactors,  this  cycle  therefore  has  little  advantage. 
However,  if  the  reactor  were  a  gaseous  nuclear  reactor,  which  virtually  pre  ¬ 
cludes  mixing  of  the  propellant  with  the  nuclear  fuel,  then  this  cycle  has  the 
advantage  that  the  fuel  and  propellant  are  kept  separate,  and  also,  much  higher 
temperatures  may  be  achievable  in  a  gaseous  nuclear  reactor  than  in  a  reactor 
which  uses  solid  or  liquid  nuclear  fuel. 

For  fuel-rod  cycles,  this  cycle  may  be  improved  by  sending  the  propellant, 

39 

after  preheat  but  before  electrical  heating,  through  the  reactor.  Then  (s)  is 
replaced  by 


P 


m  C 
P  P, 


(V) 


Following  the  same  procedure  as  in  the  previous  case  (with  T^m  0), 


T 

o 

P 


(w) 


194 


Thus,  the  stagnation  temperature  of  the  propellant  is  almost  twice  that  which 
would  have  been  obtained  if  the  propellant  were  directly  heated  by  the  nuclear 
reactor. 

As  previously  mentioned,  gaseous  nuclear  reactors  may  be  capable  of 
higher  temperatures  than  solid  fuel  rod  reactors,  but  direct  heating  of  the  pro¬ 
pellant  is  difficult,  since  the  gaseous  propellant  must  be  separated  from  the 
gaseous  nuclear  fuel  prior  to  expulsion.  Although  vortex  separation  has  been 

40 

suggested,  it  may  also  be  possible  to  use  an  MHD  generator  for  this  purpose. 
The  working  fluid  consists  of  a  mixture  of  propellant  and  gaseous  nuclear  fuel. 
The  nuclear  reaction  heats  the  entire  mixture  to  a  sufficiently  high  temperature 
to  ionize  the  gas ,  and  the  mixture  flows  through  an  MHD  generator  in  which  the 
nuclear  reaction  is  stopped  and  enthalpy  is  extracted  from  the  mixture.  The 
temperature  of  the  mixture  is  thereby  reduced  imtil  the  nuclear  fuel  condenses, 
which  is  then  extracted  from  the  flow  and  recirculated  to  the  reactor.  The 
electrical  power  which  has  been  generated  in  the  MHD  generator  is  then  added 
to  the  propellant,  which  is  now  free  of  nuclear  fuel,  and  exhausted  from  the 
electrothermal  engine.  Such  a  system  is  interesting  but  highly  speculative  at 
present. 
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14.9  MHD  GENERATOR  EXPERIMENTS 

2 

The  first  MHD  generator  experiments  were  those  of  Halasz  and  Karlovitz, 
which  were  carried  out  at  Westinghouse  up  to  about  1946,  using  electron -beam 
ionization.  Due  to  rapid  recombinaticm  and  formation  of  negative  icms,  the 
generated  power  density  was  far  below  theoretical.  No  further  experiments 
were  performed  until  about  1958,  when  experiments  were  initiated  using  thermal 
ionization  of  seeded  gases,  see  Table  14.7.  The  earliest  of  these  experiments 
utilized  plasma  jets  as  the  method  for  thermal  heating  of  the  gas,  since  gas 
temperatures  could  easily  be  obtained  which  are  much  higher  than  those  avail¬ 
able  from  combustion  flames. In  these  experiments,  the  measured  electrical 
conductivity  and  generated  power  agreed  within  a  factor  of  two  of  the  theoretically 
predicted  performance,  although  the  total  power  was  small  (less  than  12KW)  and 
the  run  times  were  less  than  10  seconds.  These  experiments  mainly  verified 
the  validity  of  the  physical  principle  of  MHD  power  generation.  A  typical  voltage - 
current  plot  is  shown  in  Figure  14.47. 

42 

The  first  combustion  experiment  was  performed  at  Westinghouse  Electric, 

in  which  over  10  KW  was  generated;  followed  rapidly  by  a  similar  experiment  At 

MHD  Research,  and  a  much  smaller  experiment  at  the  General  Electric 

44 

Research  Laboratory.  These  experiments  also  verified  the  theory  closely; 

however,  the  magnetic  interaction  parameter  for  these  experiments  was  quite 

small.  The  largest  combustion  experiment  at  the  present  time  is  the  AVCO 
45 

Mark  II,  which  generated  1300  kilowatts  in  October,  1962.  A  Hall  generator 

46 

has  also  been  successfully  operated  on  an  air -hydrogen  flame,  with  also 

excellent  agreement  with  theory.  Combustion  generation  experiments  have  also 

47  48 

been  performed  in  Poland,  and  have  been  reported  in  the  Soviet  Union,  but 

details  of  the  latter  are  not  yet  available. 

Power  generation  experiments  have  not  yet  been  performed  using  non- 

thermal  ionization  as  of  October,  1962,  but  experimental  equipment  is  under 

49 

construction  at  C.  A.  Parsons  Nuclear  Research  Centre  and  the  General 

50 

Electric  Space  Sciences  Laboratory. 
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Genera!  Electric  General  Electric  General  Electric  .  , 

AVCO  o  o  •  DU  DU  Westinghouse 

- -  Space  Sciences  Research  Research  _ _ _ 

MHD  Research  Mark  I _ Mark  n  Laboratory  Laboratory  Laboratory  Model  1  Model  2 
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14.9a  mhD  Generator  Materials 


For  the  plasma  jet  and  combustion  experiments,  some  initial  experiments 
were  performed  using  stabilized  zirconia  for  the  insulator;  however,  at  high 
temperatures  the  zirconia  generally  becomes  electrically  conducting  and  most 
experimenters  have  changed  to  magnesia  for  the  insulator.  One  unexpected 

46 

effect  is  that  the  resistivlly  of  the  ceramic  decreases  durii^  an  experiment; 
it  is  not  yet  known  whether  this  is  due  to  the  presence  of  potassium  which  could 
change  the  chemical  composition  of  the  surface  of  the  insulator.  The  presence 
of  molten  ash  will  probably  Intensify  this  problem. 

For  very  large  MHD  generators,  wall  cooling  may  be  possible  without 
adversely  affecting  the  performance,  and  experiments  have  been  performed  in 
which  water-cooled  tubes  krere  used  as  the  insulator.  These  same  tubes  can 
also  serve  as  the  electrodes  if  the  tubes  follow  equipotential  lines  as  in  the 
series -segmented  geometry,  and  if  an  insulation  is  placed  in  between  the  tubes. 
Since  the  temperature  of  the  tubes  is  low,  presumably  the  temperature  of  the 
insulation  between  them  will  also  be  low  and  hence  alumina  may  be  satisfactory. 

3  41 

The  most  commonly  used  electrode  material  has  been  graphite,  ’ 

46 

although  silicon  carbide  has  also  been  utilized.  However,  zircon  compounds 

have  also  been  utilized  because  of  their  relatively  high  electrical  conductivity, 

43 

both  with  and  without  embedded  tungsten.  In  general,  if  the  electrodes  are 

sufficiently  hot  to  emit  electrons  thermionically,  there  has  been  almost  zero 
41 

cathode  voltage  drop;  otherwise  voltage  drops  up  to  100  volts  have  been 
observed.  The  gra^Mte  and  carbide  electrodes  have  been  observed  to  erode 
very  rapidly  which  is  probably  due  to  chemical  reactions  with  water  vapor  in 
the  products  of  combustion.  For  a  continuously  operating  generator,  these 
types  of  electrodes  could  be  continuously  fed  into  the  generator,  which  may 
cause  a  leakage  and  seal  problem  if  the  gas  inside  the  generator  is  above  atmos- 
I^eric  pressure.  The  pressure  of  molten  ash  could  interfere  with  proper  oper¬ 
ation  of  the  electrodes  as  well  as  the  insulator. 


\ 
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14.9b  Seed  Material 


The  seed  has  generally  been  injected  as  a  salt,  K  CO  ,  or  as  a  hydroxide, 

it  3 

KOH.  The  salt  dissolves  readily  in  water  which  provides  a  convenient  method 

for  injection,  but  powdered  potassium  carbonate  has  also  been  used.  The 

hydroxide  dissolves  easily  in  alcohol,  which  then  can  be  mixed  directly  with  the 
43 

fuel,  eliminating  the  necessity  for  a  separate  injection  system  and  also  allow¬ 
ing  for  more  precise  control  of  the  seed  Injection  rate. 

For  closed  cycle  operation,  pure  alkali  metal  is  required  for  the  seed. 
This  can  be  provided  by  injection  of  the  liquid  metal;  use  of  a  side  stream  which 

9 

is  bubbled  through  the  liquid  metal;  a  spray  evaporator,  or  a  direct  boiler  for 
the  liquid  metal. 
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Magnetohydrodynamic  power  generaticm  has  now  been  under  activ 
development  for  over  four  years,  but  there  has  not  yet  appeared 
any  complete  description  of  the  theory.  This  report  is  intended 
to  close  this  obvious  gap.  Most  of  the  theory  presented  herein 
was  developed  by  the  author  and  personnel  at  the  General 
Electric  Company.  Some  of  this  material  has  not  been  published 
previously;  those  results  which  have  been  published  are 
referenced. 

The  topics  covered  are:  electrical  conductivity  in  MHD 
generators,  optimum  "seed"  ratio,  local  analyses  of  the 
continuous  and  segmented  electrode  geometries;  Hall  geometry, 
helical  flow  geometry;  magnetically  induced  ionization; 
polytropic  efficiencies;  compressible  analyses  of  the  ccmstant 
velocity,  temperature,  Mach  number,  pressure  and  cross- 
sectional  area  flows;  end  losses;  AC  generation;  cycle 
efficiencies ;  and  a  summary  of  generating  experiments  at  the 
General  Electric  Company  and  other  places.  Geometries  other 
than  linear  are  not  considered  herein;  the  most  important  of 
those  omitted  is  the  vortex  generator. 
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